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Abstract 


A  measurement  of  G^/G^  at  values  of  0.38  and  0.50  GeF^  was 
performed  at  Bates  Linear  Accelerator  Center  using  a  focal  plane  polarimeter 
(FPP)  to  measure  the  recoil  proton  polarization  following  elastic  scattering  of 
polarized  electrons  on  a  hydrogen  target.  These  are  the  first  measurements 
of  all  three  polarization  components  and  determination  of  the  proton’s  form 
factors  via  spin  transfer  rather  than  Rosenbluth  separation.  The  ratio  G^  / G^ 
is  obtained  in  a  nearly  systematic  free  manner.  The  measured  polarizations 
are  consistent  with  the  Plane  Wave  Impulse  Approximation  (PWIA).  This  is 
the  first  experiment  to  use  the  newly-built  BATES-FPP  and  that  device  is 
described  in  detail. 
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Chapter  1 


Introduction 


This  dissertation  describes  an  experiment  (88-21)  performed  at  the  Bates 
Lineeur  Accelerator  Center  in  Middleton,  Massachusetts  to  measure  proton 
electromagnetic  form  factors.  While  these  have  been  measured  many  times 
before  using  Rosenbluth  separation,  limitations  of  that  method  inspired  us 
to  use  an  alternative  one  measuring  polarization  transfers  in  the  ^H(e,e'p) 
reaction.  This  method  allows  the  ratio  to  be  obtained  in  a  nearly 

systematic-free  manner.  The  polarization  of  the  outgoing  protons  was  measured 
using  a  focal  plane  polarimeter  (FPP)  built  by  our  group  for  this  and  other 
experiments.  To  the  extent  that  G’^  and  G^  axe  known,  this  experiment  also 
provided  on-site  calibration  of  the  FPP.  The  experiment  was  performed  at  two 
values:  0.38  and  0.50  GeF^,  and  is  the  first  such  experiment  to  measure  all 
three  proton  polarization  components  and  determine  the  form  factors. 

The  remainder  of  this  chapter  will  discuss  the  previous  Rosenbluth  form 
factor  measurements  and  the  models  used  to  describe  the  form  factors.  The 
next  chapter  of  this  dissertation  will  develop  the  formalism  explaining  how  to 
measure  form  factors  using  polarization  transfer.  The  experimental  set-up  used 
for  this  experiment  will  be  described  in  chapter  3.  Chapter  4  overviews  the  data 
analysis  while  chapter  5  discusses  the  results.  Appendix  A  provides  detailed 
derivations  of  cross  section  relationships  used  in  this  document.  Appendix 
B  describes  the  FPP  calibration  run  performed  at  the  Indiana  University 
Cyclotron  Facility  (lUCF)  while  appendix  C  delves  into  alignment  aspects  of 
the  FPP. 
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1.1  Using  Electron  Scattering  to  Evaluate  and 

Electron  scattering  has  been  a  successful  tool  for  many  years  in  stud3ring 
the  structure  of  the  nucleus.  The  main  advantages  of  the  electromagnetic  probe 
are  the  following: 

1)  The  electromagnetic  interaction  is  well  understood  and  can  be  described 
quantitatively  very  well  due  to  the  precision  of  the  underlying  theory: 
quantum  electrodynamics  (QED).  With  electron  scattering,  the  cross 
section  can  be  directly  related  to  the  transition  matrix  elements  of  the 
charge  and  current  density  operators  of  the  target,  i.e.  the  nuclear 
electromagnetic  structure. 

2)  The  electromagnetic  interaction  is  weaJc,  on  the  order  of  the  fine  structure 

constant  a  =  This  is  part  of  the  reason  that  QED 

gives  such  precise  predictions,  since  the  interaction’s  weakness  allows 
perturbation  theory  to  be  used.  The  weakness  also  allows  the  whole 
target  to  be  probed,  whereas  the  strong-force  interactions  of  proton¬ 
scattering  tend  to  just  probe  the  nuclear  surface.  The  target  is  also 
probed  without  greatly  disturbing  it,  whereas  with  proton  scattering  it 
is  difficult  to  separate  the  scattering  mechanisms  from  target  structure 
effects.  Electron  scattering  can  usually  (for  1/Z  »  a,  where  Z  is  the 
nuclear  charge)  be  treated  in  the  first  Bom  approximation,  i.e.  as  the 
exchange  of  a  single  virtual  photon. 

3)  Electron  scattering  allows  one  to  vary  the  momentum  transfer  q  while 
keeping  the  energy  transfer  iv  fixed  (and  vice  versa).  This  is  what 
separates  electron  scattering  from  photon  scattering  since  real  photon 


T.  de  Forest  Jr.  and  J.D.  Walecka,  Adv.  Phys.  16,  1  (1966). 
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scattering  has  the  constraint  that  =  0.  Virtual  (or  “off- 

mass  shell”)  photons,  however,  have  q^  <  0.  Thus  one  is  allowed  to 
study  the  ^-dependence  of  the  transition  matrix  elements  and  use  Fourier 
transforms  to  determine  the  spatial  distributions  of  the  nucleus’s  charge 
and  current  densities. 

Figure  1.1  illustrates  the  one- photon  exchange  process  in  (e,e')  scattering 
off  a  nucleon.  The  electron  has  initial  energy  E  and  initial  momentum  k.  After 
scattering  through  an  angle  6e  its  final  energy  is  E'  and  its  final  momentum  is 
k'.  The  electron  interacts  with  a  nucleon  by  the  exchange  of  a  virtual  photon. 
The  electron’s  energy  transfer  is 

u  =  E-E'  (1.1) 

and  its  3-momentum  transfer  is 

q  =  k-k'  .  (1.2) 

The  four-momentum  transfer  is  expressed  as  9^  =  (a;,  g  );  its  squared  magnitude 
as 

-f  =m^  =  -Q^  ,  (1.3) 

which  can  be  approximated  in  the  extreme  relativisitic  limit  (mg  «  0)  as: 

gJ«-2fc-jb'«-4i;e^?eSin^(^e/2)  .  (1.4) 

Using  only  Lorentz  covariance,  gauge  invariance,  and  parity  conserva¬ 

tion;  the  doubly-differentiaJ  cross  section  for  unpolarized  (e,e')  (inclusive)  scat¬ 
tering  can  be  shown  to  depend  on  only  two  response  functions,  Rl  and  Rt, 
which  contain  all  the  target  structure  information: 

•  (1-5) 


d^a 

dildw 


M.  Gourdin,  Nuovo  Cim.  21,  1094  (1961). 
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Figure  1.1  Feynman  diagram  for  one-photon  exchange  in  (e,e')  . 


ctm  is  the  Mott  cross  section.  It  represents  the  scattering  off  a  spinless,  pointlike 
particle  and  is  given  by: 


cos^(fie/2) 
4E2  sin^(0e/2) 


(1.6) 


iZxr  and  are  the  nuclear  response  functions  to  longitudinal  and  transverse 
photons,  respectively.  For  the  specific  case  of  elastic  (e,e')  scattering  from  the 
nucleon  this  differential  cross  section  is  expressed  as; 


r  da  1 

■  da  1 

- 1 

Rosenbluth 

V 

_ 1 

J  718 


1  +  r 


+  rGUQ^)i 


— ^ - h  2  tan^  — ) 

1+r^  2^ 


(1.7) 


where  r  =  /AMp,  A  detailed  derivation  of  this  equation  can  be  found  in 
appendix  A.  The  first  term  on  the  right  side,  the  nonstructure  differential  cross 
section,  is  just  the  Mptt  cross  section  times  a  target  recoil  factor.  Ge  and 
Gm  are  called  form  factors  and  contain  the  additional  information  about  the 
finite  spatial  extent  of  the  nucleon’s  charge  and  magnetization  distributions 
respectively.  One  can  see  from  comparing  the  two  cross  section  equations  that 


M.N.  Rosenbluth,  Phys.  Rev.  70,  615  (1950). 
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Rl  relates  to  the  charge  distribution  and  Rt  relates  to  the  magnetization 
distribution.  The  Rosenbluth  formula  (1.7)  is  vaJid  for  both  protons  and 
neutrons,  but  with  different  form  factors  (denoted  as  Gg(jvf)  neutrons  and 
^E(M)  protons).  The  accurate  measurement  of  these  form  factors  has 
consistently  been  one  of  the  leading  efforts  of  experimental  nuclear  physics. 
Likewise,  any  baryon  structure  model  that  wishes  to  be  considered  successful 
must  be  able  to  predict  their  values.  Before  describing  the  existing  form  factor 
measurements  and  models  several  points  should  be  made  about  them. 

1.2  The  Ge  and  Gm  Form  Factors 


As  is  shown  in  appendix  A,  Ge  and  Gm  axe  closely  related  to  the 
nucleon’s  charge  distribution  and  magnetic  moment  distribution,  respectively. 
Hence,  as  goes  to  zero,  G^  1,  G^  0,  G^  —*  fip,  and  G^  fin-  It 
should  be  noted  that  Ge  aud  Gm  (the  Sachs  form  factors)  are  not  the  only 
possible  pair  of  form  factors.  Often  used  axe  the  Dirac  (Fi)  and  Pauli  (F2) 
form  factors.  The  two  sets  are  simply  related: 


Gf  =  Fi - — — Fs 

®  ^  4M2  ^ 

Gm  =  Fi  +  kFz  . 


As  goes  to  zero;  Ff  ^  1,  Fj*  — ^  1,  F2"  —*  1,  and  F”  —*■  0.  These  form 
factors  predate  the  Sachs  form  factors.  The  use  of  the  Sachs  form  factors  has 
a  more  readily  understood  physical  interpretation  (Fourier  transform  of  charge 
and  magnetic  moment  distributions)  and  also  avoids  having  a  cross-term  in  the 
Rosenbluth  equation.  W 


Even  without  a  cross-term,  the  Rosenbluth  equation  presents  a  problem 
with  respect  to  the  extraction  of  G^:  the  t  term  in  front  of  G\j.  This  causes  the 


W  L.N.  Hand,  D.G.  Miller,  and  R.  Wilson,  Rev.  Mod.  Phys.  S6,  335  (1963). 
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G%  term  to  contribute  less  to  the  cross  section,  making  its  extraction  harder, 
as  is  increased.  Compounding  the  problem  for  the  neutron  is  the  fact  that 
is  much  smaller  (since  it’s  chargeless)  than  making  its  extraction  even 
more  difficult. 

A  final  problem  for  determining  the  neutron  form  factors  is  the  fact  that 
there  are  no  free  neutrons  to  scatter  from,  so  instead  of  using  a  simple  hydrogen 
target,  the  neutron  form  factors  are  determined  mostly  from  scattering  off 
a  deuteron  target.  Thus  the  scattering  is  actually  quasielastic  instead  of 
elastic.  The  deuteron  is  used  for  several  reasons:  it’s  the  simplest  nuclear 
target  available,  it  is  loosely  bound,  and  its  wave-functions  are  the  best  known. 
Since  the  deuteron  is  loosely  bound,  the  neutron  is  barely  off  the  mass-sheU 
(provided  the  neutron’s  recoil  momentum  pr  is  near  zero).  This  is  important 
since  differences  between  bound  and  unbound  nucleons  are  thought  to  be 
possible.  Of  course,  the  analysis  becomes  trickier  with  the  introduction  of 
deuteron  wave  functions  and  meson-exchange  currents  in  the  deuteron,  not 
to  mention  differentiating  between  proton  and  neutron  scattering.  Since  this 
experiment  dealt  with  proton  form  factors,  we  were  spared  these  additional 
problems. 

1.3  Form  Factor  Models 


Before  looking  at  the  existing  data,  a  discussion  of  the  many  form  factor 
models  is  in  order.  While  the  models  are  usually  based  somewhat  on  physical 
arguments,  they  are  all  at  least  partially  fitted  to  the  measured  data.  The  most 
common  approach  is  to  try  to  fit  a  single  underlying  form  factor  to  all  four  of 
the  form  factors,  though  in  actuality  it  is  the  proton  data  that  most  constrains 
the  models  since  its  form  factors  have  been  measured  more  precisely. 
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1.3.1  Dipole  Fit 


The  most  commonly  used  fit  is  the  dipole  fit.  Its  form  is; 


Gd{Q^) 


1 

(1 + 


(1.9) 


A  is  a  fit  parameter  whose  value  is  0.71  (GeV/c)^  based  on  the  existent  data. 
The  dipole  fit  is  so  popular  that  data  and  other  theories  are  often  normalized 
with  respect  to  it.  Though  it  works  fairly  well,  it  is  completely  empirical.  Note 
that  Gd(0)  =  1.  In  relation  to  Gd,  the  form  factors  are: 

GW)  =  GW) 

(1-10) 

GUQ")  =  th^GoW^)  I 


hence  these  three  form  factors  seem  to  scale  in  relationship  to  each  other.  Notice 
that  they  go  to  the  correct  values  as  goes  to  zero.  The  dipole  form  is  not 
appropriate  for  G'^  as  at  =  0  its  value  must  go  to  zero.  For  the  best 
empirical  fit  seems  to  be: 

guq^)  =  ^^Gd{Q^)  .  (1.11) 

For  a  Paris  potential  deuteron  model,  the  fit  values  are  a  =  1.25  ±  0.13  and 
h  =  18.3  ±  3.4. 


1.3.2  Vector  Dominance  Model  Fits 


A  model  that  has  been  fruitful,  though  ultimately  not  definitive,  is  the 
vector  dominance  model  (VDM).  In  this  model,  the  form  factors  provide  a 
physical  understanding  of  nucleon  structure  based  on  dispersion  relations.  As 


t®]  B.  Frois  and  S.  Platchkov,  Modem  Topics  in  Electron  Scattering,  eds.  B.  Prois  and 
I.  Sick,  World  Scientific,  (1991). 
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depicted  in  figure  1.2,  the  virtual  photon  converts  to  a  neutral  vector  meson 
(which  may  be  isoscalar  or  isovector)  prior  to  coupling  with  the  nucleon.  The 
absorptive  parts  of  the  form  factors  (the  spectral  function)  are  parameterized 
by  pole  terms  at  vector  resonances,  so  that  the  form  factors  can  be  expressed 
as  a  sum  over  meson  propagators  times  coupling  strengths: 


IS(V) 


pism  ^  ^ _ ifi _ 


(1.12) 


(It  is  the  coupling  strengths  and  masses  Mjs{v)-  that  are  usually  fitted.) 

The  above  form  factors  are  related  to  the  standard  Dirac  and  Pauli  ones  by: 

Ff  =  +  F/^ 


=  -{Fi^  +  Fi^) 

Kp 

F”  =  F/-®  -  Fl'' 


(1.13) 


where  Kp  =  1.793  and  /c„  =  —1.913.  (Some  fits  use  different  normalizations.) 


The  VDM’s  need  for  resonances  led  to  the  predictions  of  the  existence 
of  the  isovector  vector  meson  p  and  the  isoscalar  vector  meson  a;,  The 

isoscalar  strange  (f)  meson  is  usually  included,  as  well.  A  good  fit  cannot  be 
achieved  using  the  vector  mesons  alone,  however,  and  all  of  the  fits  include 
additional  empirical  poles.  Classic  examples  of  this  type  of  fit  are  the  Mainz 
eight-pole  fit  and  the  fit  of  Hohler  et  aL,  which  treated  the  p  exchange 


L.  Andivahis  et  aLj  Phys.  Rev.  D  60,  5491  (1994). 

T.A.  GrifFy  and  L.I.  SchifF,  High  Energy  Physics,  Vol.  1,  ed.  E.H.  Burhop,  Academic 
Press,  (1967). 

G.G.  Simon  et  aL,  Nucl.  Phys.  A333,  381  (1980). 

[9]  G.  Hohler  ei  a/.,  Nucl.  Phys.  B114,  505  (1976). 
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Figure  1.2  Diagram  of  Vector  Dominance  Model  (VDM)  from  reference  . 

more  completely  than  a  simple  pole  term,  taking  the  effects  of  the  two-pion 
continuum  into  account. 

Some  VDM  models,  called  extended  vector  dominance  models  (EVDM), 
describe  the  effective  interaction  of  the  photon  with  the  nucleon  as  having 
two  pieces:  the  term  involving  vector  mesons  and  a  term  involving  direct 
interactions  with  the  nucleon.  In  the  time-like  region  (q^  >  0)  the  vector  meson 
contributions  dominate  near  the  meson  poles.  However,  in  the  space-like  region 
<  0)  where  the  intermediate  mesons  are  far  off-shell  the  model  must  be 
extended  to  include  the  non-negligible  direct  contributions. 

A  particularly  successful  EVDM  model  at  intermediate  energies  is  by 
Gari  and  Krumpelmann.  The  Gari-Krumpelmann  model  incorporates 

quark  dynamics  at  large  via  perturbative  QCD  (pQCD).  Notice  the  model’s 


[10]  p  Mergell  et  al,  hep-ph/95 06375  preprint  (20  Jun  95). 

M.  Gari  and  W.  Krumpelmann,  Z.  Phys.  A322,  689  (1985). 
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isovector  Fi  term; 


where  mo  =  0.776  GeV, 


Fi{Q^)  = 


A?  +  g2  Al  +  Q- 


Q^  =  Q‘ 


log(^) 

^QCD 

'o8(aT^) 


(1.15) 


(1.16) 


This  gives  a  meson  monopole  fit  at  low  which  dies  out  as  increases  to 
leave  a  direct  photon-quark  coupling  term.  Fi  gives  the  expected  behavior  at 
high  Q^,  such  as  Q~^  scaling  (which  will  be  discussed  in  the  next  section).  Both 
meson  and  quark  features  exist  for  ~  A^.  Fits  to  the  data  gave  the  values: 
Qp/fp  =  0.377,  Ai  =  0.795  GeV,  Aj  =  2.27  GeV,  and  Aqcd  =  0.29  GeV. 

1.3.3  Quark  Models 


Many  efforts  have  been  made  to  model  the  form  factors  at  high  based 
on  quark  theory.  Dimensional  scaling  laws  predicted  that 


G{Q^  oo) 


=  (^) 


(1.17) 


where  uh  is  the  number  of  valence  quarks,  Perturbative  QCD  was  later 
shown  to  give  this  scaling  result: 


GM(g=)  =  c=f!)[.n(^); 


-4/3^ 


(1.18) 


02]  s_j.  Brodsky  and  G.R.  Farrar,  Phys.  Rev.  D  11,  1309  (1975);  S.J.  Brodsky  and  B.T. 
Chertok,  Phys.  Rev.  D  14,  3003  (1976). 

03]  G.P.  Lepage  and  S.J.  Brodsky,  Phys.  Rev.  Lett.  43,  545,  1625(E)  (1979);  G.P. 
Lepage  and  S.J.  Brodsky,  Phys.  Rev.  D  22,  2157  (1980). 
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and  implied  a  decay  of  GmQ^  at  high  as  given  by  the  running  coupling 
constant  as{Q^): 

where  ^0  is  a  QCD  one-loop  function  that  equals  9  for  a  three-flavor  quark 
system.  The  rate  of  decrease  is  determined  by  the  scale  parameter  AgcD* 
C  is  a  model-dependent  constant  of  unreliable  value.  Nonetheless,  a  recent 
update  to  Hohler  et  al  incorporated  this  large-Q^  pQCD  behavior  into 
their  VDM  model.  Their  flt  claimed  to  have  poles  identifiable  with  physical 
states  in  all  but  one  case.  Those  poles  were:  a;,^,5'(1600),p  (with  two-pion 
cut),  p(1450),p(1650),  and  p(1690).  The  p(1690)  is  their  one  claimed  unphysical 
state.  According  to  the  Particle  Data  Book,  the  mass  of  the  second  p  is 
1700  MeV,  not  1650. 

Unlike  perturbative  QCD,  QCD  sum  rules  axe  able  to  determine  the 
overall  magnitude  of  Cjf  and  not  just  its  evolution.  This  method  uses  sum  rules 
to  estimate  wave  function  moments  and  then  the  wave  functions  themselves; 
after  which  perturbative  calculations  can  be  performed.  These  wave  functions 
are  quite  different  from  nonrelativistic  or  asymptotic  wave  functions  and  imply 
that  in  the  proton  the  u-quark  whose  spin  is  parallel  to  the  proton’s  momentum 
carries  around  65%  of  the  proton’s  total  momentum, 

Additionally,  the  center-of-mass  bag  model  has  also  been  used  to  predict 
the  form  factors.  In  it,  the  SU(6)  symmetry-breaking  spin-dependent  forces 
arise  either  from  pQCD  or  from  nonperturbative  instanton  effects.  The  fitted 
parameter  is  the  ratio  between  the  d-quark  and  u-quark  distribution  radii  in 

M.  Aguilar-Benitez  et  al,  Phys.  Rev  D  60,  1173  (1994). 

[15]  V.L.  Chernyak  and  A.R.  Zhitnitsky,  Phys.  Rep.  112,  173  (1984). 

X.  Song  and  J.S.  McCarthy,  Phys.  Rev.  C  46,  1077  (1992). 
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the  proton.  Small  pion-cloud  conections  are  needed  for  <  0.3GeV^/c^, 
1.4  Method  of  Rosenbluth  Separation 


An  of  the  form  factor  measurements  presented  in  the  next  section 
with  the  exception  of  one  (SLAC-NE4)  were  performed  using  the  Rosenbluth 
separation  method.  Defining  the  quantity 

e(de)  =  [l  +  2(l  +  r)tan2(<?/2)]-'  ,  (1.20) 


equation  1.7  can  be  rewritten  as 


dor 

d<r 

1 

Rosenbluth 

dCle_ 

n.l+'T  . 

(1.21) 


e  is  the  virtual  photon’s  longitudinal  polarization,  ranging  from  0  to  1  as  goes 
from  TT  ^  0®.  For  Rosenbluth  separation  it  is  convenient  to  define  a  reduced 
cross  section  that  is  a  function  of  and  c: 


dor 

_6(1  + 

't) 

da 

dfle 

reduced 

da 

dflg 

dQc 

ns 

^  J  Rosenhluih 


(1.22) 


=  rGl,{Q^)  +  eGim 


Using  this  equation  a  set  of  measurements  at  fixed  but  different  values  of 
E  and  6^  can  be  used  to  determine  the  form  factor  values.  The  reduced  cross 
sections  are  plotted  versus  e  with  tG\i  as  the  intercept  and  G\  as  the  slope. 
Figure  1.3  shows  the  Rosenbluth  plots  of  Walker  et  ah 


The  problem  of  determining  the  form  factors  from  kinematically  different 
measurements  shoidd  be  apparent.  The  extracted  form  factors  are  quite 
sensitive  to  systematic  errors  in  E^E\  and  0.  For  example,  one  can  maximize 
the  count  rate  a:nd  longitudinal  polarization  at  forward  angles,  but  there  the 


[1*^1  R.C.  Walker  et  al,  Phys.  Rev.  D  49,  5671  (1994). 
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Figure  1.3  The  Rosenblutb  plots  Irom  reference  at  (a)  =  1*0,  (b) 

=  2.0,  (c)  =  2.5,  and  (d)  =  3.0  (GeV/cf.  The  error  bars  include  statistical 

and  systematic  errors  but  not  normalization  errors. 

cross  section  is  very  sensitive  to  variations  in  0.  Additionally,  spectrometer 
acceptances  large  enough  to  get  good  statistics  might  require  model  folding 
over  the  acceptance  range.  As  mentioned  earlier,  at  high  Q^, 
making  the  Rosenbluth  fit’s  slope  hard  to  determine. 

1.5  Previous  and  Data 

Figure  1.4  shows  selected  data  from  previous  G^  and  G^  measurements 
with  particular  emphasis  on  the  intermediate- values  0.5  — 3.0  (GeV/cf). 
The  solid  line  is  the  dipole  fit  (eqs.  1.9  and  1.10).  The  data  are  from:  Walker  et 
al.  (open  diamonds),  Bartel  et  al.  (open  squares),  SLAC-NE4  (small 


[18]  Bartel  et  al.,  Nucl.  Phys.  B58,  429  (1973). 

[19]  p.E.  Bosted  et  al.,  Phys.  Rev.  C  42,  38  (1990). 
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open  circles),  SLAC-NEll  (shaded  diamonds),  Janssens  et  al.  (Xs), 

and  Holder  et  al.  (open  circles).  The  SLAC-NE4  data  are  only  for  That 
experiment  was  done  at  6^  =  180°,  where  the  cross  section  depends  entirely  on 
(as  can  be  seen  in  eq.  1.7,  taking  the  nonstructure  term’s  angle-dependence 
into  account).  Figure  1.5  shows  the  same  data  set  normalized  to  the  dipole  fit. 

For  the  Mainz  (short  dash)  ^lnd  Mergell  et  al.  (thin  line)  VDM 
fits  and  the  extended  VDM  fit  of  Gari  and  Krumpehnann  (dots)  seem  to 
work  well.  The  other  fits  are  Janssens  et  aVs  three-pole  fit  (thick  dots). 
Holder  et  ol.’s  VDM  fit  8.2  (dot-dash),  and  Blatnik  and  Zovko’s  extended 
VDM  model  with  asymptotic  constraints  (long  dash).  For  G^  the  data  from 
the  difierent  measurements  seems  to  diverge  above  1  GeV^  and  a  conclusion  is 
difficult  to  reach.  The  remainder  of  this  section  will  discuss  the  proton  form 
factors  in  the  low  and  high  regions.  For  figures  and  discussion  of  G^  and 
G'^  measurements  please  see  Markowitz, 

1.6.1  Low  Data 


In  the  low-Q^  region  much  of  the  form  factor  measurement  effort  has 
been  focussed  on  obtaining  the  proton’s  root-mean-square  (rms)  charge  radius. 
In  the  Breit  frame  (section  A.2),  and  the  form  factors  axe  the  Fourier 


P.E.  Bosted  et  al,  Phys.  Rev.  Lett.  68,  3841  (1992). 


[21]  T.  Janssens  et  al,  Phys.  Rev.  142,  922  (1966). 


[22]  S.  Blatnik  and  N.  Zovko,  Acta  Phys.  Austr.  39,  1  (1974). 


[2^]  P.  Markowitz,  Ph.D.  dissertation,  William  &  Mary  (1992),  unpublished. 
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transforms  of  the  charge  and  magnetic  moment  distributions: 

Gsif)  =  /  pi£)e^^-^d^x 

J  (1.23) 

Gm(s")(#X!T=  , 

where  G£’(Af)(0)  is  given  in  eq.  1.10.  These  spatial  distributions  are  not  Lorentz 
invariant,  though,  and  consequently  are  only  useful  in  the  nonrelativistic  limit 
at  low-Q^ .  Since  the  form  factors  are  functions  of  a  single  variable  they  can  be 
expanded  in  a  Taylor  series  at  q  ^  =  0: 


(1.24) 


Expanding  the  Fourier  integral  for  the  charge  form  factor  gives: 

Gsif)  =  /  +  iq-x-  +  ...^  p{B)d^x 


1  +  *9  •  ® 


p{x)d^x 


(1.25) 


1  -  gl?!^  <  > +•••  > 


where  <  >=  /  r^  p{r)4Trr^  dr.  Matching  terms  in  the  two  previous  series 

equations  gives  the  relationship: 


(1.26) 


Applying  this  formula  to  the  dipole  fit  gives  the  value  of  the  proton’s  electric 
rms  radius  as  0.81  fm.  Early  Stanford  data  gave  a  value  near  this  (0.805  ± 
0.011),  W  but  a  later  result  from  Mainz  determined  by  fitting  the  low-Q^  data 
to  a  polynomial  of  second  order  gave  0.862  ±  0.012  fm.  This  is  consistent 
with  new  NIST  (National  Institute  of  Standards  and  Technology)  measurements 
using  Fourier-Bessel  analysis  which  gave  a  value  of  0.865  ±  0.020  frn.  Recent 


R.G.  Sachs,  Phys.  Rev.  120,  2256  (1962). 

[25]  p_  Borkowski  ei  al.,  Z.  Phys.  A276,  29  (1975). 

M.  McCord  et  al.,  Nucl.  Inst.  Meth.  B56,  496  (1991). 
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Lamb  shift  measurements  are  consistent  with  the  larger  proton  size.  This 
discrepancy,  unfortunately,  limits  the  attainable  precision  of  some  important 
experimental  tests  of  QED. 

1.5.2  High-Q^  Data 

In  the  high-Q^  region  it  is  much  easier  to  measure  Gm  than  Ge  due 
to  the  factor  of  t  in  the  Rosenbluth  cross  section.  Above  =  5  (GeV/c)^, 
Ge  only  contributes  a  few  percent  to  the  cross  section.  Arnold  et  al, 
measured  G^  up  to  =  31.3  (GeV/c)®  at  the  Stanford  Linear  Accelerator 
Center  (SLAG).  In  their  analysis  they  presumed  that  G^  continued  to  scale  with 
since  any  small  deviation  from  this  would  have  little  effect  at  such  high 
energies.  They  found  (figure  1.6)  that  at  around  =  8  (GeV/c)^  G^  scaled 
with  respect  to  Q~^  as  the  naive  dimensional  scaling  laws  said  and  that  after 
that  there  was  a  decrease  of  Q^G^j^  with  increasing  as  predicted  by  pQCD. 
MergeU  et  aVs  fit,  however,  shows  that  sizeable  hadronic  pole  contributions 
still  exist  between  10  and  20  (GeV/c)^  and  that  more  data  is  probably  needed 
before  concluding  that  scaling  has  been  seen. 

The  harder  term  has  been  measured  up  to  =  8.83  (GeV/c)^ 
by  Bosted  et  al.  Their  results  (figure  1.7)  show  that  G^  is  most 

consistent  with  the  dipole  fit  at  high-Q^.  They  found  G^,  on  the  other  hand, 
to  be  most  consistent  with  VDM  models,  particularly  the  one  by  Gari  and 
Krumpelmann. 

In  summary,  then,  there  exists  no  model  that  is  totally  successful  at 


[27]  D.J.  Berkeland  et  al.,  Phys.  Rev.  Lett.  76,  2470  (1995). 


t28]  R.G.  Arnold  et  al.,  Phys.  Rev.  Lett.  67,  174  (1986);  A.F.  Sill  et  al,  Phys.  Rev.  D 
48,  29  (1993). 
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Figure  1.6  Scaling  of  at  bigb-QK  Tie  fits  use  pQCD  (BL) 

and  sum  rules  (CZ)  witb  Aqcd  =  100  MeV, 

explaining  all  (or  even  one)  of  the  form  factors  over  the  entire  range. 

1.6  Motivation  for  this  Experiment 

Hopefully  by  this  point  our  motivation  for  this  experiment  is  understood. 
Ge  and  Gm  contain  the  unknown  knowledge  about  nucleon  structure  and  hence 
need  to  be  known  precisely  to  constrain  nucleon  models.  At  intermediate- 
energies,  knowledge  of  the  form  factors  helps  to  model  the  crossover  area 
between  a  low-energy  description  in  terms  of  meson  degrees  of  freedom  and 
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0  2,  4  6  8 

q2  [{GeV/c)2] 


Figure  1.7  Proton  form  factors  at  higb-Q^.  Please  see  reference  for 
explanation  of  Sts. 


a  high-energy  description  in  terms  of  quarks  and  gluons.  The  systematic 
uncertainties  involved  with  the  Rosenbluth  method  and  the  limited  precision  of 
the  data  at  high-Q^  point  to  the  need  for  a  new  method  of  form  factor 
measurement.  Our  method,  polarization  transfers  in  the  (e,e'p)  reaction, 
allowed  for  the  simultaneous  measurement  of  the  form  factors,  reducing  the 
systematic  errors.  Performing  such  an  experiment  at  Bates  Linear  Accelerator 
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Center  also  served  to  lay  the  ground  work  for  similar  future  experiments  at 
the  next  generation  of  electron-scattering  labs  such  as  the  Continuous  Electron 
Beam  Accelerator  Facility  (CEBAF).* 


[29]  CEBAF  PR  93-027, “Electric  Form  Factor  of  the  Proton  by  Recoil  Polarization”,  C.F. 
Perdrisat  and  V.  Punjabi,  co-spokespersons. 


*  CEBAF  was  renamed  the  Thomas  Jefferson  National  Accelerator  Facility  on  May  24, 
1996. 
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The  (c,e'p)  Reaction 


2.1  Electron-Proton  Coincidence  Scattering 


Measurement  of  proton  polarizations  necessitates  use  of  the  coincidence 
(e,e'p)  reaction  rather  than  the  simpler  (e,e')  reaction  described  in  the  previous 
chapter.  Figure  2.1  shows  an  (e,e'p)  reaction.  (Our  reaction  was  actually  (e,e'p ) 
since  the  incoming  electron  and  outgoing  proton  polarizations  were  known,  but 
that  does  not  affect  the  coincidence  kinematics.)  The  angle  between  p'  and 
g'is  given  by 

9pq  =  cos“^  (p'  •  q)  .  (2.1) 


“Parallel  kinematics”  occur  when  p'  and  q  are  parallel  and  Opq  equals  0  or  180 
degrees.  The  angle  between  the  scattering  planes  of  the  electron  and  the  proton 
is 

^planes  ~ 


=  tan  ^ 


(AiX  A!')-(gxp')] 

L  |fc  X  X  p' 
p'-iqx  k±) 


p'  •k± 


(2.2) 


where  k±  ^  k  X  q.  The  first  equation  is  easier  to  visualize  while  the  second 
is  better  for  computational  purposes.  It  allows  the  quadrants  to  be  easily 
determined  since  the  numerator  (denominator)  is  positive  in  quadrants  1  and  2 
(1  and  4)  and  negative  in  3  and  4  (2  and  3).  If  <l>pianes  equals  0  or  180  degrees 
the  kinematics  are  called  “in-plane  kinematics”.  For  a  hydrogen  target  it  is 
always  the  case  that  9pq  —  0°  and  ^planes  —  0  . 


(35) 
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Measuring  the  outgoing  proton’s  polsirization  requires  a  second  scatter¬ 
ing.  In  our  case  this  is  done  by  scattering  the  protons  from  carbon  nuclei  in 
graphite  blocks  in  the  polarimeter.  The  angular  distribution  is 

i{e,<i>)  =  h{e)[\^Ac{e)pf^ -n]  ,  (2.3) 

where  Iq{0)  is  the  angular  distribution  for  no  polarization,  Ac{6)  is  the  analyzing 
power,  pb  is  the  transverse  polarization  of  the  proton  before  it  scatters,  and 


E.  Aprile-Giboni  ei  al,  Nucl.  Inst.  Meth.  216,  147  (1983). 
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Figure  2.2  The  angle^ averaged  (5° -20°  )  C  analyzing  power  as  a  function 

of  proton  energy. 

n(=  f  X  p"/\f  X  lies  in  the  plane  of  the  second  scattering.  For  our 
polarimeter’s  coordinate  system  (described  in  section  4.1.4)  this  leads  to: 

=  Ioi0)[l -p^^AciO)  sin  (p  +  p(PAc(6)  cos  <f>]  .  (2.4) 

The  analyzing  power  is  a  function  of  9  and  Tp,  the  proton’s  kinetic  energy. 
Figure  2.2  shows  the  angle-averaged  analyzing  power  as  a  function  of  Tp. 
Notice  its  peak  at  200  MeV.  Our  experiment’s  kinematics  were  chosen  so  that 
the  proton  kinetic  energies  would  range  near  this  value.  Figure  2.3  shows 
the  analyzing  power’s  ^-dependence  for  a  given  proton  energy.  At  small  angles 
(~<  3°)  the  ^^C(p,p')  elastic  cross  section,  as  shown  in  figure  2.4,  is  dominated 
by  multiple  Coulomb  scattering  which  has  zero  analyzing  power.  These  events. 


M.W.  McNaughton  et  al.,  Nucl.  Inst.  Meth,  A241,  435  (1985). 
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Figure  2.3  0-dependence  of  Ac{0)  for  a  given  proton  energy. 


however,  account  for  a  large  majority  of  all  the  scattering  events.  Hence  there 
is  a  need  for  the  experimental  hardware  to  reject  small-angle  events.  Since  the 
cross  section  decreases  with  the  useful  events  are  primarily  in  the  range  from 
5°-20^ 


Assuming  conservation  of  angular  momentum  and  parity,  an  electron 
that  is  longitudinally  polarized  will,  in  the  first  Born  approximation,  only 
transfer  longitudinal  and  transverse  polarizations  for  in-plane  kinematics. 
Normal  polarization  may,  however,  be  induced  by  final  state  interactions. 
Thus  normal  polarization  is  independent  of  the  electron  beam’s  helicity.  For 
hydrogen,  there  aren’t  any  final  state  interactions.  (Two-photon  exchange, 
however,  can  result  in  a  very  small  normal  polarization.) 
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Figure  2.5  Transformation  of  target  polarization  to  focal  plane  coordinates. 
X  is  the  angle  the  proton^s  spin  precesses  relative  to  its  momentum. 

Equation  2.4  allows  for  the  determination  of  only  two  of  the  three 
polarization  components.  If  precession  of  spin  occurs  due  to  magnetic  fields 
between  the  target  and  the  focal  plane,  however,  it  is  possible  to  extract  the 
third  polarization  component.  For  a  vertically-bending  proton  spectrometer, 
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one  must  taJce  into  account  the  rotation  displayed  in  figure  2.5: 

A 

“  C'OS'^Tlfgf 

Vfv  ”  ^^9^  (2.5) 

A 

i/p  =  -  sin  X  ritgt  +  cos  x  itgt 

where  x  is  the  angle  the  spin  precesses  relative  to  the  momentum,  (x  is  discussed 

A 

further  in  chapter  3.)  Since  ritgt  is  helicity-independent  while  Itgt  is  not,  they 
can  be  separated  by  making  measurements  at  both  helicity  states  using  the 
relationships: 


pT  = 


tgt 

Pi  = 


p^nn+h)+p(j>{-h) 

2cosx 

2sinx 


Equation  2.4  becomes: 


I{e,  (f>)  =  Io{e)  [l  -  (pjf*  cos X sin  x)^c(^)  sin  ^-p\^*Ac{6)  cos  ,  (2.7) 

where  the  term  should  be  near  zero  for  hydrogen.  Using  this  equation,  the 
polarizations  at  the  target  can  be  measured  if  one  knows  x  and  the  analyzing 
power.  As  is  shown  in  appendix  A,  knowing  pi  and  pt  allows  one  to  determine 
the  ratio  Ge/Gm  hy  the  relationship: 

pT  _  2M  Ge 

pfflf  -  (£;  +  £;f)tanf  Gm 

It  is  also  possible  to  measure  the  beam  helicity  with  the  FPP  using  the  target 
polarizations: 
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Experimental  Equipment  Setup 


3.1  Accelerator  Facility 

This  experiment  was  performed  at  the  Bates  Linear  Accelerator  Center 
-  a  facility  operated  by  the  Massachusetts  Institute  of  Technology  (M.I.T.) 
that  is  located  in  Middleton,  Massachusetts.  Figure  3.1  shows  a  layout  of  the 
laboratory.  The  accelerator  is  designed  to  accelerate  electrons  (either  polarized 
or  unpolarized)  in  the  100-1000  MeV  range  with  beam  currents  up  to  50  fiA. 
Beam  energies  above  approximately  490  MeV  are  achieved  by  routing  the 
beam  through  a  recirculator  for  a  second  pass  through  the  accelerator.  This 
experiment  used  polarized  electrons  of  approximately  580  MeV  and  weis  singles- 
rate  limited  to  currents  of  less  than  20  fiA.  The  beam  pulses  ranged  from  12  to 
17  fisec  long  and  the  beam’s  pulse  repetition  rate  was  600  Hz  for  a  duty  factor 
of  1%.  Bates’s  South  Hall  Ring,  a  pulse  stretching  device,  should  produce  a 
duty  factor  of  85%  but  was  not  available  for  our  use. 

This  experiment  used  the  South  Hall’s  beam  line  B,  which  is  detailed  in 
figure  3.2.  The  beam  energy  is  determined  to  about  two-tenths  of  a  percent 
using  the  calibrated  EBl  chicane  magnet.  Three  beam  line  toroids  -  BTl, 
BT2,  and  BT3  -  measure  the  charge  passing  through  the  beam  line.  BT3  was 
connected  to  a  Brookhaven  Instrument  Co.  (BIC)  current  integrator  which 
was  calibrated  by  injecting  a  known  current  into  a  single  loop  (“Q-loop”)  in 
the  toroid’s  transformer.  Charge  integration  was  started  and  stopped  by  the 
experimental  gate  and  the  integrator’s  output  was  recorded  in  a  scaler.  BTl 


(42) 


BATES  LINEAR  ACCELERATOR  CENTER 
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inear  Accelerator  Center. 
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and  BT2,  conversely,  went  to  a  charge  integrator  that  was  gated  by  the  AND 
of  the  experimental  gate  with  the  computer-not-busy.  This  integrator’s  output 
went  to  an  ADC,  allowing  the  charge  in  each  beam  pulse  to  be  individually 
measured.  Of  every  600  beam  bursts,  20  contained  no  electrons,  providing  a 
pedestal  measurement  for  the  ADC  which  was  subtracted  from  the  measured 
charge  to  determine  the  actual  charge  delivered.  BTl  and  BT2  were  calibrated 
in  the  same  manner  as  BT3  so  that  we  knew  the  nC/channel  in  the  ADC. 

Also  in  the  beam  line  are  several  beam  position  monitors.  (We  used 
two.)  These  monitors  are  actually  resonant  cavities  with  dipoles  at  both 
ends.  Depending  on  the  beam’s  position,  the  different  ends  sense  the  passing 
beam’s  electromagnetic  field  at  different  times,  allowing  the  beam  position  to 
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be  determined.  There  is  also  a  BeO  flip  target*  1  m  upstream  from  the  target 
for  viewing  the  beam’s  position.  Photomultiplier  tubes  were  used  to  monitor 
the  beam’s  halo.  One  was  located  near  the  target  while  another  was  near  the 
M0ller  polarimeter.  The  polaximeter  is  the  last  major  component  of  the  B  beam 
hne  and  will  be  described  in  the  next  section. 

3.2  Polarized  Beam 

3.2.1  Polarized  Source 

This  experiment  used  a  polarized  beam  of  25  to  35  percent  polarization  and 
currents  up  to  about  20/iA.  The  accelerator’s  polarized  injector  is  illustrated  in 
figure  3.3.  The  longitudinally  polarized  electrons  are  produced  from  a  Gallium- 
Arsenide  crystal  using  circularly  polarized  laser  light.  The  electrons  are  then 
accelerated  to  60  kV  in  the  photocathode  gun  and  then  an  additional  320  kV 
in  the  source’s  accelerating  column.  A  Wien  filter  can  be  used  at  this  point  to 
rotate  the  beam’s  polarization.  Two  45®  magnets  then  bend  the  beam  to  inject 
it  into  the  main  accelerator  line.  At  this  point  its  380  kV  injection  energy  is 
the  same  as  that  of  the  unpolaxized  beam  when  it  comes  out  of  its  thermionic 
source.  The  injector  is  differentially  pumped  using  ion  pumps  (which  have  no 
moving  parts)  with  a  vacuum  in  the  vicinity  of  the  crystal  of  «  10“^^  torr  for 
hydrogen  and  ft!  10“^®  torr  for  more  damaging  water  and  carbon. 

The  laser  light  used  in  the  injector  is  from  a  9  watt  Ti-Sapphire  laser 
which  is  optically  driven  by  a  30  watt  CW  argon  laser.  The  argon  laser  light  is 
pulsed  with  a  repetition  rate  of  600  Hz  using  a  mechanical  chopper  (a  slotted 


*  A  flip  target  is  a  target  that  can  quickly  be  put  into  the  beam  line  with  the  press  of 
a  button  or  the  flip  of  a  switch,  but  has  only  one  setting  rather  than  an  adjustable 
position  range. 

M.  Farkhondeh,  private  explanatory  tour  of  source. 
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Figure  3.3  Layout  of  polarized  source. 

wheel),  which  increases  the  power  output  of  the  Ti- Sapphire  laser.  The  argon 
laser’s  output  lies  in  the  450-500  nm  range  while  the  Ti-Sapphire  laser’s  light  is 
at  or  near  752  nm.  The  laser  beam’s  intensity  can  be  modified  from  the  control 
room  prior  to  entering  the  injector.  Pockels  cells  are  used  to  set  the  beam 
burst  size  and  to  change  the  linearly  polarized  light  to  right  or  left  circularly 
polarized.  (A  Pockels  cell  is  an  electro-optical  device  containing  a  birefringent 
crystal  whose  phase  retardation  is  proportional  to  the  applied  voltage,  allowing 
the  polarization  of  the  incident  beam  to  be  changed,  )  The  beam’s  helicity 
is  flipped  on  a  pulse-to-pulse  basis.  Each  cycle  of  the  60  Hz  line  voltage  is 
divided  into  10  time  bins.  The  hehcity  for  the  first  10  bins  is  (pseudo-)randomly 


E.  Hecht,  Optics,  Addison- Wesley,  (1987). 
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determined  while  the  next  10  bins  are  the  complement  of  the  previous  10.  Then 
there  are  random  bins  again  and  so  on. 

The  GaAs  crystal  has  a  band  gap  of  1.52  eV,  vrith  the  S1/2  state  serving 
as  the  conduction  band  and  the  P3/2  state  serving  as  the  valence  band.  For 
positive(negative)  helicity  light,  optical  selection  rules  require  Am^-  =  +1(— !)• 
For  positive(negative)  helicity  light  three  times  as  many  electrons  go  to  the 
mj  =  — l/2(+l/2)  state  as  to  the  mj  =  +l/2(— 1/2)  state,  giving  a  theoretical 
polarization  of  50%.  We  typically  had  polEurization  of  25%  to  35%. 

Strained  crystals  (thin  GaAs  grown  on  a  GaAs(i_a,)Pie(=o-»o.29)  buffer)  remove 
the  m-state  degeneracy  in  the  P3/2  state  so  that  polarizations  on  the  order 
of  80%  can  be  achieved.  This  was  tried  for  our  experiment  but  the  quantum 
efficiency  (the  number  of  electrons  produced  per  incident  photon)  was  too  low. 

The  GaAs  crystal  is  coated  with  a  thin  film  (<  1  atomic  layer)  of  Cs 
and  NF3  to  reduce  its  work  function  and  allow  electrons  in  the  conduction 
band  to  be  emitted  from  the  crystal.  The  quantum  efficiency  of  the  crystal  was 
around  1%.  The  quantum  efficiency  is  very  sensitive  to  surface  contamination  so 
occasionally  the  crystal  needs  to  be  reactivated.  This  procedure  takes  about  an 
hour.  First  the  surface  of  the  crystal  is  heated.  Then  it  is  sprayed  alternatively 
with  Cs  and  NF3  several  times. 

3.2.2  Spin  Precession  in  Beam 

After  coming  out  of  the  polarized  source,  the  longitudinally  polarized 
beam  goes  through  a  360°  bend  in  the  recirculator  (if  the  energy  desired  is  high 
enough  to  warrant  this)  and  a  90°  bend  to  get  into  the  South  Hall.  Due  to 
the  precession  of  spin  this  vnll  change  the  direction  of  the  beam’s  polarization. 


D.T.  Pierce  et  al,  Rev.  Sci.  Inst.  61,  478  (1980). 
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This  can  be  compensated  for  by  using  the  Wien  filter,  which  rotates  the  spin 
using  crossed  electron  and  magnetic  fields.  There  are,  however,  certain  “magic” 
energies  which  do  not  require  the  use  of  the  Wien  filter.  Our  energy,  574  MeV 
was  one  of  these. 


The  spin  precession  of  a  charged  particle  traveling  in  an  electromagnetic 
field  is  given  by  the  Thomas  equation: 


ds  _  e  _ 
dt  me 


In  our  case  E  =  0  and  /3  •  B  is  small  so  the  spin  precession  reduces  to 


—sx  +  =  3  XOJ, 

me  \2  'iJ 

—  +  . 

me  \2  7/ 


Likewise,  the  equation  of  motion  for  a  charged  particle  moving  in  an  electro¬ 
magnetic  field  is 

^  =  e{E+-xB)  .  (3.3) 

at  c 

Using  ^  =  0  again  and  the  relationship  p  =  mjv  gives 

^=t{-£-xB)  =  pxiS, 

dt  nfmc  /g 

Up  = - B  . 

'ymc 

The  total  spin  precession  with  respect  to  the  beam’s  direction  can  now  be 


written  as: 


A  f  If  ^  % 

Au  =  Ug  -  Up  =  — (-  -  1) 


mc^2 


- 1) 


J.D.  Jackson,  Classical  Electrodynamics,  John  Wiley  &  Sons,  (1975). 
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There  are  two  areas  in  the  accelerator  where  the  beam  is  bent:  the 
recircnlator  (if  used)  and  the  90®  entrance  to  the  South  HaU;  so  that  the  final 
spin  angle  relative  to  the  beam  is 

Axr"“'  =  (|-l)(xr'”T”“"  +  x;“"7'‘"“)  .  (3-6) 

where  The  “preaccelerator”  beam  energy 

is  20  MeV.  Therefore  the  equation  relating  the  final  beam  energy  with  the 
beam  energy  after  the  first  pass  through  the  accelerator  but  before  being  bent 
through  the  recirculator  is  -  20)  so  that  7"'=’'’‘=  can 

be  expressed  as:  =  ’'/2> 

“magic”  energies  can  now  be  determined  by  setting  =  mr.  The  South 

Hall’s  “magic”  energies  are  280  (which  doesn’t  use  the  recirculator),  574,  and 
868  MeV  for  n  =  1,2,3.  We  used  574  MeV  for  our  beam  energy. 

3.2.3  M0ller  Polarimeter 

Approximately  20m  upstream  from  the  target  on  beam  line  B  is  a  M0ller 
polaximeter  to  measure  the  beam  polarization.  It  employs  elastic  scattering  of 
the  spin  polarized  electron  beam  with  the  spin  polarized  M0ller  target  electrons. 
The  device  ^lnd  the  theory  behind  its  operation  are  described  in  Arrington  et 
ah  and  Wagner  et  ah 

The  cross  section  for  polarized,  elastic  electron-electron  scattering  may 
be  written  as 

=  +  .  (3.7) 

J.  Arrington  et  aL,  Nuc.  Inst.  Meth.  A311,  39  (1992). 

B.  Wagner  et  aL^  Nuc.  Inst.  Meth.  A294,  541  (1990). 
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where  is  the  unpolarized  cross  section,  Pb{Pt)  polarization 

components  of  the  beam(target),  and  Aij  are  the  nine  asymmetries,  z  is  the 
direction  of  the  beson  and  x  lies  on  the  detected  scattering  plane.  In  the  extreme 
relativistic  limit,  the  asymmetries  are,  to  lowest  order  in  QED:  . 


Azz  = 


Axx  —  -^yy  — 


(7  +  COS^^cm)siP^^em 
(3  +  COs2  Oemy 
sin^  0. 


'em 


Axz  - - 


•^xy  — ^yx  —  -^yz  —  —  0 


(3  +  COs2  0cmY 

2  Sm  0cm  cos  0cm 
7(3  +  COs2  0cmy 
zy 


(3.8) 


where  7  is  the  Lorentz  factor  relating  the  lab  and  center-of-mass  frames.  By 
setting  0cm  =  90°  these  reduce  to  Azz  =  —  |  and  Axx  =  —Ayy  =  —  |  with  the 
cross-terms  going  to  zero. 


By  taking  measurements  in  both  longitudinal  directions  and  assuming 
that  the  beam  and  target  are  polarized  only  along  the  beam,  the  measured 
asymmetry  should  be 

A  =  AzzPhPT  ,  (3.9) 


where  P^  =  Pt  cos  Ot.  {6t  is  the  target  angle.)  Since  the  target  can  only  be 
magnetized  in  its  plane,  it  was  typically  tilted  at  a  relatively  small  target  angle 
of  ±30°  with  respect  to  the  beam  so  that  the  longitudinal  polarization  could 
be  measured.  Assuming  cos  6b  ^  1  and  averaging  the  asymmetries  gives 


Pb  = 


9  A+aoo  ±  A_3o» 
7P7’Cos(30°)  2 


(3.10) 


Due  to  Mott  scattering  and  other  backgrounds,  the  measured  asymmetry 
is  actually 


1 
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where  B  is  the  background  rate  and  S  is  the  M0ller  scattering  rate.  Our  signal 
to  background  ratio  was  approximately  2:1.  Using  the  two  Cerenkov  detectors 
in  coincidence  mode  would  have  reduced  background  but  would  have  imposed 
more  stringent  reqidrements  on  beam  position  stability,  detector  efficiency,  and 
mechanical  alignments.  Thus  we  used  them  for  separate  measurements  (singles 
mode).  Each  M0ller  run  included  both  a  so-called  real  scan  to  determine  the 
signal  to  background  ratio  and  a  peakscan  to  accumulate  statistics  for  the 
asymmetry  measurement.  The  scans  were  performed  by  varying  the  current 
of  the  quadrupole  used  to  bend  the  scattered  electrons  away  from  the  beam 
line  and  taking  measurements  for  set  periods  of  time.  The  polarization  was 
measured  about  once  a  day  or  whenever  the  polarized  source’s  crystal  was 
reactivated. 

Figure  3.4  shows  the  polarimeter’s  layout.  There  are  two  thin  foil  targets 
(13  fim  and  25  /xm  thick)  of  Supermendur,  an  alloy  consisting  of  49%  Fe,  49% 
Co,  and  2%  “Vacoflux”.  Also  in  the  target  ladder  are  BeO  and  A1  targets 
used  for  alignments  and  an  empty  target  to  allow  the  beam  to  pass  through. 
The  target  is  polarized  by  a  magnetic  field  from  Helmholtz  coils  around  the 
target.  The  field  reaches  nearly  150  Gauss  and  the  target  polarization  measures 
about  8%.  Downstream  from  the  target  is  a  lead  disk  collimator  to  define  the 
scattering  angle  of  interest.  Our  scattering  angle  was  2.42®,  which  is  90®  in  the 
center-of-mass  system.  Behind  the  collimator  is  a  25  cm  diameter  quadrupole 
magnet  which  bends  the  scattered  electrons  away  from  the  beam  line  to  the 
Cerenkov  detectors  while  letting  the  main  beam  pass  through  unaffected. 

3.3  Target 

This  experiment  was  the  first  to  use  a  new  target  built  at  M.I.T.-Bates 
called  the  Basel  Loop  target,  pictured  in  figure  3.5.  It  had  separate  liquid 
hydrogen  and  deuterium  targets  as  well  as  solid  BeO  (34.51  mg/cm^)  and 
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Quadrupole 

Target  Collimator  Magnet  Detectors 


Figure  3.4  Miller  polarimeter. 

C  (91.46  mg/cm*)  targets  and  a  blank  target.  The  liquid  target  loops  were 
identical  and  were  cooled  in  series  to  approximately  20  K  by  gaseous  He  from  a 
200  W  refrigerator.  The  temperatures  were  measured  with  diodes  and  carbon 
glass  resistors.  Internal  heaters  were  used  to  maintain  a  constant  temperature 
as  the  beam  was  turned  off  and  on.  Internal  fans  were  used  to  lessen  local 
boiling  and  to  improve  the  heat  exchanger’s  cooling  efficiency.  The  geis-handling 
systems  for  the  two  loops  were  separate  -  the  chief  difference  being  that  the 
deuterium  loop  was  designed  to  vent  to  a  storage  tank  due  to  deuterium’s 
expense.  Pressure  trcinsducers  were  located  at  the  gas  fill  lines  and  the  target 
pressure  was  ideally  slightly  above  1  atmosphere.  Each  loop  had  a  capacity 
of  1.5  liters.  The  thin  target  cell  W2ills  were  made  of  Havar,  which  has  the 
following  composition:  42%  Co,  19.5%  Cr,  19.1%  Fe,  12.7%  Ni,  2.7%  W,  2.2% 
Mo,  1.6%  Mn,  and  0.2%  C.  Further  target  properties  are  listed  in  table  3.1.  An 
IBM-compatible  personal  computer  in  the  counting  bay  was  used  to  control  the 
heaters  and  monitor  the  temperature  and  pressure. 

3.4  MEPS 


The  scattered  electrons  were  detected  in  the  Medium  Energy  Pion 
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Target 

LHz 

LDz 

■41I.IW1W 

5 

3 

Cell  Wall  Thickness  (fim) 

25.4 

10.16 

Cell  Wall  Density  (g/cm*) 

0.042 

0.017 

Nominal  Temperature  (K) 

20.3 

23.7 

Liquid  Density  (g/cm®) 

0.071 

0.16 

Liquid  Density  (g/cm®) 

0.36 

0.48 

Table  3.1  Target  Parameters 


Spectrometer  (MEPS).  The  layout  of  the  spectrometer  is  shown  in  figure 
3.6  and  some  of  its  properties  are  listed  in  table  3.2.  MEPS  is  a  QQSD 
spectrometer  (two  quadrupoles  followed  by  a  split  dipole).  The  focusing  in 
the  vertical  bend  plane  (the  momentum  dispersion  direction)  is  point-to-point 
((a;|fl)=0).  In  the  transverse  plane,  however,  it  is  parallel-to-point  {{y\y)  =0  and 
(y|^)^0),  allowing  the  scattering  angle  {(/>)  to  be  better  determined  since  the 
well  determined  focal  plane  y  position  is  directly  proportional  to  the  scattering 
angle.  (These  TRANSPORT  matrix  elements  are  described  in  more  detail  in 
chapter  4.)  The  magnets’  power  supplies  were  set  in  the  counting  room  via  a 
computer  dedicated  to  that  purpose.  The  dipole  field  was  measured  internally 
using  a  Rawson-Lush  probe  while  the  quadrupole  fields  were  assumed  to  scale 
linearly  with  the  magnets’  currents.  A  2-in  thick  lead  collimator  was  at  the 
front  of  MEPS.  The  dimensions  of  its  opening  were  ±2.59  cm  vertically  and 
±3.78  cm  horizontally  for  an  acceptance  of  14  msr. 
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Figure  3.6  MBPS  spectrometer  layout. 

The  MEPS  focal  plane  detector  system  consists  of  three  scintillators, 
an  aerogel  Cerenkov  detector,  and  a  two-plane  vertical  drift  chamber  (VDCX, 
where  the  ‘X’  stands  for  the  “crossed”  planes).  The  aerogel  had  an  index  of 
refraction  of  1.05  so  that  prirticles  traveling  at  or  above  95%  the  speed  of  light 
would  produce  Cerenkov  radiation.  Ten  photomultiplier  tubes  (PMTs)  aligned 
transverse  to  the  momentum  dispersion  direction  detected  light  in  the  diffusion 
cavity  above  the  aerogel.  This  detector  is  used  to  differentiate  electrons  from 
pions,  but  at  our  kinematics  there  were  few  pions  so  we  did  not  use  the  Cerenkov 
detector’s  information. 
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Maximum  Solid  Angle 

20  msr 

Momentum  Acceptance 

±10% 

Maximum  Momentum 

Momentum  Resolution 

5*10-^ 

Angular  Range 

35°  -  140° 

Radius  of  Curvature 

0.75  m 

Front  Drift  Distance 

0.53  m 

Flight  Path 

4.7  m 

Bend  Angle 

H- ^ 

o 

o 

Focal  Plane  Dispersion 

Table  S.2  MBPS  Parameters 


3.4.1  MEPS  Scintillators  and  Pilot 


The  three  scintillators  of  the  MEPS  focal  plane  are  made  of  NE-110 
plastic.  MSO  is  directly  above  the  VDCX  and  has  dimensions  of  17.8  cm  by 
58.4  cm.  It  is  1/4  inch  thick.  MS3  and  MS4  are  above  and  below  the  aerogel 
detector  respectively.  Both  have  an  active  area  of  20.3  cm  by  91.0  cm  and 
are  1/8  inch  thick.  AH  three  scintillators  are  each  connected  to  two  PMTs  via 
Lucite  light  guides.  The  PMTs  on  MS3  and  MS4  are  at  the  ends  but  for  MSO 
they  are  attached  next  to  each  other  on  the  side. 

A  coincidence  of  all  three  scintillators  was  used  to  generate  the  MEPS 
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Figure  3.7  MBPS  pilot  logic. 

pilot.  *  Figure  3.7  details  the  MEPS  pilot  logic.  After  the  six  scintillator 
PMT  signals  were  discriminated  the  signal-pairs  from  MS3  and  MS4  were 

*  Filoi  is  the  French  word  for  trigger.  We  referred  to  our  individual  spectrometer 
triggers  as  pilots,  reserving  the  word  trigger  for  the  final  data-readout  signal. 
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sent  to  a  meantimer,  which  makes  the  timing  relatively  independent  of  the 
the  event  position  along  the  device.  Due  to  the  side-by-side  placement  of 
the  PMTs  for  MSO,  their  signals  were  OR’ed  instead.  MS4MT  was  widened 
using  a  discriminator  and  MS3MT  was  delayed,  so  that  it  defined  the  timing 
of  the  3-fold  coincidence  MEPS  pilot  (MPI).  MPI  was  sent  to  OHIPS,  the 
proton  spectrometer,  via  a  900  105ns  cable  for  use  in  the  trigger  logic,  as 
well  as  to  be  scaled  and  have  its  timing  information  recorded  in  a  time-to- 
digital  converter  (TDC).  The  six  analog  PMT  signals  and  the  three  logic  signals 
(MSOOR,  MS3MT,  and  MS4MT)  were  sent  upstairs  to  the  counting  bay  on  500 
cable  where  they  were  scaled  and  sent  to  a  TDC.  The  analog  signals ’s  pulse 
shape  information  was  also  recorded  on  an  analog-to-digital  converter  (ADC). 
Amplification  was  necessary  to  send  the  signals  upstairs.  An  upstairs  (counting 
bay)  version  of  MPI  was  also  created  to  check  for  discrepancies  against  the 
downstairs  (South  Hall)  one. 

3.4.2  MEPS  VDCX 


The  MEPS  vertical  drift  chamber  has  two  planes  of  128  wires  each. 
The  wire  planes  of  the  two  chambers  are  3.81  cm  apart.  They  are  oriented 
orthogonally  to  one  another  and  45®  relative  to  the  dispersive  and  transverse 
directions.  The  active  area  of  each  plsme  is  7  X  25  inches.  The  anode  wires 
cire  20  fim  diameter  gold-plated  tungsten  wire  placed  4.23  mm  apart.  50  /xm 
diameter  beryllium-copper  guard  wires  separate  the  anode  wires  forming  field 
“cells”,  as  shown  in  figure  3.8.  Aluminized  mylar  at  high  voltage  («  —9.0  kV) 
lies  on  both  sides  of  each  plane.  The  entire  chcimber  is  tilted  with  respect  to 
the  spectrometer’s  central  momentum  ray  at  an  angle  of  approximately  45°  so 
that  a  particle  track  will  generate  signals  in  several  “cells”.  1®®J 


[3®]  D.  Caditz,  B.S.  Thesis,  M.I.T.  (1983),  unpublished  (Bates  Internal  Report  #85-04). 
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Figure  3.8  VDCX  layout  with  Reid  lines. 

The  chambers  operate  with  a  gas  mixture  of  roughly  50%  argon  and 
50%  isobutane.  When  a  charged  particle  travels  through  the  chamber  it  ionizes 
argon  atoms.  The  released  electrons  are  accelerated  by  the  electric  field  towards 
the  anode  wire  and  quickly  achieve  a  terminal  drift  velocity  of  fs5  cm//xs.  The 
electrons  also  ionize  other  argon  atoms  producing  an  avalanche.  The  isobutane 
is  used  to  absorb  photons  radiated  by  the  excited  argon  atoms.  Otherwise 
the  photons  can  cause  additional  avalanches.  Near  the  anode  wire,  however, 
the  rapidly  increasing  electric  field  overcomes  the  quenching  and  the  avalanche 
created  is  detected  as  a  voltage  pulse  in  the  wire,  The  maximum  drift  time 
for  the  MEPS  planes  is  around  250  ns.  The  particle  trajectories  can  then  be 
determined  by  using  the  drift  times  recorded  from  the  different  2inode  wires. 
(This  is  discussed  in  the  next  chapter.) 


[39]  W.R.  Leo,  Techniques  for  Nuclear  and  Particle  Physics  Experiments,  Springer- 
Verlag,  (1987). 
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The  MBPS  VDCX  is  instrumented  with  the  LeCroy  4290  Drift  Cham¬ 
ber  Operating  System  (DCOS),  consisting  of  16  channel  LeCroy  2735  amplifi¬ 
er/discriminator  (a./d)  cards,  32  channel  4291B  TDCs,  a  4298  system  controller, 
and  a  4299  databus  interface.  Figure  3.9  displays  this  system.  Chamber  art¬ 
work  (applied  circuit  traces)  carries  signals  from  the  VDCX  wires  to  connectors 
for  16  channel  twisted-pair,  shielded  ECL  cable.  (We  typically  referred  to  the 
shielded  ECL  cable  as  “black”  cable  due  to  its  color  and  the  unshielded  cable 
as  “ribbon”  cable.)  These  cables  go  to  the  a/d  cards,  whose  threshold  can  be 
set  from  upstairs.  (We  set  them  at  3  volts.)  The  signals  then  leave  the  VDCX 
hut  on  ribbon  cables  to  the  TDCs  in  a  CAMAC  (Computer  Automated  Mea¬ 
surement  And  Control)  crate  on  the  back  of  MBPS.  The  TDCs  were  operated 
in  common  stop  mode  and  drift  times  for  all  the  wires  that  fired  during  an 
event  were  read,  allowing  multiple  particle  trajectories.  The  system  controller 
receives  data  from  the  TDCs  through  the  CAMAC  crate’s  back  plane  and  gen¬ 
erates  a  stream  of  packed  module  and  wire  data  words.  The  packed  data  is 
then  sent  upstairs  through  50-channel  fiat  ECL  cable  to  a  databus  interface. 
This  interface  is  necessary  for  the  CAMAC  readout  system  (which  is  described 
in  section  3.8).  The  DCOS  system  can  be  calibrated  using  software  to  gener¬ 
ate  starts  and  stops  within  the  system  controller  and  to  adjust  the  slopes  zind 
pedestals  of  the  TDCs. 

3.5  OHIPS 


The  proton  arm  for  our  experiment  was  the  One-Hundred  Inch  Proton 
Spectrometer  (OHIPS)  to  which  we  added  a  newly-built  focal  plane  polarimeter 
(FPP).  OHIPS  is  a  QQD  spectrometer  whose  magnets  are  controlled  upstairs 
in  the  counting  bay  and  whose  dipole  field  is  measured  with  a  Rawson-Lush 
probe  just  like  MBPS.  Also,  like  MBPS,  the  focusing  is  point-to-point  in  the 
momentum  dispersion  direction  ((®|5)=0),  but,  unlike  MBPS,  the  focusing  in 
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On  the  back  of  MEPS 


Upstairs  in  SHCB 


4291B  TDCs 


4298  Controller 


Figure  3.0  DCOS  system  for  the  MEPS  VDCX. 


the  transverse  direction  is  point-to-point  ((y|^)=0)  as  well.  This  allows  for 
better  resolution  of  the  scattering-point  along  the  beam  in  an  extended  target 
and  also  allows  the  target  walls  to  be  excluded.  The  first  quadrupole  in  OHIPS 
had  a  vertical  focus  and  the  second  a  horizontal  focus.  This  focus  is  referred 
to  as  a  “VH”  or  “reverse”  focus  and  provides  a  large  vertical  angle  acceptance 
but  a  smaller  horizontal  angle  acceptance  than  a  “HV”  or  “normal”  focus. 
The  collimator  at  the  front  of  OHIPS  was  5.1  cm  thick  and  was  beveled  with 
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Figure  3.10  CHIPS  layout. 


frontface  dimensions  of  ±4.129  cm  horizontally  by  ±12.263  cm  vertically  and 
backface  dimensions  of  ±4.271  cm  by  ±12.687  cm.  The  drift  distance  at  the 
front  of  OHIPS  can  be  adjusted  to  allow  for  increased  angular  acceptance  at 
the  expense  of  angular  resolution,  or  vice  versa.  Figure  3.10  gives  a  sideview 
of  OHIPS  and  table  3.3  lists  some  of  OHIPS’s  properties.  Please  note  that  by 
equation  3.5,  OHIPS’s  90®  bend  precesses  the  spin  relative  to  momentum  by 
the  angle  x  =  7161.4®. 
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Solid  Angle 

7.5  msr 

Momentum  Acceptance 

±5% 

Maximum  Momentum 

Momentum  Resolution 

Design 

Best  Obtained 

1-10“® 

1.4-10-® 

Angular  Range 

19°  -  140° 

Radius  of  Curvature 

100  in 

Front  Drift  Distance 

1.6  m 

Flight  Path 

9.3  m 

Bend  Angle 

90° 

Focal  Plane  Dispersion 

Table  3.3  OHIPS  Parameters 


Many  modifications  were  made  to  OHIPS  in  the  fall  of  1993  so  that 
the  FPP  cotdd  be  placed  on  top  of  it.  The  shielding  hut  was  redesigned 
and  made  larger.  To  keep  the  weight  under  control  the  shielding  walls  were 
made  thinner  but  it  was  still  necessary  to  add  cross  struts  to  the  shielding  hut 
platform  supports.  An  air  conditioner  and  a  dehumidifier  were  placed  inside 
the  shielding  hut  to  make  the  climate  suitable  for  operating  the  FPP  whose  a/d 
cards  generated  a  lot  of  heat.  Unfortunately,  the  air  conditioner  often  dripped 
water  on  the  FPP’s  discriminator  card  power  supply  so  a  gutter  system  had 
to  be  rigged  up.  A  new  support  structure  for  the  OHIPS  VDCX  was  built  to 
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facilitate  the  support  of  the  FPP  and  also  to  allow  for  better  alignment  between 
the  spectrometer,  the  VDCX,  and  the  FPP.  A  new,  larger  electronics  platform 
was  put  on  the  back  of  OHIPS  to  handle  the  FPP’s  additional  electronics. 
The  old  wooden  ladders  on  OHIPS  were  replaced  with  safer  molded  aluminum 
ladders  and  OHIPS  got  a  fresh  paint  job. 

3.5.1  OHIPS  Focal  Plane  and  Pilot 

The  proton  arm  focal  plane  detector  system  consisted  of  five  scintillators, 
a  VDCX,  and  the  focal  plane  polarimeter.  Figure  3.11  shows  the  focal  plane 
including  detail  of  the  polarimeter.  Scintillator  OSO  was  directly  above  the 
VDCX  while  FSl  was  right  below  the  first  FPP  wire  chamber.  FS2  consisted  of 
three  scintillators:  FS2AB,  FS2CD,  and  FS2EF;  which  were  positioned  above 
the  large,  top  FPP  wire  chamber  and  overlapped  slightly.  Each  scintillator  had 
two  PMTs  attached  to  it  via  Lucite  light  pipes.  Table  3.4  gives  the  dimensions 
of  the  scintillators. 


Scintillator 

Width  (cm) 

Length  (cm) 

Thickness  (mm) 

OSO 

20.3 

65.0 

5.0 

FSl 

37.0 

71.0 

6.4 

FS2(AB,CD,EF) 

145.0 

32.0 

13.0 

Table  3.4  OHIPS  Scintillator  Dimensions 


As  detailed  in  figure  3.12,  the  OHIPS  pilot  was  generated  by  a 
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Figure  3.11  OHIPS  focal  plane  including  FPP  detail. 


coincidence  of  the  the  three  scintillator  planes.  This  can  be  thought  of  as  the 
first  (and  simplest)  level  of  a  three-level  FPP  trigger  (which  will  be  described 
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fully  in  section  3.7).  The  ten  scintillator  PMT  signals  were  fanned-out  to 
ADCs  and  discriminators.  The  two  discriminated  signals  for  each  scintillator 
were  meantimed,  producing  five  meantime  signals  that,  together  with  the  ten 
discriminated  signals,  were  sent  to  a  Phillips  761  NIM/ECL  converter  before 
being  sent  through  ribbon  cable  and/or  a  LeCroy  4418  Programmable  Delay 
Module  (PDM)  to  a  LeCroy  4434  Scaler,  a  LeCroy  4300B/4303  FERA/FERET 
TDC,  and  a  LeCroy  2365  Programmable  Logic  Unit  (PLU).  Since  the  full  trigger 
electronics  take  over  900  ns  to  generate  the  TDC  start  signal,  the  logic  signals 
were  delayed  using  the  PDMs  and  long  rolls  of  ribbon  cable. 

In  the  PLU  the  three  FS2  meantimed-signals  were  OR’ed  to  produce 
FS20R.  This  was  looped  back  through  the  PLU  so  that  it  could  be  AND’ed 
with  OSOMT  and  FSIMT  to  produce  the  CHIPS  pilot  (OPI).  The  timing 
was  determined  by  OSOMT.  Using  software,  it  was  easy  to  change  the  PLUs 
programming  so  that  the  pilot  logic  could  be  set  up  differently.  With  all  of  the 
looping  back  through  PLUs,  delaying  through  ribbon  cables,  and  fanning-out 
through  PDMs  and  ECL/NIM  converters;  our  three-level  electronics  became 
rather  complicated.  A  more  in-depth  description  is  given  in  McIntyre  and 
Warren, 

3.5.2  CHIPS  VDCX 


The  CHIPS  VDCX  is  very  much  like  MEPS’s.  It  also  has  two  crossed 
planes  separated  by  3.81  cm.  The  chamber  is  tilted  approximately  45°  with 
respect  to  the  central  momentum  ray.  Each  plane  has  128  wires,  though  we 
only  read-out  110  of  them.  The  chief  difference  between  the  two  is  the  readout 
system.  CHIPS  uses  a  delay-line  system  rather  than  DCCS. 


I'*®]  J.  McIntyre  and  G.  Warren,  Bates  Technical  Note  #94-01. 
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Figure  3.12  CHIPS  pilot  logic. 


Figure  3.13  shows  the  CHIPS  delay-line  system.  Each  wire  from  the 
CHIPS  VDCX  is  connected  to  a  MVLIOO  a/d  card.  For  the  two  separate  planes, 
every  fourth  card  is  daisy-chained  together  using  fixed-length  lemo  cables  of 
time  T  (2.2  ns),  so  that  each  delay  line  has  input  from  27  or  28  wires.  With 
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this  system  only  four  wires  per  plane  can  be  read  out,  unlike  MEPS.  Both  ends 
of  the  delay  line  are  TDC’ed  and  have  distinct  and  different  times.  The  left 
and  right  times  for  a  signal  occurring  on  anode  wire  n  out  of  a  total  number  of 
wires  N  on  the  delay  line  is: 


=  (w  -  iK  +  Id 

ti  =  {N  -  n)T  +  td  , 


(3.12) 


where  td  is  the  drift  time  associated  with  the  anode  wire  (as  was  pictured  in 
figure  3.8).  By  taking  the  sum  and  difference  of  and  td  and  n  can  be 
solved  for,  respectively,  as 


{tR+tL)-{N-l)T 


(3.13) 


and 

Using  the  wire  numbers  and  drift  times  the  particle  trajectories  can  then  be 
determined  in  the  software. 


(tH-<L)  +  (iV  +  l)r 


3.6  FPP 


This  experiment  inaugurated  a  focal  plane  polarimeter  (FPP)  on  OHIPS.  The 
FPP  was  built  by  a  collaboration  of  the  University  of  Virginia,  M.I.T.,  and 
William  &  Mary  over  about  a  three  year  period.  Figure  3.11  shows  the  FPP 
in  detail.  It  consists  of  a  graphite  analyzer,  four  multiwire  proportional  wire 
chambers  (MWPCs),  and  two  scintillator  planes  (mentioned  in  the  OHIPS 
discussion).  The  graphite  lies  above  a  pair  of  MWPCs  and  below  a  pair  of  larger 
MWPCs,  allowing  the  particle’s  track  to  be  determined  both  before  and  after 
it  scatters  in  the  analyzer.  The  FPP  was  calibrated  at  the  Indiana  University 
Cyclotron  Facility  (lUCF)  in  February  1993.  This  is  described  in  appendix  B. 
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Wire  #  in  VDCx  reference  plane 


MBD  :  Microprogrammed  Branch  Driver 
DL#  :  Delay  Line  # 

Figure  3.13  CHIPS  delay  line  system. 

The  FPP  was  designed  to  satisfy  four  basic  requirements: 

1)  Resolution  of  scattering  angle  in  graphite  to  <  1®. 

2)  Able  to  reject  small  scattering  angles  (w<  5°)  quickly  (substantially  less 
than  Ifis). 

[41]  Lourie  et  a/.,  lUCF  Scientific  and  Technical  Report,  135  (May  1992  -  April 

1993). 
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3)  Complete  azimuthal  scattering  angle  coverage  for  10%  momentum 
acceptance  and  20®  scattering  in  the  FPP. 

4)  Suitable  for  extended  targets  up  to  10  cm  in  length. 

Requirements  ^3  and  ^4  set  the  sizes  of  the  MWPCs  and  requirement  ^  1 
set  their  wire  spacing.  Requirement  ^2  is  met  by  the  FPP’s  fast,  smaU-angle 
rejection  system  described  later  in  this  section.  As  was  explained  in  section 
2.2,  this  system  was  necessary  as  approximately  95%  of  the  events  undergo 
only  smaR-angle  multiple  Cotilomb  scattering  and  are  not  useful  for  measuring 
polarizations.  This  also  means  that  the  graphite  analyzer  must  be  thin  enough 
so  that  the  multiple  Coulomb  scattering  angles  are  smaller  than  those  angles 
which  are  not  rejected.  The  graphite  analyzer  is  composed  of  slabs  of  many  sizes 
from  0.5  cm  to  10  cm  that  can  be  exchanged  and  stacked  on  top  of  each  other  so 
that  up  to  30.5  cm  of  graphite  are  available,  allowing  the  FPP  to  operate  over 
the  full  range  of  possible  Bates  proton  energies.  We  used  thicknesses  of  5  to 
10  cm.  The  FPP  support  structure  has  three  different  settings  for  supporting 
the  large  chamber  directly  above  the  graphite  so  that  the  separation  between 
the  two  top  chambers  can  be  as  large  as  possible  for  better  scattering  angle 
determination. 

3.6.1  FPP  MWPCs 


MWPCs  are  similar  to  VDCs.  The  main  difference  is  that  no  drift 
time  information  is  recorded,  just  the  presence  of  a  hit,  making  the  position 
resolution  for  individual  chambers  less  accurate.  The  FPP  MWPCs,  also,  do 
not  have  ground  wires  between  the  anode  wires  like  the  CHIPS  VDCX  does. 
There  are,  however,  ground  wires  at  the  ends  of  the  wire  planes  so  that  the 
fields  at  the  plane  edges  are  relatively  uniform.  The  chambers’  gas  mixture  was 
50%  argon  and  50%  isobutane  bubbled  through  refrigerated  isopropyl  alcohol. 
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Outflow  occurred  through  silicon  oil  to  prevent  backflow  and  to  allow  the  flow 
rate  to  be  monitored. 

The  anode  wires  are  20  /xm  gold-plated  tungsten  wire  soldered  and  epoxied 
under  tension  to  copper  artwork  on  G-10  flberglass  frames.  The  artwork  carries 
the  signals  to  tabs  on  the  exterior  of  the  chamber.  Each  chamber  has  both 
X  and  y  planes,  which  were  positioned  normal  to  the  spectrometer’s  central 
momentum  ray.  Between  and  on  both  sides  of  the  wire  planes  are  planes  of 
stretched  1/4  mil  aluminized-mylar.  We  operated  the  MWPCs  at  4.7-4.8  kV. 
As  dictated  by  the  thicknesses  of  the  G-10  frames,  the  x  and  y  planes  are  1/2” 
apart  in  the  small  chambers  and  1”  apart  in  the  large  frames.  Table  3.5  lists 
dimensions  and  parameters  of  the  FPP  MWPCs. 


Chamber 

1 

2 

3 

4 

38.0 

38.0 

93.0 

93.0 

Length  (cm) 

74.0 

74.0 

148.0 

148.0 

Wire  Spacing  (in) 

IBS 

1/12 

1/6 

1/6 

Number  of  X  Wires 

336 

336 

336 

336 

Number  of  Y  Wires 

176 

176 

208 

208 

0 

14.6 

42.2 

87.1 

CO 

1-H 

1 

13.3 

39.7 

84.5 

Table  3.5.  FPP  MWPC  Parameters,  XI  is  73.98  cm  above  the  CHIPS  VDCX 
origin  and  22.3  cm  below  the  bottom  face  of  the  graphite  analyzer.  Chamber  3  has  2 
other  possible  settings. 
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3.6.2  FPP  Electronics 

3.6.2. 1  MWPC  Readout  System 

The  FPP  MWPCs  were  read  out  much  like  the  MEPS  VDCX.  Instead  of 
DCOS,  however,  it  used  Proportional  Chamber  Operating  System  (PCOS)  III. 
Shielded  twisted-pair  (black)  ribbon  cables  took  the  signals  to  LeCroy  2735PC 
a/d  cards.  The  card  thresholds  could  be  set  via  computer  in  a  range  of  1-7.5 
volts.  We  typically  had  the  threshold  at  about  5  volts  but  had  to  set  many 
higher  to  avoid  “hot”  or  “ringing”  wires,  a  frequent  problem  that  we  were 
often  battling.  The  cards  were  plugged  into  custom-made  artwork  boards  that 
could  handle  16  cards  each.  These  boards  were  attached  to  racks,  which  were 
attached  to  the  FPP  support  structure.  The  rack  bottoms  contained  three 
cooling  fans  each.  Ribbon  cables  took  the  discriminated  signals  from  the  cards 
to  the  PCOS  crates  outside  of  the  shielding  hut.  With  so  many  cables  and  racks 
things  were  quite  crowded  in  the  shielding  hut.  The  cables  often  became  loose 
or  the  connections  did  not  line  up  properly  resulting  in  “dead”  wires,  another 
problem  we  were  frequently  fighting. 

The  signals  next  went  into  LeCroy  2731A  32-channel  delay  latch  modules 
(DCLs).  There  were  six  CAMAC  crates  containing  all  these  DCLs,  each 
controlled  by  a  LeCroy  2738  PCOS  III  system  controller  in  the  rightmost  slot. 
Planes  XI,  X2,  Yl,  and  Y2  each  had  their  own  crate,  whereas  planes  X3  and 
X4  shared  a  crate,  as  did  Y3  and  Y4.  The  DCLs  had  variable  delays  between 
300  and  600  ns  that  were  set  via  software.  Ours  was  set  at  339  ns.  This 
gives  the  coincidence  electronics  enough  time  to  generate  a  gate  defining  the 
time  interval  over  which  valid  wire  chamber  signals  may  be  recorded.  When  the 
coincidence  electronics  sends  the  controller  this  gate,  called  “El”,  the  controller 
sends  this  information  on  to  the  DCLs  which  “latch”  any  signals  that  fall  within 
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the  El  gate  and  set  their  LAM  (Look-At-Me).  It  is  important  to  note  that  by 
adjusting  the  time  of  the  El  gate  you  can  effectively  adjust  the  wire-cell  size  of 
the  MWPCs  since  a  drift  time  that  is  too  long  (or  too  short  if  you  mess  up!) 
win  not  fall  within  the  gate.  Our  El  width  was  75  ns. 

Subsequently,  the  controller  reads  the  information  in  modules  whose  LAMs 
are  set  and  repackages  the  data.  It  takes  the  controller  100  ns  to  look  for  LAMs 
and  100  ns  per  module  to  read  out  the  data.  The  controller  has  a  look-up 
memory  which  allows  logical  addresses  to  be  assigned  within  a  crate.  This 
is  particularly  important  when  more  than  one  plane  is  in  a  crate.  The  six 
controllers  are  daisy-chained  together  to  a  LeCroy  4299  databus  interface  for 
data  readout.  Figure  3.14  schematically  shows  the  PCOS  system. 

3.6.2.2  Small- Angle  Rejection  System 

The  FPP’s  MWPC  readout  electronics  were  equipped  with  a  small-angle 
rejection  system  to  quickly  reject  the  approximately  95  percent  of  events  with 
scattering  angles  too  small  to  be  useful.  Figure  3.15  gives  the  electronic  diagram 
for  this  system,  which  is  described  in  detail  in  Lourie  et  al.  Each 

PCOS  system  crate  controller  outputs  one  16-bit  data  word  for  each  hit  in  its 
plane(s).  Bits  1-9  contained  the  wire  and  slot  numbers  while  bit  11  distinguished 
between  planes  3  and  4.  The  wire  number  information  went  to  two  LeCroy 
2378  Arithmetic  Logic  Units  (ALUs).  These  units  can  calculate  sums  and 
differences  of  two  data  words.  The  ALUs  were  strobed  by  a  four-fold  coincidence 
of  the  Data  Ready  signals  from  each  of  the  three  PCOS  crates  and  bit  11  from 
crate  X3/X4,  which  meant  that  the  hit  is  in  plane  4.  (The  setup  for  Y  was 
identical.)  One  of  the  ALUs  calcidated  (plane  i  wire^  —  plane  1  wire#)  while 
the  other  calculated  (plane  4  wire#  —  (plane  1  wire#)/2).,  giving  scattering 


R.W.  Lourie  et  al,  Nuc.  Inst.  Meth.  A308,  83  (1991). 
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Figure  3.14  PCOS  III  system. 

angle  information  before  and  after  the  graphite.  (The  factor  of  2  was  due  to 
the  fact  that  the  wires  in  the  lower  planes  are  twice  as  dense  as  the  wires  in  the 
upper  planes).  These  differences  were  fed  into  a  LeCroy  2372  Memory  Lookup 
Unit  (MLU)  and  combined  to  form  the  address  of  a  single  bit  in  the  MLU. 
The  MLU  was  preprogrammed  so  that  bits  corresponding  to  good  scattering 
angles  were  set.  The  output  bit  was  AND’ed  with  the  strobe  to  handle  the 
problem  of  performing  logic  with  ECL  levels  rather  than  NIM  pulses.  Finally, 
a  coincidence  (FPASS)  checked  to  see  if  either  X  or  Y  passed  the  smaU-angle 
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test.  Figure  3.16  shows  the  carbon-scattering  angular  distribution  measured 
during  our  experiment  by  the  FPP  with  and  without  the  small-angle  rejection 
system,  showing  the  two  orders  of  magnitude  enhancement  in  the  proportion  of 
large  angle  events  using  this  system. 


FPP  •  Small  Angle  Rejection 

Jti»7,1994 


Figure  3.15  FPP  small-angle  rejection  system. 
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Figure  3.16  FPP  angular  distribution  with  (solid  line)  and  without  (dashed 
line)  the  small-angle  rejection  system. 

3. 6.2.3  Multi-hit  Rejection  System 

In  addition  to  small-angle  events,  the  FPP  electronics  were  also  able  to 
reject  events  which  had  more  than  two  hits  in  any  single  plane  (which  occurred 
about  three-fifths  of  the  time).  It  used  a  LeCroy  4532  ECL  Majority  Logic  Unit 
(MALU).  This  module  works  by  summing  the  inputs  (in  this  case  Data  Ready 
signals  from  the  various  planes)  and  seeing  if  they  pass  a  user-set  threshold. 
The  MALU  was  gated  using  a  delayed  strobe  so  that  the  MALU  woTild  not 
sum  its  inputs  until  after  the  first  two  hits  had  occurred.  Since  planes  3  and  4 
were  in  the  same  PCOS  crate  a  plane  bit  was  used  to  distinguish  between  plane 
3  and  plane  4  events,  which  were  read  out  first.  This  also  meant  that  plane  3’s 
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multi-hit  events  could  not  be  tested  for  rejection. 

3.7  Coincidence  Trigger  Electronics 

Figure  3.17  gives  a  functional  overview  of  this  experiment’s  electronics. 
Level  1  used  scintillators  to  generate  pilots,  as  was  discussed  earlier.  Level  2 
included  the  coincidence  logic  between  MEPS  and  CHIPS  and  the  generation 
of  “prescales”  (events  randomly  accepted  for  taping  regardless  of  subsequent 
trigger  logic  for  diagnostic  purposes).  In  either  case,  data  from  one  (or  both) 
of  the  spectrometers  were  “latched”  pending  level  3.  Level  3  used  the  results  of 
the  FPP  rejection-systems  described  in  the  previous  section  to  generate  either 
CAMAC  enables  and  the  event  data  CAMAC  read  signal  (PROMPTS),  or  reset 
the  spectrometers’  readout  electronics.  Inhibits  were  used  to  make  sure  that 
data  was  not  interfered  with  by  other  data  during  decision  and  readout  times. 

3.7.1  Level  2  Trigger  Electronics 

Level  2  started  with  the  pilots  OPI  and  MPI  generated  in  level  1.  (Figure 
3.18  shows  this  level’s  electronics  logic.)  These  were  AND’ed  to  form  the 
coincidence  pilot  (CPI).  OPI  was  much  narrower  (around  10  ns)  and  delayed  so 
it  determined  the  timing.  In  order  to  study  non-coincidence  electron  and  proton 
events,  and  also  coincidence  events  that  did  not  have  good  FPP  events,  prescale 
events  were  made  by  AND’ing  the  pilots  with  pulses  from  signal  generators. 
These  signal  generators  were  in  the  counting  bay  so  that  they  could  be  easily 
controlled  by  the  experimenters.  OPI,  MPI,  and  CPI  were  used  to  create  the 
prescales  OPS,  MPS,  and  CPS,  respectively;  with  the  timing  set  by  the  pilots. 
CPI  was  also  AND’ed  with  helicity  plus  and  minus  signals  from  the  control 
room,  with  the  results  scaled  so  that  any  helicity-dependence  could  be  noted. 


Each  of  the  latches  OLATCH  and  MLATCH  were  formed  in  NIM  MALUs 
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Figure  3.17  Overview  of  experiment  electronics. 

requiring  a  level  of  3.  The  MALU  inputs  for  OLATCH  were  OPS,  CPI,  and 
two  copies  of  OPI  (making  it  required);  which  set  the  timing  and  reduced  the 
output  signal’s  jitter.  MLATCH  was  made  the  same  way.  A  copy  of  MLATCH 
was  sent  back  to  MBPS  and  used  as  a  DCOS  stop.  A  copy  of  OLATCH  was 
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Figure  3.18  Level  2  trigger  electronics. 


used  as  the  El  signal  for  latching  the  FPP  wire  signals  in  PCOS.  All  8  signals 
from  this  level  needed  in  the  next  (3  pilots,  3  prescales,  and  2  latches)  were 
delayed  over  650  ns  using  a  PDM  and  ribbon  cable.  This  time  was  needed  for 
the  PCOS  readout  and  small-angle  rejection  test  to  be  completed. 
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FPP  Trigger;  Level  3 


Figure  3.19  Level  3  trigger  electronics. 
3.7.2  Level  3  Trigger  Electronics 


Figure  3.19  shows  the  third  and  last  level  of  the  trigger  electronics.  Level 
3  used  the  8  signals  from  level  2  and  FPASS.  CPI  and  FPASS  were  AND’ed  to 
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make  COINJEVENT.  This  was  wired  to  a  Schlumberger  bit  register,  as  were 
OPS,  MPS,  and  CPS  so  that  the  event  type  could  be  identified  by  which  bits 
were  set.  FPASS,  CPS,  and  OPS  were  OR’ed  in  a  PLU,  the  result  of  which 
was  then  AND’ed  in  a  subsequent  PLU  with  OLATCH  and  a  logically  inverted 
multi-hit  signal.  (In  other  words,  if  the  event  was  a  multihit  and  passed  the 
MALU,  this  AND  will  fail.)  This  made  up  the  OHIPS  CAMAC  Enable  (OCE) 
signal,  which  was  used  to  start  the  OHIPS  TDCs  and  ADC  gates.  The  AND 
of  a  logically-inverted  OCE  and  a  delayed  OLATCH  was  used  to  reset  PCOS 
system  crate  controllers  and  databus  interfaces.  A  similar  procedure  made  MCE 
and  the  MEPS  resets.  OCE  and  MCE  were  OR’ed  to  create  PROMPTS,  which 
meant  that  an  event  occurred  that  we  wished  to  read  out.  Upstairs,  an  AND  of 
MPI  and  MCE  started  a  TDC  that  was  stopped  by  an  AND  of  OPI  and  OCE. 
This  gave  us  a  coincidence  time-of-flight  (CTOF)  spectrum. 

3.7.3  Experimental  Control  Electronics 


When  the  electronics  was  deciding  whether  it  was  interested  in  an  event  or 
not,  or  when  it  was  reading  out  an  event,  new  data  needed  to  be  inhibited  so 
as  not  to  corrupt  the  old  data.  Likewise,  once  the  electronics  had  decided  that 
an  event  was  of  interest,  it  needed  a  way  to  trigger  the  data  acquisition  system 
(described  in  section  3.8)  to  readout  the  data.  These  jobs  were  handled  by  our 
experimental  control  electronics,  which  are  shown  in  figure  3.20 

There  were  two  types  of  trigger  inhibits:  Hardware  Blank  and  Frontend 
Inhibit.  The  Hardware  Blank  signal  was  used  to  inhibit  new  data  at  the  second 
and  third  levels  of  the  trigger  while  a  decision  was  made  whether  or  not  to 
read  out  an  event.  (The  first  level  was  not  inhibited  due  to  the  fact  that  since 
charge  is  counted  for  a  whole  beam  burst  the  number  of  pilot  scalers  also  needs 
to  be  counted  for  a  whole  beam  burst.)  Several  such  decisions  could  be  made 
within  a  burst  but  once  a  decision  was  made  to  record  an  event  the  Hardware 
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Figure  3.20  Experimental  control  electronics. 

Blank  would  inhibit  the  rest  of  the  beam  burst.  The  Front  end  Inhibit,  on  the 
other  hand,  inhibited  aU  three  trigger  levels  while  an  event  was  being  read  out  to 
prevent  new  data  from  corrupting  the  event.  (It  was  also  used  for  other  reasons, 
which  are  mentioned  later.)  This  typically  inhibited  the  next  two  bursts  entirely 
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while  an  event  was  being  recorded. 

The  Hardware  Blank  signal  inhibited  data  by  inhibiting  the  discriminator 
outputs.  It  was  created  from  an  OR  between  El  and  PROMPTS  signals 
that  had  been  widened  in  gate-delay  generators  (GDGs).  The  El  gate-delay 
generated  signal  was  1.1/is  long  while  PROMPTS’s  was  24/is.  The  two  signals 
overlapped.  The  El  signal  came  first  and  inhibited  new  data  while  the  PCOS 
system  was  read  out,  the  small-angle  rejection  decision  was  made,  and  resets 
were  made  (if  necessary).  If  the  coincidence  electronics  decided  to  record  the 
event,  PROMPTS  triggered  event  S  (which  did  this)  and  inhibited  the  rest  of 
the  beam  burst. 

Upstairs,  a  Timeslot  Trigger  signal  was  received  from  the  control  room 
when  beam  was  delivered.  This  signal  was  input  to  the  experimental  gate 
controller.  The  experimental  gate  controller  received  start  and  stop  signals 
from  an  output  register  that  was  controlled  by  the  data  acquisition  computer. 
It  could  also  be  disabled  (i.e.,  stopped  due  to  a  problem)  with  signals  from  the 
output  register  or  from  trip  signals  from  OHIPS  or  MEPS  wire  chambers.  The 
gun  controller  started  and  stopped  the  BT3  current  integrator. 

The  GUN  signal  was  used  with  some  gate-delay  generators  (GDGs)  to 
create  the  event  4  (scaler-read)  trigger  every  12  seconds.  Trigger  signals  were 
input  to  the  trigger  module  described  in  section  3.8.  The  trigger  module  also 
put  out  a  computer  busy  signal.  The  GUN  was  AND’ed  with  the  computer  busy 
signal.  An  OR  of  this  signal  with  the  not-GUN  created  the  Prontend  Inhibit 
signal,  which  was  sent  downstairs  to  inhibit  the  OHIPS  and  MEPS  ScintUlator 
and  OHIPS  VDCX  discriminators,  thereby  inhibiting  data  collection  when  the 
computer  was  busy  or  when  there  was  no  GUN  (i.e.,  run  and  timeslot  trigger). 
The  computer-busy  generated  by  an  event  8  (physics  event)  was  about  4  msec 
long  and  thus  inhibited  the  two  beam  bursts  subsequent  to  the  event  burst. 
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The  GUN  was  also  AND’ed  with  the  computer-not-busy  signal.  This  signal 
started  two  GDGs:  one  whose  output  gated  the  BTl  and  BT2  ADCs  and  a 
second  which  was  stopped  by  a  PROMPTS  signal.  The  output  signal  of  this 
second  GDG  was  AND’ed  with  a  pulse  generator.  The  output  of  this  AND  was 
scaled  and  served  as  a  monitor  of  when  events  occurred  within  beam  bursts. 
In  addition,  the  signal  which  started  these  two  GDGs  also  triggered  event  10, 
which  read  out  the  beam  position  monitors  and  toroid  ADCs,  as  well  as  the 
helicity  bit  and  a  timeslot  scaler  received  from  the  control  room. 

3.8  Data  Acquisition  System 

Of  the  many  crates  of  electronics  that  this  experiment  used,  only  five 
contained  modules  such  as  scalers,  TDCs,  ADCs,  and  databus  interfaces  that  we 
wished  to  read  data  from  or  write  to.  (Of  course,  additional  crates  were  reached 
through  the  databus  interfaces.)  These  five  CAMAC  crates  had  CAMAC  crate 
controllers  in  their  right-most  slots  that  were  daisy-chained  together  using 
Branch  Highway  Cable.  This  cable  was  terminated  at  one  end  and  connected 
to  a  Microprogrammable  Branch  Driver  (MBD)  at  the  other.  The  MBD  was 
connected  with  flat  cable  to  a  card  in  a  Digital  Equipment  Corp.  Micro Vax 
II  computer,  on  which  our  data  acquisition  software  ran.  The  data  acqtusition 
system  used  in  this  experiment  was  called  ‘Q’.  It  was  developed  at  LAMPF  (Los 
Alamos  Meson  Physics  Facility)  to  be  used  with  VAX’s  VMS  (Virtual  Memory 
System)  computer  operating  system  and  ran  as  a  subprocess  in  DCL  (Digital 
Command  Language).  The  two  main  tasks  of  Q  axe  the  data  aquisition  and  the 
analysis. 

Q  acquired  data  on  an  event-by-event  basis.  Events  were  triggered 
either  via  software  or  by  input  signals  to  the  trigger  module  located  in  the 
first  CAMAC  crate.  Events  0-3  and  13  were  software  triggered.  (Event 
13  was  actually  triggered  by  the  cryotarget  PC  via  a  Kinetic  Systems  3340 


Chapter  S:  Experimental  Equipment  Setup  85 


communication  interface.)  Lower  numbered  triggers  take  precedence  over  higher 
numbered  ones.  Event  0  was  triggered  at  the  end  of  a  run  to  disable  the  run 
controller.  Event  1  disabled  the  controller  if  a  run  was  suspended.  Event  2 
was  triggered  to  reset  the  electronics  and  enable  the  run  controller  if  a  run 
was  resumed.  Event  3  was  triggered  at  the  start  of  a  run  to  initialize  the 
electronics  and  enable  the  run  controller.  Event  13  dealt  with  reading  out 
target  parameters.  Event  8  was  the  physics  data  readout  event.  Event  4  read 
out  scalers.  Event  6  read  the  beam  profile  monitor  while  event  10  read  the 
beam  current  monitor  and  related  information. 

When  the  Q  subprocess  is  initialized,  a  special  QAL  (Q  Acquirition 
Language)  code  is  loaded  into  the  MBD.  The  QAL  code  defines  the  various 
CAMAC  modules  emd  teUs,  for  each  particular  event-type,  which  ones  are  to 
be  written  to,  read,  cleared,  LAM  enabled,  etc.  The  read-out  order  defines  the 
structure  of  the  raw  data  arrays.  The  raw  data  structure  for  the  main  data 
events  (8)  and  the  scalers  (4)  are  given  in  tables  3.6  and  3.7  respectively.  By 
bit  testing  the  Schlumberger  register,  readout  for  single  arm  and  prescale  events 
was  sped  up  by  setting  the  appropriate  data  words  to  zero  in  software  rather 
than  actually  reading  the  irrelevant  modules. 

Data  acquired  by  Q  was  immediately  written  to  disk.  A  batch  process 
continuously  searched  for  completed  runs  and  copied  them  (using  Q)  to  tape. 
Other  computers  were  used  for  offline  analysis.  Online  analysis  was  performed 
when  the  computer  was  not  busy  acquiring  data.  User-supplied  FORTRAN 
codes  analyzed  the  raw  data  into  useful  information  such  as  particle  kinematics 
and  measured  polarizations.  Q  allowed  the  user  to  create  histograms  of  the  data 
(raw  or  calculated),  subject  the  data  to  various  tests,  and  set-up  and  change 
parameter  values.  Figure  3.21  is  a  fiowchart  of  the  experiment’s  analysis  code. 
Details  of  the  analysis  are  described  in  the  next  chapter.  One  of  my  main 
activities  in  preparation  for  this  experiment  was  writing  and  modifying  the 


Chapter  3:  Experimental  Equipment  Setup  86 


many  different  codes  used  by  this  system  and  combining  the  many  pieces  into 
one  complete  package. 
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Element 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 

13 

14 

15 

16 

17 

18 

19-20 

21 

22-31 

32-41 

42-46 

47 

48-63 

64 

65-70 

71-80 

81 

82-87 

88 

89-90 

91 

92-100 

101 

102 

.103-127 

128 

129-130 

131-136 

137 

138-183 


Word  Description 

_ Event  Bit  Register _ 

_ Trigger  TDC  Flag _ 

OPI  TDC 
MPI  TDC 
OPS  TDC 
MPS  TDC 

_ CPI  TDC _ 

FPASS  TDC 
OLATCH  TDC 
OR  (CPS,FPASS,OPS)  TDC 
MLATCH  TDC 
OR  (CPS,FPASS,MPS)  TDC 
OLATCH  #2  TDC 
OCE  TDC 
MCE  TDC 
MLATCH  #2  TDC 

_ CPS  TDC _ 

Hardware  Blank  TDC 
Beam  Position  Monitor 
OHIPS  Delay  Line  TDC  Flag 
OHIPS  Scintillator  ADCs 
OHIPS  Scintmator  TDCs 
OHIPS  Meantimer  TDCs 
OHIPS  F2MT  OR  TDC 
OHIPS  Delay  Line  TDCs 
MEPS  Scintillator  ADC  Flag 
MEPS  Scintillator  ADCs 
MEPS  Aerogel  ADCs 
MEPS  Aerogel  Sum  ADC 
MEPS  Scintillator  TDCs 
MEPS  SO  OR  TDC 
MEPS  Scintillator  MT  TDCs 

_ CTOF  TDC _ 

MEPS  Aerogel  TDCs 
MEPS  Aerogel  Sum  TDC 
MEPS  DCOS  Flag 
25  MEPS  DCOS  Words 
FPP  TDC  Flag 
FPP  MLU  Data  Ready  TDCs 
FPP  Crate  Data  Ready  TDCs 
FPP  PCOS  Flag 
46  PCOS  Wire  Chamber  Words 


Table  3.6  Event  8  Data  Array  Structure. 
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8 


9 


10 


11 


12 


13-22 


23 


24 


25 


26 


27 


28 


29 


30-39 


40-44 


45 


46 


47 


48 


49 


50 


51 


52 


53 


54 


55 


56 


57 


58 


59 


Word  Description 


MBPS  Scintillators 


MBPS  SO  OR 


MBPS  S3  MT 


MBPS  S4  MT 


BT3 


MPI  upstairs 


blank 


MBPS  Aerogels 


MBPS  Aerogel  Sum 


Gun 


Gun  and  Computer  Not  Bus’ 


Gun  and  Computer  Not  Busy  and  H 


Gun  and  Computer  Not  Busy  and  H" 


Gim 


Gun  and  Computer  Not  Bus’ 


OHIPS  Scintillators 


CHIPS  MTs 


OHIPS  FS2  OR 


CPI 


CPS 


OPS 


OLATCH 


CPS 


MPS 


MLATCH 


MPI 


MBPS  DCOS  Sto 


PCOS  B1  Start 


CPI  and  H" 


CPI  and  H- 


CPI  and  H- 


Table  S.7.  Event  4  Scaler  Data  Array  Structure  (cont.  on  next  page).  Header 
word  is  not  included. 
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1 

Word  Description 

60 

Hardware  Blank 

61 

OPI 

62 

CPI  (1/burst) 

63 

FPASS 

64 

FPASS 

65 

CPS  (1 /burst) 

66 

67 

FPASS  1 

68 

69 

MPS  (1 /burst) 

70 

OLATCH  (1/burst) 

71 

MLATCH  (1/burst) 

72 

73 

OCE  (1 /burst) 

74 

MCE  (1/burst) 

75 

MLATCH  #2  (1/burst) 

76 

OR  (CPS, FPASS, OPS)  (1 /burst) 

77 

OR  (CPS,FPASS,MPS)  (1/burst) 

78-93 

OHIPS  Delay  Lines 

94 

PROMPTS 

95 

COINJIVENT 

96-101 

PCOS  Crate  Data  Readies 

102-103 

PCOS  MLU  Data  Readies 

104 

PCOS  SmaU-angle  Pass 

105 

PCOS  Multi-hit  Pass 

106 

OHIPS  Reset 

107 

MEPS  Reset 

108 

OPIr 

109 

OPI  2t 

110 

MPIt 

111 

MPI  2t 

112 

CPIt 

113 

CPI  2r 

Table  3.7.  Event  4  Scaler  Data  Array  Structure  (cont.). 


8821  Software 


Figure  3.21  Computer  analysis  Sowcbart, 


(90) 
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Chapter  4 
Data  Analysis 


At  the  end  of  the  previous  chapter  an  overview  and  flowchart  of  the 
experiment’s  data  analysis  code  was  given.  In  this  chapter  the  specifics  of 
the  anzdysis  will  be  discussed.  First  comes  the  extraction  of  information  from 
the  VDCXs  of  OHIPS  and  MBPS,  as  well  as  coincidence  information.  Next  the 
FPP  MWPC  analysis  is  described,  showing  how  the  second-scattering  angles  are 
determined.  Last  is  the  extraction  of  polarizations  at  the  focal  plane  and  their 
subsequent  determination  at  the  taxget,  taking  spectrometer  magnet  effects  into 
account. 

4.1  OHIPS  VDCX  Analysis 

As  was  described  in  section  3.5.2,  the  OHIPS  VDCX  delay  line  signals  from 
both  the  right  and  left  sides  can  be  added  and  subtracted  to  give  the  drift  times 
and  wire  numbers  for  an  event.  These  quantities,  modified  by  offsets,  are  shown 
in  figure  4.1.  They  are  formed  for  each  of  the  eight  delay  lines  (4  per  plane). 
While  converting  from  wire  numbers  to  wire  positions  is  quite  straightforward, 
converting  from  drift  times  to  drift  distances  is  a  little  more  involved. 

4.1.1  Determination  of  Drift  Distances 

As  shown  in  figure  3.8  and  described  in  section  3.4.2,  a  released  electron 
in  a  VDC  travels  at  a  constant  velocity  until  very  near  the  anode  wire.  Since  a 
spectrum  of  drift  times,  i.e.,  the  number  of  events  per  time  bin,  is  proportional 
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■10 


TSl  -60  11:45:45  AD0=  1:  1  FAL=  9977  LSI =  0 


ISI  0011:45:58  Ann=  I:  I  FAU  9977  LST=  0 


Figure  4.1  Raw  VDCX  spectra.  The  top  figure  shows  wire  numbers  while  the 
second  figure  shows  the  drift  times  for  the  same  delay  line. 


to  the  the  drift  velocity  (vj): 


dN  _  ^ 

dt  ds  dt  ds  ’ 


figure  4.1b  has  a  long  flat  region  formed  by  those  events  whose  tracks  passed 
through  the  constant  parallel  field  line  section  of  the  ceU.  The  spectrum’s  peak 
results  &om  those  tracks  passing  through  the  radial  field  line  area  of  the  cell 
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near  the  anode  wire.  Though  the  actual  drift  velocity  does  not  increase  much 
until  very  near  the  wire,  is  effectively  increased  due  to  the  shortest  field 
line  now  being  a  radial  one  perpendicular  to  the  inclined  track  instead  of  a 
parallel  one  near  the  edge  of  the  cell,  The  gradual  fall-off  at  the  end  of 
the  spectrum’s  flat  region  reflects  the  rising  probability  at  large  drift  times  that 
another  wire  on  the  same  delay  line  could  have  a  shorter  drift  time  and  thus  be 
the  “hit”  for  that  delay  line  instead. 


Our  method  for  converting  the  drift  times  to  drift  distances  for  the  different 
regions  is  based  on  the  work  of  Bates-ELSSY  spectrometer  users  who  spent 
much  time  investigating  this  problem,  Starting  with  a  set  of  high-quality 
drift  time  histograms,  each  histogram  is  divided  into  four  regions:  before  the 
peak,  after  the  peak  but  before  the  flat  region,  the  flat  region,  and  the  fall-off 
tail.  By  positioning  channel  1  on  the  first  significant  chaimel  of  the  peak’s 
leading  edge,  the  drift  distance  for  a  particular  drift  time  histogram  channel  is: 


drift  dist.  =  drift  dist.  of  prev.  channel  -f  Vr  •  drift  time/  histogram  channel  , 

(4.2) 

where  Vr  is  the  effective  drift  velocity  for  the  region  the  channel  is  in.  In  the  flat 
region  and  the  fall-oflf  tail,  the  effective  drift  velocity  is  just  the  terminal  drift 
velocity  (ut)  of  cm/fis.  For  the  peak’s  trailing  edge,  however,  the  eflfective 
drift  velocity  is: 


Vr  =  Vi 


^  of  counts  in  channel 


avg.  ^  of  counts/ channel  in  flat  region 

while  for  the  leading  edge  it  is  taken  to  be: 

_  #  of  counts  in  peak  channel 

avg.  of  counts/channel  in  flat  region 


(4.3) 


(4.4) 


[43]  Bertozzi  ei  aL,  Nuc.  Inst.  Meth.  141,  457  (1977). 
D.  Jordan  et  al..  Bates  Internal  Report  #92-03. 
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Figure  4.2  VDCX  drift  distance  spectra  calculated  from  the  previous  figure’s 
drift  time  spectra. 

These  drift  distance  values  for  particular  drift  time  histogram  channels  axe  then 
made  into  lookup  tables,  a  separate  table  for  each  VDC  plane.  Once  the  lookup 
table  has  been  made,  any  similar  drift  time  histogram  can  be  converted  to  drift 
distances.  Figure  4.2  shows  a  drift  distance  histogram  made  in  this  manner. 

4.1.2  Determination  of  Wire  Plane  Coordinates 


The  next  step  in  the  VDCX  analysis  is  to  convert  the  drift  times  and  wire 
positions  to  a  track  trajectory.  An  important  part  of  doing  this  is  deternaining 
which  drift  times  (which  are  all  positive)  act  as  positive  coordinates  and  which 
act  as  negative  coordinates  in  relation  to  the  wire  plane.  The  track’s  slope  and 
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wire-plane  intercept  were  determined  using  a  least  squares  linear' fit.  Due  to  the 
large  error  possible  in  very  small  drift  distances,  those  wire-hits  with  drift 
distances  less  than  1  mm  were  not  used.  AU  tracks  with  four  consecutive  hits 
or  three  hits  with  no  more  than  one  wire  gap  (which  could  have  been  because 
of  the  previously  mentioned  very  small  drift  distance)  were  fit  and  those  with 
<  0.2  were  accepted. 

The  analysis  code  was  also  designed  to  handle  an  OHIPS  VDCX  problem 
that  has  been  noticed  by  previous  experiments:  namely,  a  delay  time 

difference  jitter  that  is  on  the  order  of  one  delay  line  spacing,  causing  the  wire 
position  to  be  shifted  by  four  wires.  For  events  which  the  algorithms  could  not 
handle  or  returned  high  values,  analysis  attempts  were  made  by  shifting  each 
wire  position  by  four  in  either  direction  and  examining  if  a  better  fit  could  be 
obtained.  Given  the  adjacent  wire  and  x^  requirements,  it  is  extremely  unlikely 
that  this  coidd  produce  false  data.  This  increased  the  number  of  analyzed 
OHIPS  events  by  about  2%. 

For  an  event  to  be  considered  good  its  track  inclination  had  to  be  within 
an  acceptable  range  determined  by  the  physical  orientation  of  the  chambers. 
This  test  was  performed  by  putting  a  box  cut  on  the  2-D  histogram  of  slope  1 
vs,  slope  2. 

4.1.3  Determination  of  Focal  Plane  Coordinates 


Once  the  intercept  and  slope  are  known  for  each  of  the  two  planes,  the 
focal  plane  coordinates  can  be  calctilated.  Figure  4.3a  shows  the  orientation  of 
the  two  planes  to  each  other.  First,  the  top  and  bottom  planes  are  projected 


S.D.  Penn,  Ph.D.  dissertation,  MIT(1993),  unpublished. 
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onto  a  single  intermediate  one  via  the  relationships: 

,  mTDch 

uj.  =  -UT - 2 — 

,  ,  mBDch 

Us  =v.B  +  — 2 —  » 

where  m,-  are  the  track  slopes  (run  over  rise),  u,-  are  the  wire-plane  intercepts, 
and  Dch  is  the  separation  between  the  two  planes.  The  next  step  is  to  rotate  the 
projected  coordinates  to  a  system  with  x'  towards  the  direction  of  increasing 


momentum: 


(x' \  /  COS7  sin7\  f  u't>\  „  ,  ^ 

)  =  (  .  )  (  f  -»  =  «'  ,  (4.6) 

Vf  J  \-s1n7  COS7/  \uj<J 

as  shown  in  figure  4.3b.  Finally,  a  rotation  about  yj  to  account  for  the  VDCX’s 
tut  with  respect  to  the  central  momentum  ray  is  needed: 


(Xf\  /  cos  fl  —smCl\/x'\ 
Zf  )  Vsinfi  cos  ft  /  \z' ) 


ft  =  45.3' 


as  shown  in  figure  4.3c.  Since  z'  =  0,  a/  =  •  cos  ft .  In  the  past  Bf  and  <l)f 

were  also  calculated  using  the  slopes  and  rotation  angles  but  with  the  addition 
of  the  FPP  it  is  more  accurate  to  perform  a  linear  fit  using  the  VDCX  position 
coordinates  and  those  of  the  two  FPP  chambers  in  front  of  the  carbon  block. 


4.1.4  Determination  of  Target  Coordinates 


Determining  target  coordinates  from  the  focal  plane  coordinates  requires 
an  understanding  of  the  spectrometer’s  optics.  These  are  approximated  by  a 
second  order  TRANSPORT  matrix  that  relates  the  target  coordinates  xt  = 
to  the  focal  plane  coordinates  Xf  =  by  the 

equation 

®/  =  -^'^Tijux^x)  ,  (4.8) 

i  i.fc 

where  Mij  (Tijk)  is  the  first(second)  order  TRANSPORT  matrix.  Figure  4.4 
shows  the  two  coordinate  systems  as  they  relate  to  each  other.  Both  coordinate 
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Figure  4.3  OHIPS  VDCX  coordinate  systems.  4.3a  shows  the  orientation  of 
the  two  planes  with  respect  to  each  other.  4.3b  and  4.3c  show  the  rotations  about  J 
and  y/,  respectively. 


systems  have  i  along  the  spectrometer’s  central  momentum  ray  and  x  in  the 
dispersion  direction  in  the  magnet’s  bend  plane.  0  and  ^  are  the  angles  relative 
to  z  in  the  x  and  y  direction,  respectively.  S  is  the  relative  momentum  from 
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Transverse  Plane 


Figure  4.4  The  target  and  spectrometer  focal  plane  coordinate  systems. 


which  the  proton’s  momentum  can  be  determined: 


P  =  P,  +  SP,  , 


(4.9) 


where  Pc  is  the  central  momentum  in  the  spectrometer  (as  set  by  the  dipole’s 
magnetic  field). 


Keeping  the  dominant  terms  in  equation  4.8,  the  OHIPS  focal  plane 
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coordinates  are 
Xf  =  {x\x)xt  + 


Vf  =  Voff  +  {y\y)yt  +  {y\yS)ytOt  +  {y\yh)ytS  +  {y\<l>6)4>tS  +  {y\^6)<t>t0i 
Of  =  9off  +  {e\e)et  +  {e\0S)9t6  +  {9\8)8  +  {9\8^)8'^ 


(4.10) 


=  ^o//  +  {^\y)yt  + 

+  {4>\<l>8)<l>t8  +  {^\y8)yt8  +  {^\y9)yt9t  . 

Since  there  axe  four  equations  but  five  unknowns;  ®f ,  the  out-of-plane  coordinate 
at  the  target,  is  not  determined  but  assumed  to  be  zero  since  the  beam  spot  is 


small.  This  makes 


1/ = , 


(4.11) 


which  we  used  in  our  actual  target  coordinate  determination  equations: 


e,  -  e.„  -  («!<)<  - 

‘  (9|«)  +  {«!««)« 

“  =  {y\y)  +  {y\yS)8  +  {y\9y)9t 
b  =  {y\<l,9)9t  +  {y\<f,8)8 
c  =  {^\y)  +  {4>\9y)9t  +  {<f>\y8)8 
d={<l,\<l>)  +  {<l>\<l>8)8+{<l>\(f>9)9t 

^  <^(y/  -  yoff)  -  K^f  -  4>off) 

ad-ch 

1  ^  -  ^o//)  -  c(y/  -  yoff) 

*  ad  —  cb 


(4.12) 


The  values  for  these  matrix  elements  and  offsets  (listed  in  table  4.1)  were 
determined  by  putting  a  sieve  slit  collimator  of  known  hole  spacings  in  the  front 
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of  OHIPS  and  quasielasticaJly  scattering  electrons  through  the  spectrometer. 
Then  by  measuring  at  the  focal  plane  the  positions  of  the  peaks  associated 
with  the  different  sieve  slit  holes,  the  matrix  elements  in  the  above  equations 
could  be  adjusted  so  as  to  best  reproduce  the  physical  placement  of  the  sieve 
slit  holes.  Those  matrix  elements  involving  8  were  determined  by  taking  elastic 
data  at  different  magnet  settings,  thus  moving  the  location  of  the  elastic  pealc 
at  the  focal  plane.  (This  is  known  as  a  focal  plane  scan.)  Matrix  elements 
involving  yt  were  determined  from  runs  using  a  thin,  tilted  BeO  target  so  that 
the  y  position  was  both  known  and  variable. 


Element 

Value 

Element 

Value 

{x  6) 

4.149 

(«i<) 

10.16 

Voff 

0 

-0.0656 

{y\y) 

-1.345 

^off 

0 

{y\yS) 

-0.0533 

{<t>\y) 

-0.380 

{y\yO) 

0.00195 

WM) 

0.000538 

(ym) 

-0.000523 

-0.011768 

(y|#) 

0.00879 

m) 

-0.756 

0off 

9.527407 

(m) 

-0.00149 

m 

-0.4963 

ms) 

0.04257 

im 

-0.01529 

Table  4.1.  OHIPS  1st  and  2nd  order  optics  matrix  elements  for  a  1.60  m  first 
drift  space,  VH  focus. 
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4.2  MBPS  VDCX  Analysis 

4.2.1  Determination  of  Analyzable  Events 


Though  the  MBPS  VDCX  chamber  is  quite  similar  to  OHIPS’s,  the  read¬ 
out  system  is  very  different.  Rather  than  a  delay-line  system,  MBPS  uses  the 
DC  OS  system  in  which  drift  times  for  all  the  fired  wires  are  read.  Since  one 
is  not  limited  to  four  wires  as  in  OHIPS,  multiple  particle  trajectories  can  be 
read.  This  makes  the  analysis  more  complicated  but  also  less  corruptible. 

In  our  analysis,  MBPS  events  with  more  than  24  wire  hits  were  not 
accepted.  Since  DCOS  was  operated  in  common  stop  mode,  all  of  the  drift 
times  needed  to  be  “mirrored”  by  subtracting  them  from  a  maximum  time  and 
then  using  an  offset  to  make  the  drift  time  spectra  start  at  zero,  as  shown  in 
OHIPS  figure  4.1a.  The  drift  times  were  converted  to  drift  distances  in  the 
same  manner  as  in  OHIPS. 

The  analysis  code  next  looked  for  active  wire  clusters.  Starting  with  the 
lowest  active  wire  number,  the  first  cluster  went  until  a  second  hole  (no  active 
wire)  was  found  or  the  highest  allowed  active  wire  number  was  reached.  The 
next  cluster  then  started  at  the  next  active  wire  and  so  on.  A  maximum  of 
eight  clusters  were  accepted,  though  the  cluster  distribution  (shown  in  the  top 
of  figure  4.5)  peaked  strongly  at  one.  The  maximum  number  of  wires  in  a 
cluster  that  the  code  would  accept  was  sixteen.  The  peak  number  was  four  or 
five  (as  is  seen  in  the  bottom  of  figure  4.5).  To  be  analyzable,  a  cluster  had  to 
have  at  least  three  hits. 

For  each  analyzable  cluster,  the  wire  with  the  minimum  drift  time  was 
chosen  to  be  the  pivot  wire.  Two  situations  were  tried  to  determine  the  cluster’s 
track  (slope  and  intercept):  1)  the  drift  times  of  active  wires  below  the  pivot 
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Figure  4.5  MBPS  cluster  spectra.  The  top  figure  shows  the  number  of  clusters 
per  event  and  is  strongly  peaked  at  one.  The  second  figure  shows  the  number  of  wires 
per  cluster. 

were  taken  to  be  negative  or  2)  the  drift  times  of  active  wires  below  the  pivot 
and  the  pivot  were  taken  to  be  negative.  (The  negative  drift  times  had  to  be 
in  the  low  wires  for  the  required  positive  slopes.)  The  situation  that  gave  the 
lowest  S'  linear  fit  was  then  used. 

All  the  analyzed  clusters  were  then  compared  to  see  which  had  the  smallest 
zero  crossing  time,  i.e.,  which  required  the  smallest  timing  offset  to  make  the 
drift  distances  above  and  below  the  pivot  point  line  up.  If  both  planes  had 
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pivots,  good  clusters,  and  good  slopes  an  event  was  analyzable. 
4.2.2  Determination  of  Focal  Plane  Coordinates 


Determination  of  the  focal  plane  coordinates  worked  much  as  in  CHIPS. 
Due  to  a  different  orientation  of  the  wire  planes  with  respect  to  each  other,  the 
equivalent  of  equation  4.6  for  MEPS  is: 


sin  7 
—  cos  7 


7  =  45*^ 


As  shown  in  figure  4.3c,  d/  is  related  to  the  track  slopes  by 


{4.13) 


Bf=9-D, 

= tan~^ 

/  dx'\ 

\d^)  " 

.  —1 

•1 

duTf 

=  tan 

cos7^^,  -i-sm7^^, 

COS  7  sin  7 

=  tan 

— -  +  — - 

-D 

rriB  rriT 

=  tan“^  1 

rriT  COS7  +  ttib  sin7l 

TTIb^TIt 


(4.14) 


while  is  given  by 


tan 


-1 


=  tan  ^ 


dyfjd^ 

dzf/dz^ 


=  tan  ^ 
=  tan“^ 


sin  'y  _  cos  7  -i 

Tns 

cosD  +  (^  +  ^)sinn 

*■  V  rriB  wir  ' 

rriT  sin  7  —  ms  cos  7 

msmT  cos  fl  +  (mr  cos  7  +  ms  sin  7)  sin  D 


(4.16) 
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4.2.3  Determination  of  Target  Coordinates 


For  MBPS  we  used  reverse  matrix  elements  as  given  by  the  inverse  of  eqn. 
4.8: 

x\  =  >  (4*16) 

i  j,k 

so  that  the  target  coordinates  were  determined  by 

Sf  =  0 

St  =  (^1®)®/  +  {S\x^)x'f 

Ot  =  Ooff  +  {0\x)xf  +  {e\e)0f  ^.{0\x^)x)  +  {0\0x)0fXj  +  {0\y^)y)  +  {0\0^)0) 

<Pt  =  <f>off  +  (^1®)®/  +  {^\y)yf  +  +  (^1®^)®/ 

+  {<l>\xy)xfyf  +  {<f>\04>)0f(f>f  +  {^\0y)0fyf 

yt  =  yoff  +  {y\x)xf  +  {y\^)<l>f  +  {y\y)yf  +  {y\4>'^)4>}  • 

(4.17) 


Table  4.2  gives  the  relevant  MEPS  reverse  matrix  elements. 

4.3  Scattering  Angles 

MEPS  and  OHIPS  both  detect  particles  in  the  same  plane.  This  plane  is 
horizontal  and  at  the  height  where  the  beam  hits  the  target.  From  fig.  4.4  it 
should  be  clear  that  <f>i  is  in-plane  while  Of  is  out-of-plane.  Therefore,  the  in- 
and  out-of-plane  proton  scattering  angles  (denoted  as  horizontal  and  vertical, 
respectively)  are  defined  as: 

nh  _  nOHIPS  •  OHIPS 
“  ^hall  ~  9t 

nv  _  nOHIPS 


(4.18) 
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Element 

Value 

Element 

Value 

(^l®> 

0.5521 

(m) 

-0.34877 

0.000777 

(^|y) 

-10.756 

^off 

-9.7003 

(^|®2) 

-0.0153835 

m 

0.26617 

{(f>\xy) 

0.19069 

m 

-0.38908 

{m) 

-0.0021425 

{e\x^) 

-0.02387 

{4>\6y) 

0.010918 

{e\ex) 

0.0019076 

Voff 

3.1830 

mn 

0.0003005 

(y|®) 

-0.24770 

W) 

0.091117 

{y\<l>) 

0.83792 

4>off 

0 

{y\y) 

-6.3145 

(^1®) 

0.10358 

{y\<f>^) 

-0.0025595 

Table  4.2.  MBPS  reverse  matrix  elements. 


and  the  corresponding  electron  angles  as: 


gv  _  _ qMEPS 


(4.19) 


and  dKS'-*  are  the  positive  in-plane  angles  between  the  beam-dump 
pipe  {6  =  0)  and  the  spectrometers.  The  electron’s  actual  scattering  angle  is 
thus: 


6e  —  cos  ^(cosflg  COS^g)  . 


(4.20) 
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4.4  Coincidence  Spectra 

4.4.1  Time-of-Flight  Spectra 

Any  coincidence  experiment  is  going  to  have  both  true  coincidences  and 
accidentals.  The  rate  of  these  accidentals  is  proportional  to  the  individual 
spectrometer’s  singles  rates  and  the  coincidence  resolving  time  (A^t): 

Ra  =  ,  (4.21) 

where  d.f.  is  the  accelerator’s  duty  factor.  Proof  of  the  occurrence  of  true 
coincidence  events  is  provided  by  the  coincidence  time-of-flight  (CTOF)  TDC 
spectra.  The  accidentals  that  get  into  this  spectrum  will  have,  in  general,  a 
flat  time  distribution  while  the  time-correlated  coincidence  events  will  result  in 
a  narrow  spike  on  top  of  the  flat  background.  An  example  of  this  spectrum 
for  this  experiment  is  shown  in  figure  4.6.  For  experiments  with  appreciable 
background,  one  normally  has  to  correct  this  spectrum  for  pathlength  and 
velocity  variations  in  order  to  improve  the  signal  to  noise  ratio.  Such  was 
not  the  case  for  us,  given  our  kinematics  and  simple  hydrogen  target. 

4.4.2  Missing  Energy  Spectra 

Two  quantities  that  are  useful  for  determining  if  the  coincidence  kinematics 
axe  properly  understood  are  the  missing  energy  and  momentum.  For  elastic 
scattering  from  hydrogen,  these  should  both  be  zero.  The  missing  energy  is 
equal  to: 

Em=t>}-Tp  =  E-E'  -Tp  .  (4.22) 

The  missing  momentum,  also  referred  to  as  the  recoil  momentum,  is  given  by: 

Pr  =  k  —  it'  —  p'  .  (4.23) 
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Figure  4.6  CTOF  TDC  spectrum  with  no  corrections. 

In  order  for  these  spectra  to  peak  at  zero,  the  input  quantities  have  to  reflect 
what  happened  at  the  scattering  point  inside  the  target.  This  means  that 
corrections  have  to  be  made  for  energy  losses  in  the  target. 


There  are  two  basic  types  of  electron  energy  losses  in  targets:  radiation  and 
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ionization  losses.  At  our  experiment’s  energies,  radiation  losses  cause  a  larger 
average  energy  loss  than  ionization  losses  in  hydrogen.  However,  the  energy 
loss  probability  goes  like  l/E^  so  that  the  most  probable  energy  loss  is  zero. 
The  most  probable  ionization  losses,  however,  are  nonzero  so  that: 

E  =  Eo-  eion  .  (4.24) 


The  most  probable  energy  loss  due  to  ionizing  atomic  electrons  (Landau 
straggling)  is  given  by: 


Cton  =  ^  In  ^  -I-  0.37 

(4.25) 

where 

^  2irNaeHp  Z 
meC^  A  ’ 

(4.26) 

where  p  is  the  material  density  and  t  is  its  thickness,  and 

2me 

(4.27) 

I  =  Z  •  1.35  X  10“®  MeV  and  is  the  average  atomic  ionization  potential.  Losses 
in  the  target  walls  were  insignificant.  Roughly  speaking,  the  incoming  and 
scattered  electron  both  lost  about  0.9  MeV  from  ionization  losses. 


The  proton’s  energy  losses  were  calculated  using  RASP  -  a  program 
based  on  nuclear  data  tables,  The  energy  losses  were  dominated  by 

those  occurring  in  the  hydrogen  itself,  rather  than  in  the  target  walls.  Total 
losses  in  the  target  were  approximately  2  MeV. 

I'*®!  H.  Uberall,  Electron  Scattering  from  Complex  Nuclei,  Vol.  B,  Academic  Press  (1971). 

C.S.  Whisnant,  private  communication. 

[48]  Barkas  and  M.  Berger,  Tables  of  Energy  Losses  and  Ranges  of  Heavy  Charged 

Panicles,  NASA-SP-3013  (1964). 

J,  Janni,  At.  Data  and  Nucl.  Data  Tables  27,  nos.  4-5  (1982). 
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Figure  4.7 


Missing  energy  (top)  and  momentum  (bottom)  spectra. 


With  the  energy  losses  calculated,  the  true  momenta  seen  by  the  spectrom¬ 
eters  were  known.  Using  the  measured  relative  momentum  S,  the  actual  centr^ 
momentum  of  the  spectrometers  could  then  be  determined  using  equation  4.9. 
Figure  4.7  shows  a  typical  missing  energy  (top)  and  missing  momentum  (bot¬ 
tom)  from  our  experiment.  The  missing  energy  is  accounted  for  to  within  1  MeV 
and  the  missing  momenta  peaks  near  zero,  as  expected.  (The  reason  it  does 
not  peak  at  zero  is  that  its  yield  integral,  in  addition  to  containing  a  peak  at 
zero  widened  due  to  resolution,  multiple  scattering,  etc.;  contains  a  phase  space 
term  p^dp  that  goes  to  zero  as  p  goes  to  zero.)  The  tails  in  both  spectra  are 
due  to  radiative  losses  by  the  electron,  both  before  and  after  the  scattering. 
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4.5  FPP  MWPC  Analysis 

Our  largest  analysis  effort  was  devoted  to  polarization  determination.  The 
information  from  each  of  the  four  MWPCs  had  to  be  combined  to  determine 
the  particle  tracks  before  and  after  the  carbon  and  then  the  scattering  angles 
in  the  carbon.  Whereas  our  MBPS  and  OHIPS  analysis  codes  were  based  on 
the  works  of  previous  experimenters,  our  FPP  analysis  code  was  entirely  our 
own  creation.  Its  description  makes  up  this  section. 

4.5.1  Determination  of  Carbon  Scattering  Angles 


The  first  step  in  the  analysis  is  to  convert  the  raw  PCOS  data  into  wire 
numbers.  The  wire  numbers  were  converted  into  wire  positions  with  the  origin 
defined  as  the  middle  of  the  chamber.  Hit  wires  in  the  same  plane  that  were 
adjacent  were  clusterized  into  a  “single”  hit.  Unlike  drift  chambers,  position 
resolution  in  MWPCs  cannot  be  better  than  a  wire  spacing.  The  chambers  were 
aligned,  both  internally  and  with  respect  to  the  VDCX,  using  software  offsets. 
The  determination  and  consequences  of  these  offsets  are  described  in  appendix 
C. 


In  order  for  track  reconstruction  to  be  possible,  there  had  to  be  one  and 
only  one  hit  in  each  plane.  If  one  of  the  planes  had  two  nonconsecutive  hits 
and  the  other  seven  had  single  hits,  analysis  was  performed  using  both  cases,  in 
the  hope  that  cuts  performed  after  the  trajectories  had  been  determined  would 
eliminate  one  of  them. 

Figure  4.8  illustrates  the  track-reconstruction  problem.  From  the  observed 
chamber  event  positions  (represented  by  Ac,  where  c  is  the  chamber  number  1 
through  4)  and  the  known  distances  between  the  chambers  one  must  determine 
the  initial  and  final  trajectories  (5,-  and  Sf),  and  then  the  scattering  angles 
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from  those.  The  polarimeter’s  coordinate  system  is  defined  like  OHIPS’s  to 
be  right-handed  with  z  along  the  central  ray  normal  to  the  chambers  and  x 
in  the  direction  of  the  increasing  momentum  dispersion.  The  direction  i{f)  of 
the  proton’s  initial(final)  trajectory  is  determined  entirely  from  the  positional 
information  from  the  front(rear)  chambers: 


i  oc  ( Aa52i ,  A2/21 ,  A2:2i  ) 

=  {^X2l/Az21,Ay2l/Az21,l) 

=  (tan  d,-,tan  1) 


(4.28) 


f  Oc(Ax43,  A2/439  AZ43) 
=  (tan  6f,taja.  1) 


(4.29) 


The  trajectories  Si  and  5/  are  simply  straight  lines  of  slope  i  and  /: 


Si  —  Ai  -j-  iti 

Sf  =  j4j  +  ftf  , 


(4.30) 


where  Ai  is  5j’s  intercept  in  plane  1  and  is  5/’s  projected  intercept  in  plane 
1.  To  solve  for  A\  note  that 


Sf{z^)-z  =  '  z  +  f  -  Ztf{z4) 

Sf{z^)  •z  =  f-  ztf{zi) 

Z4  =  tf{Zi)  , 


(4.31) 


A'y  =Sf{Zi)  -  fZ4 

—s  — ♦ 

=-4.4  —  /Z4  . 


(4.32) 


Once  the  basic  trajectory  vectors  are  set  up  the  scattering  angles  can  be 
determined.  The  polar  scattering  angle  is  simply  the  angle  between  the  two 
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proton  trajectories: 


6  —  cos 


-1  »•/ 


(4.33) 


Determination  of  the  azimuthal  scattering  angle  <f>  is  more  complicated.  First, 
let’s  define  an  event  coordinate  system  in  which  z'  lies  along  the  proton’s 
momentum r.  Each  event’s  z'  will  be  different,  though  the  variation  will  be  small 
since  the  vertical  angle  0foc  is  small  (<~  11°).  These  primed  coordinates  must 
transform  to  the  spectrometer’s  when  the  vertical  angle  goes  to  zero  (i.e.,  the 
proton  goes  straight  up).  By  definition,  <l>  equals  0  along  the  x'  axis.  Expressing 
eq.  4.28  as  t  =  (i®, »»,!),  the  event  unit  vectors  are: 


y/a 


x'  =y  X  z'/\y  X  z'\  = 


l  +  i| 


^(IjO,  ix) 


(4.34) 


Af  a/ 

y  =z  X  X 


(— iutj,  1  +  ij.,  iy) 


y/a{l+iiy  "  " 

where  a  =  i  -  i.  In  these  coordinates,  the  relationships 

f  .  x'  =1/1  sin  d  cos  ^ 

.  ,  (4.35) 

f  -y  =1/1  sin d sin ^  , 

lead  to  the  azimuthal  scattering  angle  . 

f  ■  *'  (4.36) 

_ -fxixiy  +  /p(l  +  il)  -  iy 

y/a{f:,-ix) 

(At  this  point  some  computational  advice  is  in  order.  Remember  that  the 
arctan  function  will  only  return  values  between  —f  and  f  [quadrants  I  and 
IV].  The  two  other  quadrants  [II  and  III],  however,  have  negative  cos<f>  [the 
previous  equation’s  denominator]  values  instead  of  positive  ones.  Of  these  two 
quadrants,  quadrant  II  has  a  positive  sin^  [the  previous  equation’s  numerator] 


(4.36) 
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value  and  quadrant  III  has  a  negative  one.  Thus  a  correct  four-quadrant  ^ 
value  can  be  determined.) 


In  addition  to  the  scattering  angles  we  also  used  the  before-and-after 
trajectories  to  determine  the  distance  of  closest  approach  of  the  two  tracks. 
We  define  the  vector  A  as  the  difference  between  the  two  trajectories: 

A  =  Si-Sf  =  AA  +  iti-ftf  ,  (4.37) 

where  AA  =  Ai  —  A'^.  To  get  the  vector  coordinates  at  the  intersection  point 
(more  correctly,  the  point  of  closest  approach  of  the  two  trajectory  vectors), 
A  needs  to  be  minimized  with  respect  to  t,-  and  tf.  It  is  actually  easier  to 
minimize 

A^  =  (AA)^  +  aif  +  ^t^f  —  2‘ytitf  —  26tf +  2€ti  ,  (4.38) 

where 

a  =i-i 


p  =/ •  / 

7=?-/  (4.39) 

S  =  -f‘AA 
e  =  —  i  •  Ai4  . 


Performing  the  minimization  gives: 

S  , 

— A^  =2ati  —  2'ytf  —  2e  =  0 

A^  =2/St f  -  2‘yti  -26  =  0 

otf 


(4.40) 


Solving  these  two  equations  simultaneously  gives  ti  and  tf  at  the  moment  of 
closest  approach: 

.c  ^  iS-pe 

‘  7=  -  a/3 

^  aS-^e 

^  7^  — 


(4.41) 
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Using  eq.  4.37  the  distance  of  closest  approach  dc  between  the  two  reconstructed 
trajectories  can  be  determined: 

dx  =AAx  +  ixtf  -  /*</ 


dy  =AAy  +  iyt^  -  fy^ 

dz  =<“  — 


(4.42) 


dc  =V^  ■ 

Given  the  wire  spacings  (which  define  the  chsumber  resolutions  for  MWPCs) 
and  plane  distances  in  table  3.5,  the  maximum  error  in  the  determined  location 
of  either  the  front  or  back  track  in  the  center  of  the  carbon  block  is  about  3  mm. 
Not  taking  dz  into  account,  this  implies  that  the  distance  of  closest  approach 
should  generally  be  less  than  6  mm.  The  top  histogram  in  figure  4.9  shows  our 
dc  spectra,  which  indeed  falls  in  line  with  this  estimate  of  the  resolution. 


It’s  also  useful  to  check  that  the  trajectory  intersection  occurred  in  the 
carbon  slab.  The  coordinates  of  the  closest  approach  point  can  be  determined 
by  averaging  the  two  trajectory  vectors  at  the  point  of  closest  approach: 

f  =  s/2  =(5i((a  -  SXl»)/2 

=(i4i  +  i4j)/2  +  (iij  + /<^)/2  .  (4*43) 

={Ai  +  J4  —  fzi)l2  +  {it\  +  ft^f)!^ 

Using  this  and  the  two  trajectories,  the  pathlength  the  proton  travels  in  the 
carbon  before  scattering  can  also  be  determined.  The  bottom  histogram  in 
figure  4.8  exhibits  our  pathlength  spectra,  which  is  basically  the  width  of  the 
carbon  block  (7  cm  in  this  case). 

Once  the  track  (or  possibly  two  tracks)  had  been  reconstructed,  it  was 
tested  to  see  if  it  had  acceptable  values  for  the  x  and  y  coordinates  of  the  closest 
approach  (c),  the  initial  trajectory  angles  (d,-  and  and  the  polar  scattering 
angle  (0).  These  values  were  determined  by  stud3dng  events  that  only  had  one 
hit  (or  cluster)  per  plane.  Hopefully,  a  two-track  event  had  one  track  eliminated 
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Figure  4.9  Spectra  showing  the  FPP  MWPC^s  ability  to  determine  the 
scattering  point  from  the  measured  trajectories.  The  top  figure  shows  the  closest 
distance  between  the  two  trajectories  and  the  bottom  figure  shows  the  calculated 
distance  of  the  scattering  vertex  from  the  front  of  the  carbon  block. 


by  these  cuts  so  that  it  could  be  used.  At  this  stage  the  remaiiung  events  were 
considered  to  be  correctly  reconstructed  and  analyzable.  Using  this  method, 
we  were  able  to  analyze  approximately  55%  of  the  FPP  events  instead  of  the 
45%  that  were  single-cluster  in  each  plane. 

Since  there  is  little  information  to  be  gleaned  from  small  angle  events  and 
the  analyzing  powers  {Ac{0))  above  20  degrees  are  not  well  known,  we  next 
limited  our  “good”  events  to  those  with  scattering  angles  between  7  and  20 
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CARBON  DETECTOR 

PLATE 


Figure  4.10  Diagram  of  a  cone-test. 


degrees.  (See  appendix  C  for  a  discussion  of  the  lower  veilue.)  Furthermore, 
we  also  required  events  to  pass  a  cone-test,  shown  in  figure  4.10.  The 

difference  between  the  location  of  an  actual  event  in  the  reeir  detector  (plane  4) 
and  where  it  would  have  been  had  the  proton  not  scattered  serves  as  a  radius 
of  a  cone’s  base  for  all  possible  ^  scatterings  from  a  given  point  with  a  given 
initial  trajectory.  If  part  of  this  cone  lies  off  plane  4,  such  as  for  track  B  in 
the  figure,  this  event  must  be  thrown  out  to  avoid  making  a  biased  sample  of 
the  scattering  data.  In  our  case,  due  to  the  size  of  our  detectors,  such  events 
occurred  rarely.  An  event  passing  this  final  test  was  considered  to  be  a  good 
event  useful  for  extracting  physics  measurements. 

4.5.2  Determination  of  Spin  Precession  Angle 

Once  it  was  decided  that  an  event  was  good  it  was  also  necessary  to  determine 
its  spin  precession  angle  x  in  the  spectrometer.  The  proton’s  initial  directional 


J.A.  Edgington,  Nuc.  Inst.  Meth.  164,  175  (1979). 
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vector  at  the  entrance  to  the  spectrometer  is  given  by 


^tan  Ofgf^  tan  1) 


while  its  final  directional  vector  at  the  focal  plane  is 

Pfp  =  ^  (tang/p,tan^/p,l) 

y'l  +  tan^5/p  +  tan^^/p 


(4.44) 


(4.45) 


We  wish  to  know  the  final  vector  as  expressed  in  the  initial  target  coordinate 
system.  This  is  simply  a  coordinate  transformation: 


p;\  /o  0  -1 

P}\=  0  1  0 

Pi)  Vl  0  0 


(4.46) 


Since  the  bend  is  only  in  the  xz  plane  the  bend  angle  is  given  by: 

^hend  —  COS  (J^i  *  P f) 

From  eqn.  3.5,  then,  the  spin  precession  angle  is: 


(4.47) 


X  =  ^6enrf7(|  -  1)  > 


(4.48) 


where  7  =  +  1. 


4.6  Extraction  of  Polarizations 


The  last  section  described  the  event  selection  and  trajectory  determination 
procedures.  The  interesting  “physics”,  however,  is  not  in  the  individual  events 
but  in  the  overall  collection  of  events.  Recall  from  chapter  2  that  our  purpose 
for  using  the  FPP  was  to  measure  the  outgoing  proton  polarization  in  an  (e,e'p ) 
reaction.  This  involves  measuring  thep-^^C  angular  distribution  (eq.  2.4)  from 
the  scattering  that  occurs  in  the  FPP’s  carbon: 

I{6, 4)  =  h{9)  [1  -  p{^Ac{9)  sin  ^  +  p^^Ac{0)  cos  .  (4.49) 
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or,  in  terms  of  the  measured  asymmetries  ea(fl)  and 

I{0, 4>)  =  h{6)  [l  +  ej,  cos  ^  +  e,.  sin  .  (4.50) 

The  polarizations  are  just 

fp  _  ^z(g)(^) 

Ac{0) 

(with  the  proper  sign,  of  course).  Recall  also  from  chapter  2  that  the  focal  plane 
polarizations  are  related  to  the  initial  target  polarizations  by  the  spectrometer’s 
spin  precession  angle  x* 


Equation  4.50  includes  physical  asymmetries  but  not  intrinsic  systematic 
ones.  Their  effects  on  the  polarimeter’s  azimuthal  response  can  be  represented 
by  a  FoTirier  series: 

oo  oo 

=  1+5^  em(^)cos(m^)  +  ^  em(^)sin(m^)  ,  (4.52) 

m=l  m=l 

where  ^(0,  </>)  is  normalized  to  one  in  the  absence  of  any  systematic  azimuthal 
asymmetry.  Thus  the  actual  measured  distribution  is  the  physical 

distribution  modified  by  the  systematic  effects: 


N{e,<j>)  =  i{eA)  i{6A) 

=  h{9)R{e,<t>)  , 


(4.53) 


where 


R{6^  ^)  =  [  1  +  (cy  cos  ^  +  e*  sin  ^)  ] 

OO  OO 

X  [  1  +  ^  em{0)  cos(m<^)  +  em{6)  sin(m</>)  ]  . 

m=l  m=l 


(4.54) 


D.H.  BarkhufF,  private  communication. 
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Using  simple  trigonometry  addition  formulas,  this  is: 

R{6f  ^)  =  1  +  ^  em(9)  cos(m^)  +  ^  €^(6)  em{m<f>) 

+  ej,|cos  ^  em{9)[  +  cos(m  -  1)^  +  cos(m  +  1)^  ] 

+  ^  ]} 

+  ^  ^  -  sin(m  -  1)<^  +  sin(m  +  1)^  ] 

+  ^  +cos(m  -  1)^  -  cos(m  +  1)^  ]|  . 

(4.55) 

Many  different  sums  involving  good  events  were  created  so  that  overall 
averages  could  be  determined.  For  example,  sums  of  the  spin  precession  in  the 
spectrometer  %:  S®^^X  were  needed.  Various  kinematic 

ev  ev  ev 

factor  sums  were  needed  for  the  determination  of  Ge /Gm  and  the  beam  helicity. 
Most  importantly,  the  good  events  were  summed  five  different  ways  by  (helicity 
state)  X  (^-bin).  These  five  different  sums  were:  1,  X)  ^  sin 

ev  ev  ev  ev 

and  '^sid.2^.  These  were  used  to  determine  the  focal  plane  asymmetries 

ev 

(real  and  systematic),  using  the  Fourier  analysis  method  developed  by  fellow 
graduate  student  David  Barkhuff  and  described  in  the  following  paragraphs. 


The  asymmetries  ex(j,)  of  eq.  4.50,  assuming  full  azimuthal  acceptance,  are 
determined  by  integrals  of  the  distribution  weighted  by  cos^,  sin^,  and  1: 


i*27r  i*27r 

/  d<^  1(0,  <^)  =  2  /  #  1(0,  (t>)  cos  (f> 

Jo  Jo 

f2n  M2n 

/  d(f>  1(0,  ^)  =  2  /  d<t>  1(0,  <l>)  sin  <f> 
Jo  Jo 


(4.56) 


For  the  situation  where  intrinsic  systematic  asymmetries  are  present  the 


D.  Besset  et  al,,  Nuc.  Inst.  Meth.  166,  515  (1979). 
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weighted  integrals  Uj{0)  and  Vj{6)  are  defined  as: 

i»27r  f2n 

Uj  I  d(}>  N{9,4>)  =  /  d<f>N{0,<f>)  cos{j<f>) 

Jo  Jo 

A2ir 

vj  /  d<t>N{e,<i>)  = 

Jo 

These  integrals  axe  estimated  by  sums  over  the  observed  events  so  that 

“i  “  =  cos(j^) 

_  (4.58) 

~  V^l5^sin(j^)  =  sin(j»  , 

ev 

where  =  5^  1* 

ev 

It  should  be  pointed  out  that,  although  we  initially  binned  our  sums  by  0 
in  addition  to  helicity  and  we  ultimately  determined  that  it  did  not  affect 
the  final  result.  Thus  the  measured  asjrmmetries  e*  and  ty  and  the  weighted 
integrals  Uj  and  Vj  were  integrated  over  6.  Also,  the  analyzing  power  used  in 
equation  4.49  was  an  average  analyzing  power  Ac  rather  than  a  ^-dependent 
one.  The  average  came  from  a  sum  of  analyzing  powers  calculated  for  each 
good  event  using  that  event’s  0  and  kinetic  energy.  The  analyzing  power  used 
was  an  average  of  four  different  models.  The  four  models  were  of  two  slightly 
different  model  forms  (the  SIN  and  LAMPF  forms),  fit  both  with  and 
without  our  lUCF  200  MeV  data.  All  four  fits  used  the  data  from  SIN  and 
LAMPF.  The  standard  deviation  around  the  mean  of  the  models  was  a  third 
of  one  percent;  better  than  the  overall  2%  uncertainty  of  the  models. 

If  there  are  no  systematic  asymmetries,  evaluation  of  the  eq.  4.57  integrals 
shows  the  only  nonvanishing  averages  to  be  u\  and  vi,  which  equal  ej,/2  and 
ea,/2,  respectively.  In  this  case  and  vi  are  directly  proportional  to  the 
physical  asymmetries.  Conversely,  if  there  are  no  physical  asymmetries,  Uj  and 
Vj  equal  ej/2  and  ej/2,  making  them  directly  proportional  to  the  systematic 


I  d<f>N{0,<f>)  sin(j^) 
Jo 


(4.57) 
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(4.59) 


asymmetries.  More  generally,  however,  when  there  are  both  physical  and 
systematic  asymmetries,  one  gets: 

By  (2  +  62)  +  e*  62  +  2ei 

til  =  - 

4  +  2e„  61  +  2e,  61 

_  c®  (2  —  62)  +  Cj  62  +  26i 

4  +  2ej,  61  +  2ej.  ei 

This  becomes  a  pair  of  inhomogeneous  linear  equations: 

ttyCp  +  l^yCx  —  6y  =0 

(4.60) 

“1“  hx  “0  , 

where  the  coefficients  are  given  by 


(4.60) 


a»  =  1  +  2®2  -«iei 

=  ^C2  -  uiei 

Sy  =  2ui  —  61 


62  —  •W161 


/3x=l-  ■^e2-viei 


2t>i  —  61  . 


(4.61) 


Thus  the  observed  sums,  the  physical  asymmetries,  and  the  systematic 
asymmetries  are  related  by  the  matrix  equation: 


W-^e  =  i  , 


(4.62) 


where 


wnerc 

y  ■■(::)  -(y  ■  <“> 

From  this  equation  it  is  clear  that  the  physical  and  systematic  asymmetries 
cannot  be  simultaneously  determined.  However,  when  unpolarized  electrons 
scatter  on  protons  there  are  no  physical  asymmetries  so  the  systematic 
asymmetries  ej  and  ej  can  be  determined  as  described  earlier.  These  can  then 
be  used  to  extract  the  physical  asymmetries  &om  the  data. 


Equation  4.62  can  be  rewritten  as 


\e®/  \  UlCl  -  ^62  1  +  062  -  WlCl  /  V2vi  -  61  / 
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where  the  determinant  is  given  by: 

\W\~^  =  +  §2  —  61  e2)+-^(ei 62+6162)  — -(e2  +  e|)  .  (4.65) 

If  the  systematic  asymmetries  are  small,  a  good  estimate  of  the  statistical 
uncertainty  in  the  physical  asymmetries  can  be  obtained  by  keeping  only  the 
first  order  terms  in  Cj  and  ej.  Then  the  determinant  is  approximated  by: 

\W\(=il+viei+Uiei  ,  (4.66) 


SO  that 

/  ej,  \  ^  / 1  -  |e2  +  uiei  ttici  -  §62  \  / 2ui 

\Cj!/  V  U161  -  |e2  1  +  |62  +  viei  /  \2i;i 

2tii(l  +  uiBi  +  wiei)  —  61  —  ^162  —  ^162 
2wi(l  +  ^161  +  viei)  —  61  —  ^162  —  W162 

(to  first  order).  Provided  they  are  small  compared  to  those  in  Ui  and  vi,  one 
can  ignore  the  statistical  uncertainties  in  ej  and  ej.  By  rewriting  the  event  sums 
of  eq.  4.58  as  sums  of  different  <f>  bins: 

N 

cos(i^)  =  ^  n,-  cos(jV,)  ,  (4.68) 

€V  i 

where  n,-  is  the  number  of  events  in  the  tth  of  N  phi  bins  so  that  H  =  5^,-  n,-, 
the  statistical  imcertainties  are  given  by: 

(<'l\ ^ 

Note  that  =  n,-.  Performing  the  differentiations  and  using  eq.  4.58,  the 
statistical  uncertainties  are: 


_  /(l  +  2uiei +t;iei  -  |e2)^  («i6i  -  §62)^  \ 

<rl)  V  (vi6i-|e2)2  (1+uiei +2t;iei  -  |e2)V  , 

2fl  +  U2-2ul\ 

\  1  —  142  —  2-1*1  / 
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Using  these  results  the  average  focal  plane  polarizations  for  each  helicity 
state  can  be  determined  (remember  the  minus  sign  for  pl^): 


—Ip 

~Ac 


(4.70) 


iK'i) 

<'(?()  =  ^  • 


(4.71) 


4.7  Quadrupole  Depolarization  Effects 

The  method  of  determining  the  target  polarizations  and  GeIGm  using  the 
helicity  states  and  the  spin  precession  in  OHIPS’s  dipole  was  described  in  section 
2.2.  Precession  also  occurs  in  the  magnetic  fields  of  the  quadrupoles,  however, 
and  failure  to  account  for  this  results  in  the  determined  target  polarizations 
being  lower  than  they  should,  or  “depolarized”.  The  polarizations  are  lowered 
because  the  polarization  is  in  actuality  an  average  over  many  different  spin 
rotation  directions  2aid  failure  to  account  for  any  effect  that  randomizes  the 
spin  rotation  directions  will  lower  this  average.  Quantifying  this  effect  is 
more  difiicult  for  quadrupoles  because,  unlike  dipoles,  their  magnetic  fields  are 
inhomogeneous.  Our  solution  to  this  problem  is  based  on  the  work  of  Nurushev. 

[53] 


The  Thomas  equation  (3.1)  reduces  to 

dp  [ff^  /  g -2 


=  —pxA 


7 )  Dt 


(4.72) 


in  the  absence  of  an  electric  field.  Bl  and  Bt  are  the  magnetic  field  components 
parallel  and  perpendicular  to  the  particle’s  velocity,  respectively.  If  the  particle 


S.  Nurushev,  Nuc.  Inst.  Meth.  141,  417  (1977). 
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(in  this  case  a  proton)  moves  in  the  ii-direction  with  a  constant  velocity  the 
polarization  components  are: 

Integrating  this  gives: 

Pk{z)  =  pjfe(O)  +  —Cijk  f  PtiOa-iiO^C  •  (4-74) 

pc  Jo 


Assuming  the  fields  in  a  quadrupole  can  be  modeled  as  Bl  =  0  and 
Bt  =  G{yx  +  xy),  the  equations  of  motion  for  a  particle  in  a  quadrupole  are: 


x"  ±K^x  =  0 
y"  T  K^y  =  0 


(4.75) 


where  ^  and  the  upper(lower)  signs  refer  to  vertical  and  horizontal 

focusing  quadrupoles,  respectively.  (Recall  the  TRANSPORT  coordinates  of 
figure  4.4.)  This  results  in  eq.  4.74  becoming: 

Pt{z)  =  P<(0)  -  A  /*  piiC)y(C)dC 
Jo 

Pn{z)  =  P„(0)  +  A  /  pj(C)®(C)«^C  (4-76) 

Jo 

pi{z)  =  pj(o)  +  A  /  [pt(C)y(C)  -  Pn(0®(0]‘^C  , 

where  A  =  (1  +  Solving  these  equations,  taking  the  above  equations 

of  motion  into  account,  quadrupole  matrix  elements  (as  given  in  ref.  ) 

can  be  determined  for  relating  pi{z)  at  the  focal  plane  to  pjfe(O)  at  the  target. 
This  was  done  with  a  Monte  Carlo  calculation  where  a  realistic  population  of 
the  aperture  was  used  so  that  average  matrix  elements  could  be  obtained.  For 
hydrogen,  the  kinematic  correlation  between  momentum  and  scattering  angle 
was  enforced  in  the  calculation. 
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Kinematics 

Apt/pt 

^Pn/ Pn 

Api/pi 

Q2 

0.065 

2.26 

0.123 

Q3 

0.075 

2.61 

0.096 

Table  4.3.  Change  in  target  polarizations  after  taking  quadrupole  effects  into 
account. 


For  our  spectrometer  setup  of  two  quadrupoles  and  a  dipole  we  have: 

P/p  ~  MDMQaMQiPijf  .  (4.77) 

We  measured  four  things  at  the  focal  plane:  px  and  py  for  both  helicities. 
Writing  these  in  terms  of  sums  and  differences  (divided  by  2),  the  measured 
quantities  are  related  to  the  target  polarization  vector  via: 


0 

<  Qi2  > 

0  \ 

/Sy\ 

Ap 

<  Qii  > 

0 

0 

<  -Q22COSX 

<  Qi3  > 

0 

Jftst 

>  M 
H  H 

<  -Q21  cosx 
V-Qsisinx  > 

-<^33  sin  X  > 

0 

<  -Q23COSX 
-Qissinx  >  / 

Fn 

Ui"'  / 

(4.78) 

where  the  quadrupole  matrix  elements  Qij  represent  the  combined  effects  of 
both  quadrupoles.  Using  the  pseudo-inverse 

the  least  squares  estimates  of  the  target  polarization  are  obtained  as 


(4.79) 


Table  4.3  lists  the  increases  in  our  target  polarization  values  once  we  took  the 
quadrupole  effects  into  account. 


(p\“\ 

p‘r 

= 

/Sy\ 

Ay 

Sr 

W,’') 

\  Az/ 
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Results 


Table  5.1  lists  the  experimental  kinematics.  Measurements  were  taken  at 
two  different  values  (0.38  and  0.50  GeV^).  The  initial  beam  energy  was 
approximately  580  MeV;  or  579.05  MeV  at  the  hydrogen  target  center.  The 
momenta  in  the  table  include  the  effects  of  target  losses. 

5.1  Systematic  Errors 

There  are  three  contributors  to  the  systematic  errors  of  our  experiment:  the 
beam  helicity  uncertainty,  the  analyzing  power  uncertainty,  and  the  uncertainty 
in  the  depolarization  effects  of  the  spectrometer  magnets.  (The  last  one  includes 
kinematic  uncertainties.)  The  helicity  uncertainty  was  4%  and  contributed 
to  the  error  of  p/  and  p*.  Table  5.2  gives  the  individual  contributions  to 
this  systematic  error.  i4j,(fl)’s  absolute  error  of  ±0.02  gives  a  relative  error  of 
about  0.04  that  affects  all  of  the  polarizations.  Tables  5.3  and  5.4  detail  the 
systematic  error  from  the  uncertainty  in  the  magnetic  depolarization  effects  for 
=  0.38  GeV^  and  =  0.50  GeV^,  respectively.  For  p/  the  uncertainty  in 
beam’s  vertical  position  at  the  target  dominated  the  error  while  p* ’s  error  was 
dominated  by  the  uncertainties  in  <j>igt  and  the  beam’s  horizontal  position  at 
the  target.  The  systematic  uncertainties  for  p„  were  minimal.  Note  that  the 
magnet  depolarization  effect  systematic  errors  associated  with  pt  and  p/  are  not 
correlated.  This,  means  that  both  must  be  included  in  the  systematic  error  of 

K.  Joo,  private  communication. 
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0.38 

0.50 

82.72 

113.87 

376.3 

310.2 

PM  EPS  (MeV/c) 

375.3 

309.2 

(deg) 

35.22 

22.03 

202.8 

268.9 

645.9 

756.9 

Carbon  thick,  (cm) 

7 

9.5 

174.0 

237.0 

1.6 

1.0 

30.9 

28.2 

Total  good  counts 

1935342 

754067 

X 

195.87 

207.47 

Table  S.l  Experiment  Kinematics 


Ge/Gm’  Tables  5.5  and  5.6  summarize  the  systematic  errors  for  the  two 
measurements. 

The  helicity  uncertainty  is  common  to  both  measurements  and 
is  therefore  an  overall  systematic  error.  Likewise,  the  analyzing  power 
uncertainties  axe  common  to  both  pt  and  pi,  though  they  are  different  for  the 
separate  measurements,  as  those  were  done  with  different  proton  kinetic 
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Description 

Beam  position  fluctuation 

2.5 

Signal  to  background  ratio  uncertainty 

1.5 

Target  polarization  uncertainty 

1.3 

Target  angle  uncertainty 

1.5 

Systematic  diff.  between  2  detectors 

1.8 

Helicity  correlated  uncertainty 

0.2 

Intra-atomic  motion  of  bound  electrons 

0.3 

Total  Systematic  Error 

4.0 

Table  5.2  Miller  Systematic  Errors 


energies  (Tp)  and  different  carbon  analyzer  thicknesses.  The  uncertcdnties  in 
the  magnet  depolarization  effects  are  unique  to  each  measured  polarization  for 
both  measurements. 

5.2  Results 


Figure  5.1  displays  our  experiment’s  measured  polarizations  versus  their 
theoretical  values  from  the  plane  wave  impulse  approximation  (PWIA)  using 
a  dipole  fit.  Since  p*  and  p/  scale  with  helicity  the  polarization  transfer 
coefficients  Djj(=  Pjfh)  are  given  instead  to  provide  a  basis  of  comparison 
with  other  experiments.  The  inner  error  bars  are  statistical  while  the  outer 
are  statistical  with  the  systematic  magnet  depolarization  effects  added  in 
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Description 

beam  energy:  d=  1  MeV 

0 

0.0049 

beam  horiz.  tgt.  pos.:  ±1  mm 

0.0070 

0.0025 

beam  vert.  tgt.  pos.:  ±1  mm 

0 

0.0250 

^tgt  :  ±0.9  mrad 

0.0066 

0.0043 

quad.  1  to  quad.  2  vert,  align.:  ±  0.35  mm 

0 

0.0012 

quad.  1  to  quad.  2  horiz.  align.:  ±  0.61  mm 

0.0021 

0.0015 

quad,  to  dipole  vert,  align.:  ±  0.32  mm 

0 

0.0010 

Total  Systematic  Error 

0.010 

0.026 

Table  5.3  Magnet  Depolarization  Systematic  Errors  for  =  0.38  GeV® 


quadrature.  In  addition,  there  are  systematic  uncertainties  firoin  the  beam 
helicity  and  analyzing  power  which,  as  they  are  common  overall  uncertainties, 
are  not  included  in  the  displayed  error  bars.  The  size  of  the  overall  systematic 
error  for  each  polaxization  is  noted  at  the  bottom  of  figure.  AU  of  the  measured 
polarizations  are  near  their  theoretical  value,  though  Du  for  Q*  =  0.38  GeV^ 
is  a  little  low.  Tables  5.7  and  5.8  list  the  measured  results  for  =  0.38  GeV^ 
and  =  0.50  GeV^,  respectively. 

Also  included  in  the  figure  axe  preliminary  polarizations  measured  from 
deuterium  during  runs  interspersed  with  the  hydrogen  ones  during  this 
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Description 

beam  energy:  ±  1  MeV 

0 

0.0020 

beam  horiz.  tgt.  pos.:  ±1  mm 

0.0124 

0.0030 

beam  vert.  tgt.  pos.:  ±1  mm 

0 

0.0143 

•  ±0-9  mrad 

0.0127 

0.0030 

quad.  2  to  quad.  1  vert,  align.:  ±  0.35  mm 

0 

0.0006 

quad.  2  to  quad.  1  horiz.  align.:  ±  0.61  mm 

0.0035 

0.0010 

quad,  to  dipole  vert,  align.:  ±  0.32  mm 

0 

0.0006 

Total  Systematic  Error 

0.038 

0.022 

Table  5.4  Magnet  Depolarization  Systematic  Errors  for  =  0.50  GeV^ 


Api/pi 

^Pt/Pt 

Ap„/p„ 

^Grai/ Grai 

Ah/h 

0.04 

0.04 

— 

— 

A  Ay /Ay 

0.04 

0.04 

0.04 

— 

Magnet  depol.  effects 

0.038 

0.024 

0.0022 

0.044 

Table  5.5  =  0.38  GeV^  Systematic  Errors 
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Api/pi 

Apt/pt 

Ap„/p„ 

^Grat/ Grat 

Ah/h 

0.04 

0.04 

— 

— 

AAyfAy 

0.04 

0.04 

0.04 

— 

Magnet  depol.  effects 

0.022 

0.038 

0.0033 

0.044 

Table  5.6  =  0.50  GeV®  Systematic  Errors 


Q2(GeV2) 

0.38 

vT 

-0.0034±.0020±.00001±.0001 

-0.468±.006±.005±.027 

dT 

0.609±.023±.016±.035 

pGe/Gm 

0.941±.038±.027 

Table  5.7.  Experiment  Results  for  =  0.38  GeV^.  The  first  errors 
are  statistical.  The  second  errors  are  systematic  uncertainties  due  to  magnet 
depolarization  effects  and  the  third  are  correlated  overall  uncertainties  from  the  beam 
belicity  emd  Ay. 


experiment,  Their  consistency  with  the  hydrogen  residts  shows  that 

polarization  transfer  coefficients  on  the  deuteron  are  not  greatly  a^ected  by 
nuclear  binding  and  interaction  effects.  This  iitiplies  that  polarization  transfer 


J.  McIntyre,  private  communication. 
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Q2(GeV2) 

0.50 

pr 

0.0003±.0034±.000001±.00001 

D|f 

-0.434±.011±.008±.026 

Df 

0.882±.023±.013±.050 

(iGe/Gm 

0.996±.036±.023 

Table  6.8.  Experiment  Results  for  Q®  =  0.50  GeV®.  The  first  errors 
are  statistical.  The  second  errors  are  systematic  uncertainties  due  to  magnet 
depolarization  effects  and  the  third  are  correlated  overall  uncertainties  from  the  beam 
belicity  and  Ay. 


measurements  hold  great  promise  for  determining  G^. 

Finally,  Figure  5.2  displays  the  t^pG^fG^  ratios  that  we  obtained  (dark 
diamonds),  along  with  previous  data  from  Rosenbluth  separation  experiments. 
Again,  the  inner  error  bars  are  statistical  while  the  outer  ones  are  statistical 
and  systematic  added  in  quadrature.  Unlike  the  polarizations,  the  ratios  axe 
independent  of  uncertainties  in  the  beam  helidty  and  Ay.  The  fits  are  described 
in  section  1.5.  The  previous  data  are  from:  Baxtel  et  al.  (open  squares), 
Janssens  et  al.  (Xs),  and  Holder  et  al.  (open  circles).  The  dark  circles 
are  our  preliminary  data  from  the  deuteron  (offset  for  clarity). 

5.3  Conclusion  and  Recommendations 


This  experiment  has  successfully  measured  polarization  transfers  in  the 
^H(e,e'p)  reaction.  While  the  transverse  polarization  transfer  has  previously 
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Figure  5.1  Measured  polarizations.  The  inner  error  bars  are  statistical  while 
the  outer  error  bars  are  statistical  errors  added  in  quadrature  with  the  uncertainties 
due  to  magnet  depolarization  effects.  The  overall  uncertainties  from  the  beam  helicity 
and  analyzing  power  are  not  shown. 

been  measured  in  this  manner  at  Mainz,  this  is  the  first  experiment  to 
measure  all  three  polarization  components  and  determine  the  ratio 
The  ratio  measured  is  consistent  with  the  dipole  model  and  has  error  bars 
comparable  to  previous  measurements  using  Rosenbluth  separation.  Thus  the 

[3]  D.  Eyl  et  al,  Z.  Phys.  A352,  211  (1995). 
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Figure  5.2  ppG^/G’^  vs.  Q^.  The  symbols  and  Sts  are  listed  in  the  text. 

groundwork  has  been  laid  for  an  alternative  method  of  determining  form  factors, 
particularly  at  higher  Q^-values  where  Ge  can  be  difficult  to  extract  using 
Rosenbluth  separation. 

As  a  proof-of-principle  experiment,  this  experiment,  though  successful,  has 
revealed  unexpected  problems  that  hopefully  can  be  avoided  in  future  similar 
experiments.  For  example,  I  believe  there  is  room  for  a  significant  improvement 
in  the  data  rate  since  nearly  half  of  our  acquired  data  had  to  be  discarded  due 
to  irreducible  multi-hits.  Given  the  large  number  of  wire  chambers  needed 
for  this  experiment,  shielding  is  probably  more  important  than  in  the  typical 
(e,e^p)  experiment.  The  use  of  redundant  chambers  (i.e.,  more  than  two)  on 
both  sides  of  the  analyzer  would  also  help,  though  it  might  necessitate  a  more 
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creative  small-angle  rejection  system. 

I  would  also  recommend  that  future  polaximeters  be  built  to  malce 
alignment  easier  than  ours  and  that  experimenters  treat  its  determination  with 
a  higher  priority  than  we  did.  If  we  had  had  to  replace  a  wire  chamber  later  in 
our  experiment,  we  could  have  been  in  trouble  concerning  our  alignment  offsets, 
as  we  did  not  do  a  “straight-through”  run  until  well  after  the  experiment  was 
over.  Such  a  run  should  be  the  first  thing  done.  Then,  with  the  chambers  and 
supports  designed  to  allow  known  adjustments  to  be  made  easily,  one  could 
ensure  that  the  small  angle  rejection  system  did  not  bias  the  data  nearly  as 
much  as  ours,  thus  raising  the  useful  percentage  of  data  acquired. 

Finally,  this  experiment  shows  that  spin  transfer  physics  requires  well 
understood  magnetic  fields  and  magnet  alignments.  This  experiment  was 
done  on  a  spectrometer  with  unmapped  fields.  Knowledge  of  the  actual 
fields  as  a  function  of  position  would  have  allowed  us  to  better  determine 
the  spin  transport  effects  in  the  magnets,  as  well  as  the  focal  plane  to  target 
coordinate  transformations.  Better  knowledge  of  the  relative  alignments  of  the 
spectrometer  magnets  themselves  and,  more  importantly,  better  stability  and 
knowledge  of  the  beam  position  at  the  target  would  have  also  reduced  systematic 


errors. 


Appendix  A 


Derivation  of  Cross  Section  Relationships 


This  appendix  provides  detailed  derivations  of  differential  cross  section 
and  polarization  equations  used  in  this  dissertation  for  scattering  electrons 
off  protons.  The  first  section  details  the  necessary  formalism  to  start  from 
transition  currents  and  the  Dirac  equation  and  arrive  at  a  differential  cross 
section  in  terms  of  an  invariant  scattering  amplitude.  The  next  section 
introduces  a  reference  frame  -  the  Breit  frame  -  which  wiU  be  used  in  evaluating 
the  scattering  amplitude.  The  third  section  describes  the  proton  transition 
current  in  this  frame,  while  in  the  last  two  sections  the  cross  section  relations  for 
both  unpolEuized  electrons,  giving  the  Rosenbluth  equation,  and  for  polarized 
electrons,  are  derived. 

A.l  Cross  Section  Formalism 


In  this  section  the  formalism  leading  up  to  the  electron-nucleon  scattering 
cross  section  is  developed  as  in  Halzen  and  Martin  and  Chinitz.  The 

basic  ingredient  is  the  current-current  interaction  as  depicted  in  figure  A.l  in 
which  the  electron’s  transition  current  interacts  vrith  the  proton’s  transition 
current  J^*.  The  first  step  in  calculating  the  cross  section  is  to  determine  just 
what  these  currents  are. 

[57]  Halzen  and  A.D.  Martin,  Quarks  and  Leptons,  John  Wiley  &  Sons,  (1984). 

L.M.  Chinitz,  Ph.D.  dissertation,  U.Va.  (1990),  unpublished. 
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Figure  A.l  Transitions  currents  in  electron-nudeon  scattenng. 
The  relativistic  Dirac  equation  for  spin-|  particles  is 

Hip  =  (at  •  P  +  Pm)ij)  , 
which  can  be  written  covariantly  as 

(t**?/!  -  =  0  • 

For  free  electrons,  the  wave  function 

^  =  u(p)e->-  , 


(Al) 


(A.2) 


(A.3) 


where  ti(p)  is  a  four-component  spinor,  satisfies  this  equation.  Using  (A.2) 
and  the  continuity  equation,  dfij^  =  0,  gives  the  electron’s  transition  current 
between  its  initial  and  final  states  as: 


jfi  =  -eipf'i^^ipi 

=  -cu/7'*u,-e‘(*’^ . 


(A.4) 


Since  the  nucleon  is  not  a  point-particle,  cannot  be  written  in  the 
same  way  as  The  vertex  factor  7^  must  be  replaced  by  a  term  r^(p/,p,) 
that  contains  the  nucleon’s  structure  information.  This  malfps  the  nucleon’s 
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transition  current: 

.  (A5) 

If  one  were  at  the  quark  level,  however,  a  point-particle  vertex  could  stiU  be 
used.  Then  you  would  have 

=  ,  (^-6) 

i 

where  ej  represents  the  fractional  charge  for  each  quark  and  qj  represents  the 
quark  wave  functions.  Instead,  however,  must  assume  the  most 

general  four-vector  form  that  can  be  constructed  with  p/,  pi,  g(=  p/  —  p,), 
and  the  7  matrices  without  violating  parity  conservation.  This  is: 


r^(p/,Pi)  =  7^,Kl  +  i<r^vq’'K2  -I-  io-fiviPf  +  PifKz  +  (p/  -|-  pi)^K5  . 

(^.7) 

The  Gordon  decomposition, 


(A8) 


shows  that  the  {pf  +  pi)y,  terms  aie  really  just  linear  combinations  of  the  7^ 
and  a/ivQ''  terms,  so  that 


^^iPf,Pi)=7,^Fi  + —(r^^q’'F2  +  q^Fa  ,  (A9) 

where  M  is  the  nucleon  mass.  Current  conservation  at  the  vertex  gives  that 

z/c 

q^J^  =  eu/(g'‘7^f’i  +  ■^q'^cr^uq’' F2  -|-  g^F3)u,-  =  0  .  (.4.10) 

The  first  term  in  this  equation  is  zero  by  the  Dirac  equation.  The  second  is  zero 
since  it  is  a  contraction  of  symmetric  and  antisymmetric  tensors.  This  means 
that  jPs  must  be  zero  as  well.  Thus  one  is  left  with 


r„{p/,pi)  = -FiW'hM  +  •  U-ii) 

This  equation  is  valid  for  both  protons  and  neutrons,  though  the  form  factors  for 
each  are  different.  Fi  is  called  the  Dirac  form  factor  and  F2  is  called  the  Pauli 
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form  factor,  k  is  the  anomalous  magnetic  moment.  It  is  related  to  the  magnetic 
moment  fi  by  fiN  =  +  knF^{0).  For  protons  f’f(O)  =  Ff  (0)  =  ^ 

fip  =  2.79,  whereas  for  neutrons  F’"(0)  =  0,F2”(0)  =  1  and  /i„  =  —1.91.  Hence 
in  the  reed  photon  limit  »  0),  one  basically  probes  a  paxticle  of  magnetic 
moment  (1  +  k)^  and  charge  e  (for  the  proton,  0  for  the  neutron). 

The  electron  scattering  off  a  nucleon  in  figure  A.l  is  not  a  free  electron 
but  rather  is  traveling  in  an  electromagnetic  field  A^  caused  by  In  classical 
electrodynamics  the  motion  of  an  electron  through  an  electromagnetic  potential 
brings  about  the  substitution 

Pfi  ^  Pfi  eAp  .  (A. 12) 

(As  it  turns  out,  requiring  local  gauge  invariance  on  the  electron  Lagrangian 
mandates  this  substitution.)  This  changes  the  Dirac  equation  (A.2)  to 

=  -^'I'^Ap  .  (A.13) 

The  field’s  potential  is  defined  as 

/F  =  -e7'‘A^  .  (A.14) 

By  first  order  perturbation  theory,  the  transition  amplitude  of  a  particle  from 
an  initial  state  to  a  final  state  in  a  potential  is 

Tfi  =  y  V'/(®)^(®)V’i(®)<^®  •  (A.15) 

Using  jf  -  (A.4),  where  3®d  (A.14)  transforms  this  to 

Tfi  =  -i  j  j^Apct^x  .  (A.16) 

The  relationship  between  and  J^is  given  by  Maxwell’s  equations,  which  are 
written  covariantly  in  position-space  as 


(A.n) 
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Fourier  transforming  this  leads  to  the  result  that: 


Mi)  = 


This  makes  the  transition  amplitude 


=»  J 


Upon  integration,  Tfi  can  be  written  as 

____  -i(27r)4  4.- 


S^{k'  +  p'  —  k  —  p)M. 


(i4.18) 


(A.19) 


(i4.20) 


where  V  is  the  normalization  volume  and  M.  is  the  invariant  scattering 
amplitude  and  equals 


M  =  -e^li((!/)7'‘u(i,)^u(7/)r„u(p,) 


(X.21) 


For  the  scattering  process  k  +  p  — >  Aj'  +  p',  the  transition  rate  per  unit 
volume  is  defined  as 


Wr-^^ 
yyf*-  j,y 


(A.22) 


where  T  is  the  interaction’s  time  interval.  This  is  related  to  the  cross  section 


(j  =:  ^  — T (number  of  final  states) 

(imtial  flux) 


From  phase  space, 


number  of  final  states  =  . 

(27r)32E'  (27r)32£;; 

The  initial  flux  in  the  lab  frame  is  expressed  as 

initial  flux  =  I  Ve  I  . 

For  m  <<  E,  Ve  =  is  w  1.  Thus  the  lab  frame  cross  section  is 


da  = 


(2'7r)^|A^P 

AEM 


S^{k-k' +p-p') 


,,  d^k'  d^p' 


(2Ty2E'  {2Try2E'  ‘ 


(.4.23) 


(.4.24) 


(.4.25) 


(A26) 
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Note  that  \M.\^  has  been  used  rather  than  \M.^  since  an  unpolaxized  cross 
section  needs  to  be  averaged  over  the  initial  states  and  summed  over  the  final 
states.  The  energy  equation  E'^  =  gives  the  relationship  E'dE'  =  k'dk' 

making 


d^k'  =  k'^dk'dQ  =  E'^dE'dQ 


For  the  lab  &ame  one  can  also  write 


f  ^6(p  +  q-p')  =  —6((v  +  i) 

J  2E'  2M  ^  2M’ 


2MA 


6{E'  -  E/A) 


where 


and  the  relationships 


^  =  ¥7  =  1  +  17®“  (l) 


u)  =  E-E'  =  ^ 

2M 


(^.27) 


(i4.28) 


(i4.29) 


(i4.30) 


Pd  2k- k'  =  4EE'sm^{-) 

A 


(^.31) 


have  been  used.  Putting  in  these  values  makes  the  laboratory  frame  cross 
section  equation 

^  =  U32) 

dn  647r2M2  \  e)  '  ^  ^  ^ 


(^.32) 


It  is  usefiil  to  write  |Atp  in  terms  of  two  tensors  representing  the  two 
interacting  particles: 


(^.33) 


=  5  E  Kp')r»<iWI[<i(p')ri.>i(i>)] 


(A.34) 
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for  the  electron  and  proton,  respectively.  This  mahes 

W  .  (A35) 

To  evaluate  these  tensors  note  that  the  terms  in  brackets  are  1x1  matrix 
elements.  The  consequence  of  this  is  that  the  terms’s  complex  and  hermitian 
conjugates  are  the  same.  Using  =  7°7''  gives 


[u(fc')7‘'u(i)]*  =  [u(fe)7''u(fc')]  . 

By  writing  out  now  in  terms  of  explicit  matrix  elements  and  using  the 
completeness  relation 


u(*)(p)u('^(p)  =  +  TO 

«=1,2 

(where  electron  tensor  becomes: 

i<‘-  =  lTr{(jS'+TOh»(jS  +  m)7-)  . 

Using  trace  theorems  and  neglecting  the  mass-squared  term  gives 

1'“' =  2(jfe"‘ib'' -b  A!'*'*"  -  (jb  .  ib  V") 


=  2(ib"*fc‘'  -h  A:"'ib'*)  -I-  , 


(A36) 

(^.37) 

(A38) 


where  (A.31)  was  used.  To  evaluate  the  proton  tensor,  let’s  define  a  slightly 
different  version  of  the  nucleon  transition  current  (A.5)  without  the  electron 
charge  and  the  exponent 


=  u/r^(p/,pi)uf  .  (A.39) 

This  allows  the  proton  tensor  to  be  written  as 

=  ^Tr(  Jt)  .  (^40) 

AU  that  is  left  is  to  further  evaluate  the  transition  currents  and  put  in  the 
appropriate  momentum  four- vector  values.  Since  |Af|^  is  an  invariant  this  can 


J.D.  Bjorken  and  S.D.  Drell,  Relativistic  Quantum  Mechanics,  McGraw-Hill,  (1964). 
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Figure  A.2  Breit  reference  frame. 

be  done  in  whatever  frame  one  chooses.  The  reference  frame  that  we  will  use 
is  the  Breit  frame,  which  is  described  in  the  next  section. 

A.2  The  Breit  Reference  Frame 

A  useful  reference  frame  for  the  evaluation  of  transition  currents  is  the 
Breit  frame.  This  is  also  called  the  brick  wall  frame  since  in  this  frame  the 
proton  and  electron  are  essentially  “bouncing”  off  a  wall,  as  can  be  seen  in 
figure  A.2.  In  this  frame  p  =  — p'  and  There  is  no  energy  transfer 

(a;  =  0)  so 

=  l9Bp  • 

Since  qs  =  2p  this  means  that 

9|  =  -4|p1"  .  (AA2) 


If  the  z-axis  is  chosen  along  gs  with  scattering  in  the  x-z  plane  the  electron 


R.  Hagedorn,  Relativistic  Kinematics,  W.A.  Benjamin,  (1973). 
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four  momenta  are: 


*5  =  -^(“c 

*' = -^(“cyjcot  • 


{AAZ) 


Since  the  Lorentz  transformation  from  the  laboratory  frame  to  the  Breit  frame 
is  entirely  along  the  z-axis,  the  fact  that  x  is  the  same  in  both  frames  can  be 
used  to  get  the  relationship: 


X  =  ZCOt(^)  =2BCpt(^) 


(.4.44) 


cot  (y)  =  ^cot  (-)  . 

z  and  Zb  are  related  by  the  Lorentz  transformation  as: 


z^  =  {l-l3^)z\ 


where  =  u/c, 


(.4.45) 


(.4.46) 


All  that  is  left  is  to  establish  a  relationship  between  0  and  9b  is  to  determine 
the  value  of  This  can  be  done  by  looking  at  the  Lorentz  transformation  from 
the  lab  to  the  Breit  frame  and  using  the  fact  that  in  the  Breit  frame  p  =  —p'. 
The  four  momenta  in  the  lab  frame  are  p,-  =  (M,  0)  and  pf  =  (M  +  a;,^), 
whereas  in  the  Breit  frame  they  are  p,B  =  {EB,PiB)  and  pfB  =  (£?b,P/b)- 
Using  the  Lorentz  transformations  gives 

a:) =(-;?)(:)  • 


(,AA7) 


which  leads  to 


PiBz  =  -fijM  ; 


(A.48) 


N.  Dombey,  Rev.  Mod.  Phys.  41,  236  (1969). 


A.I.  Akhiezer  and  M.P.  R,ekalo,  Fiz.  Elem.  Chast.  Atom.  Yad.  4,  662  (1973)[Sov.  J. 
Part.  Nucl.  4,  277  (1974)]. 
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(i4.55) 
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A.3  Proton  Transition  Current  in  the  Breit  Frame 


Putting  (A.ll)  into  the  modified  proton  transition  current  (A.39)  gives 
i7/i  =  u_^(Fi'y/i  +  )'iii  .  (A.56) 

Using  Gordon  decomposition  (A.8)  this  can  be  written  as 

=  u/[7„(F,  +  kFz)  -  •  (A.67) 

To  determine  the  different  there  are  three  quantities  that  need  to  be 
evaluated:  u/u,-;  u/7oUi;  and  u/7fcU,-,  where  k  =  1,2,3.  Note  that  in  the 
Breit  frame,  where  p  =  —p\  the  second  term  in  (A.57)  is  zero  for  ,7*. 

For  positive  energy,  the  four-component  spinor  ‘tt(p)  is: 

=  (^68) 

' E+m^  ' 

while  the  corresponding  ii(p')  is: 

^  •  (A.59) 

\  +  TO  J 

5  =  1,2  and  x  is  a  two- component  spinor  with 

XiXj  =  Sij  •  (A.60) 

Using  these  values  in  the  Breit  frame  gives: 

u,u,  =  (E  +  M)(l  +  ^^^) 

_E^  +  M^  +  2EM  +  1^2  (A.61) 

E  +  M 

=  2E 


where  the  relationship  (or  •p)^  =  has  been  used.  Similarly,  u/7oUi  =  2M. 
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Using  these  results  makes  the  modified  charge  density: 

Tp2 

Jo  =  2M[Fi+(l-— )/ci?’2] 

=  2M[Fi  -  (^-62) 

=  +  ^kF,] 

where  (A.42)  has  been  used  in  the  last  step.  By  defining 

Ge  =  F.  +  (^.63) 

one  gets 

Jo  =  2MGe  .  (AM) 

The  remaining  term  to  evaluate  is  u/7;iUi.  Using  the  four-component 
spinors  (A.58  and  A.59)  here  one  gets 

u/7jfeUi  =  •  P)  +  •  p')^k]x^‘^  •  (A.65) 

Using  the  relationship 

(a  •  a)(d^  •  6)  =  a  •  6  -f  •  (o  X  6)  (A.66) 

gives 

=  X^*  •  (p  X  fifc)  +  p'k  +  •  {p'  X 

=  X^''^^*^[(P '  -  P )  X  (A.67) 

=  X  q)kX^‘'>  . 

This  makes  the  modified  proton  transition  current 

=  .  (A68) 

J.J.  Sakurai,  Advanced  Quantum  Mechanics,  Addison- Wesley,  (1967). 

[64]  Y  3^  Berestetskii,  E.M.  Lifshitz,  and  L.P.  Pitaevskii,  Relativistic  Quantum  Theory, 
Addison- Wesley,  (1971), 
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As  in  the  case  of  it  is  convenient  to  use  a  linear  combination  of  and  F2] 
in  this  case, 

Gm  =  “I*  /ci'2  •  (A..69) 


This  gives 

Jk  =  X  q)kGMX^''^ 


(A.70) 


From  this  current  and  (A.64)  it  can  be  shown  that  in  the  Breit  frame  Ge  and 
Gm  (often  called  the  “Sachs  form  factors”)  are  the  Fourier  transforms  of  the 
charge  and  magnetic  moment  distributions  of  the  proton. 


By  defining  the  modified  proton  transition  current  as 


[61]  [62] 


(A.71) 

one  can  clearly  write  the  proton  tensor  (A.40)  as 

(A.72) 

using  (A.60).  The  Fe  are  then  simply: 

To  =  2MGe 

(A.73) 

and 

Fk  =  i{c  X  q)kGM  . 

(A.74) 

More  specifically,  since  qs  lies  entirely  along  the  z-axis  in  this  frame,  the 
components  of  J-k  are: 

=  i\qB\GM<^2 

F2  =  —i\qB\GMO'i  (A.75) 

:F3  =  0  . 


The  formalism  necessary  to  evaluate  the  cross  sections  is  now  in  place. 
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A.4  Rosenbluth  Formula 


Putting  the  above  Ts  into  the  proton  tensor  (A.72),  the  proton  tensor 
components  are: 

Woo  =  AM^Gl 

(A.76) 

Wii  =  W22  =  Q^Gh  , 


with  all  of  the  others  equal  to  zero.  Likewise,  using  the  Breit  frame  four 
momenta  (A.43)  in  the  electron  tensor  (A.38),  gives  the  nonzero  electron  tensor 
components: 


L^°  =  Q^icsc^^-l)  =  Q^coe^ 

£i  £d 


=  Q‘‘  cot  ^  CSC  ^ 

2  2 

£"  =  QVt’'y  +  l)  =  i?’'csc^Y 


Ob  Ob 


{A.77) 


(Please  note  that  the  electron  tensor  is  symmetric  in  this  case.)  Therefore  the 
electron-proton  tensor  product  is 


cot^  I 


L'^’'W,.  =  AM^G%Q^^^  +  Q^GI{2  +  —^)  ,  (A.78) 

1  T  1  +  T 

where  (A. 55)  has  been  used  to  convert  from  the  Breit  frame  to  the  lab 
frame.  Putting  this  into  the  cross  section  equation  (A.32  and  A.35)  gives  the 
Rosenbluth  cross  section 


da 

dn  ^  647r2M2g4 


(A.79) 

It  can  be  insightful  to  write  this  equation  with  the  nonstructure  cross  section, 


a^cos^(0e/2)  .E' 

^Wsm\0,/2yW  ’ 


(A.80) 
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separated  out.  Using  (A.29-31),  it  is  straightforward  to  do  this  and  get 

d(T  1  _  cos^(5e/2) 

I  "  4E2siii‘(l!l./2)^¥' 

A.5  Polarizations 

In  this  section  polarizations  induced  in  protons  scattered  by  polarized 
electrons  will  be  derived,  as  is  discussed  in  Akhiezer  and  Rekalo.  These 
polarizations  are  in  terms  of  the  Sachs  form  factors  and  allow  for  the  the  ratio 
of  Ge  to  Gm  to  he  determined.  For  derivation  of  polarizations  using  a  response 
function  formalism  please  see  Donnelly  and  Raskin. 

For  a  polarized  electron  in  the  extreme  relativistic  limit,  the  term  |(l+h7®) 
must  be  introduced  into  the  electron  tensor  (A.37)  to  project  the  different 
electron  spin  states,  h  represents  the  electron  beam’s  helicity,  which  is  ±1 
for  a  purely  polzirized  beam  and  |h.|  <  1  for  a  partially  polarized  beam,  as  in 
this  experiment.  This  adds  an  antisymmetric  part  to  the  previously  discussed 
symmetric  part  of  the  electron  tensor 

This  antisymmetric  part  is  (ignoring  the  mass  terms): 

[65]  Donnelly  and  A.S.  Raskin,  Ann.  Phys.  169,  247  (1986). 

R.G.  Arnold,  C.E.  Carlson,  and  F.  Gross,  Phys.  Rev.  C  23,  363  (1981). 


(A.82) 


(A.83) 


Appendix  A:  Derivation  of  Cross  Section  Relationships  152 


As  was  done  with  the  symmetric  part,  the  Breit  frame  four-momenta  (A.43) 
are  used  to  evaluate  this  tensor’s  components,  the  results  being: 


cot 


rl2  _  7-21  _  •/rj2  ''O 

Da  -  -Da  =  -^Q  CSC  —  , 


(A.84) 


with  the  other  components  being  zero. 


Of  course  polarization  also  affects  the  proton’s  tensor  (A.72).  Specifically, 
the  proton  tensor  denoting  the  final  polarization  in  the  direction  n  of  a  proton 
which  was  initially  unpolarized  is  written  as: 


(A.85) 


Evaluating  the  tensor  components  is  simply  a  matter  of  putting  in  the  !F  values 
(A.73  and  A.75).  For  example. 


=  -Q^Gl,\Tr( 


I  \  /  riz  Ux  —  iuy  \  ' 
/  \nx+iny  -Uz  /. 
Vfhx  “{“  \ 

iuz  ) 


This  gives  while  =  0  since  their  traces  are  zero.  A 

complete  listing  of  the  nonzero  proton  tensor  components  is: 

'pS  =  -^^20^  =  -Pif  =  =  2iMmGEGM  ,  . 


7?(*)  _  'p(*)  _ 

r  19  —  — r  91  — 


(A.86) 


Note  that  if  one  were  interested  in  an  experiment  using  a  polarized  target 
instead  of  measuring  the  final  polarization  of  the  protons,  the  {B  •  h)  term  in 
(A.85)  would  have  been  sandwiched  between  the  ^’s  rather  than  at  the  end. 
The  tensor  components  in  this  case  would  have  been  identical  except  for  a  sign 
difference  in  the  z-direction.  Hence  both  options  give  the  same  information. 
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Since  the  axe  antisymmetric,  and  the  contraction  of  symmetric  and 

antisymmetric  tensors  is  zero,  only  the  antisymmetric  part  of  the  electron  tensor 
is  important  for  determining  the  polarized  cross  section  Again  using 

(A.32  and  A.35)  one  gets: 


dfle 

do-W  ha^  On 


(A.87) 


dU.  2M2 '  E 


(-Er)G'MCsc—  . 


The  polarized  cross  section  divided  by  the  unpolarized  cross  section  gives  the 
polarization: 

ds^ 

=  ^  •  (^-88) 

dSle 

I  wiU  give  the  values  of  the  polarizations,  as  stated  by  Arnold  et  al.  ,  after 
going  through  some  of  the  algebra. 


The  unpolarized  cross  section  can  be  written  as: 

dcr  da  Iq 

dCle  dQ  1  +  T 


where 


Jo  =  +  rGlt  +  2r(l  +  T)Glt  tan^  1 


(A.89) 


(A.90) 


and  T  =  By  rewriting  the  nonstructure  cross  section  (A.80)  as 


d<7l  _  a2  j;'  2  tfe 


(A.91) 


and  using  (A.55),  the  Breit-lab  relationship,  it  is  simple  to  get  the  relationship: 


PxIq  —  — vt(1  +  T)tan  — 


(A.92) 
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Somewhat  more  algebra  is  needed  to  get  Arnold  et  aVs  other  polarization 
equation: 

pJo  =  2hrG'M sec  ^  tan  ^[(1  +  t)(1  +  r  sin^  ^)]  * 

By  taMng  the  ratio  of  these  two  measured  quantities,  one  obtains  the  ratio  of 
Ge/Gm  independent  of  the  beam’s  helicity: 

Px  _ 2M  Ge 

Pz  (E  +  £?')  tan  I  (rAf 

Note  that  the  coordinate  frame  in  which  this  was  evaluated  has  z  along  q  and 
X  trcinsverse  to  it  in  the  scattering  plane  (as  opposed  to  our  spectrometer’s 
coordinate  frame,  depicted  in  figure  4.4,  in  which  x  is  normal  to  the  scattering 
plane).  Thus  to  avoid  any  misunderstandings,  I  rewrite  the  previous  equation 

pt  _  2M  Ge 

Pi  {E  +  E')  tan  |  Gm  . 


as: 


(A.94) 
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FPP  Calibration  at  lUCF 


The  FPP  was  calibrated  at  lUCF  during  February  1993.  This  appendix 
gives  a  brief  description  of  the  calibration  run  and  results,  based  on  the  more 
detailed  description  by  our  group  in  Lourie  et  al. 

At  lUCF  the  FPP  was  directly  exposed  to  a  low-intensity  (fs  10®/s)  proton 
beam  in  the  QQSP  area.  This  low-intensity  was  achieved  using  many  slits, 
collimators,  and  a  lead  degrader.  The  FPP  was  mounted  on  its  side  so  that 
its  bottom  chamber  faced  the  beam-pipe.  The  table  on  which  the  FPP  was 
mounted  was  wheeled  and  had  a  hydraulic  lift  so  that  the  FPP  could  be  moved 
in  all  three  dimensions  with  respect  to  the  beam-pipe.  This  allowed  us  to 
illuminate  many  different  regions  of  the  FPP  and  at  diflFerent  angles.  The  FPP- 
table  assembly  had  to  be  craned  into  the  QQSP  area  through  its  removable 
roof. 


Beam  energies  of  120,  150,  180,  and  200  MeV  were  used  and  the  beam 
polarization  measured  by  the  FPP  was  compared  to  the  lUCF  beamline 
polarimeters.  By  reversing  the  polarization,  instrumental  asymmetries  were 
canceled  out.  The  ratio  of  the  FPP’s  measured  polarization  magnitude 
at  200  MeV  to  lUCF.’s  was:  0.976  i  0.006.  More  important  than  the 
magnitude  measurement  is  the  ratio  of  spin-up  to  spin- down  polarizations, 
where  asymmetries  are  not  canceled.  This  ratio  was  also  measured  by  an  lUCF 
low-energy  polarimeter  using  (p,a)  scattering.  The  FPP  ratio  to  lUCF’s  was: 
1.018  ±0.011.  Figure  B.l  shows  our  measured  asymmetry  compared  to  previous 


(166) 


(0)Vd 
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fifpp  (deg) 


Figure  B.l  FPP  asymmetry  vs.  scattering  angle  at  lUCF  (solid  circles ).  Other 
data  are  from  SIN  (open  diamonds)  and  LAMPF  (open  squares).  The  two  Sts  are 
&om  the  SIN  data  set. 


data  from  SIN  and  LAMPF. 
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FPP  Alignment 


The  track  reconstruction  and  caxbon  scattering  angles  determination 
described  in  section  4.5.1  require  that  the  four  FPP  wire  chambers  be  aligned 
well,  both  with  respect  to  each  other  and  to  the  VDCX.  When  the  chambers 
were  made,  marks  were  placed  on  the  outside  that  corresponded  to  partictilar 
wires  on  the  inside.  Before  beginning  our  experiment,  much  time  was  spent 
surveying  the  FPP/VDCX  setup.  This  was  all  for  naught,  however,  as  various 
problems  caused  the  different  chambers  and  the  FPP  itself  to  be  taken  out  of 
OHIPS  prior  to  the  production  data  acquisition.  Also,  about  a  tenth  of  the  way 
into  this  experiment  a  wire  broke  in  chamber  3,  requiring  the  complete  removal 
of  the  FPP  and  the  quick  replacement  of  a  chamber.  (Good  thing  we  built 
spares!)  Hence,  we  only  had  a  crude  understanding  of  the  alignment  during  the 
experiment.  It  was  not  until  our  experimental  program  ended  months  after  this 
experiment  (see  Epilogue)  that  the  FPP  (still  in  the  same  configuration  since 
the  chamber-3  replacement)  was  surveyed.  The  walls  of  the  OHIPS  hut  were 
removed  and  concrete  blocks  were  stacked  so  that  a  theodolite  could  survey  the 
FPP  while  it  was  stiU  on  OHIPS. 

The  survey  was  concerned  primarily  with  the  two  laxge  back  chambers,  so 
that  the  absolute  angular  alignment  of  the  FPP  could  be  established.  The  two 
chambers  were  found  to  be  weU-aligned  except  for  a  -1.1  mm  offset  needed  in 
chamber  3’s  y-axis.  The  survey’s  uncertainty  was  dominated  by  the  uncertainty 
in  the  fiducial  mark  location  (0.1  mm).  This  leads  to  a  systematic  uncertainty 
of  0.3  mrad  in  the  FPP’s  ability  to  measure  an  absolute  angle. 


(157) 
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Chamber 

XI 

0.151 

-0.009 

X2 

0.103 

-0.057 

X3 

0 

-0.160 

X4 

0 

-0.160 

Y1 

-0.573 

-0.419 

Y2 

-0.356 

-0.202 

Y3 

-0.111 

0.043 

Y4 

0 

0.154 

Table  C*l.  FPP  Chamber  Offsets.  The  internal  offsets  refer  to  the  amount 
the  individual  chamber  axes  bad  to  be  shifted  relative  to  chamber  4.  The  external 
offsets  refer  to  the  amount  the  individual  chamber  axes  bad  to  be  shifted  relative  to 
the  VDCX. 


Next  it  was  necessary  to  align  the  front  two  chambers  relative  to  the  back 
ones.  During  sieve-slit  runs  a  month  after  the  completion  of  our  experiment, 
the  carbon  was  taken  out  of  the  FPP  so  that  “straight-through”  tracks  could 
be  recorded.  Using  these  data,  we  compared  actual  positions  in  the  front  two 
chambers  to  those  predicted  using  the  back  two  chambers.  The  differences 
between  them  were  called  X1(2)DIF  and  Y1(2)DIF  and  were  correlated  with 
0fp  and  respectively.  (The  focal  plane  angles  were  determined  by  the 
VDCX  in  this  case.)  These  correlations  were  exploited  to  determine  the  offsets 
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Y  IND  246  NO  TITLE  IN  FILE 
YLO=  2183  YHI=  2199 


TH=  45 


30 


TILDA2  BLK  16  X  IND  245  NO  TITLE  IN  FILE 

RUN  230  17-JAN-96  SUM=  21825.  XLO=  3208  XHI=  3165  CORE 

TST  464  17:22:33  ADD=  1:  1  FAL=  1064  LST=  0 


Figure  C.l  x  =  flcos(^)  versus  y  =  0sin(^)  showing  bow  the  FPP^s 
misalignment  effects  the  small-angle  rejection  system. 

necessary  for  the  DIP  spectra  to  center  at  zero.  This  task  was  not  as  easy 
as  one  would  think  due  to  the  fact  that  the  DIF  spectra  were  not  broad  peaks, 
but  composed  of  several  discrete  peaks  due  to  the  “integerization”  of  the  wire 
chamber  positions  by  the  wires  themselves.  Table  C.l  lists  these  “internal” 
offsets.  Rotation  angles  in  the  x-y  plane  (yaw)  were  also  determined  but  did 
not  have  much  effect  and  were  ignored.  Pitch  and  roU  are  second  order  effects 
and  were  ignored. 

G.  Warren,  Bates  Technical  Note  #96-01. 
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TST  469  14:53:43  A0D=  2:  1  FAL=  0  LST=  0  29-JAN-96  13:26:34 


TST  468  14:53:33  Ar)D=  2:  1  FAL=  0  LSI-  0  29-JAN-96  13:26:34 


Figure  C,2  x  =  ^cos(^)  cut  on  — 2  <  3/  =  ^sin(^)  <  2  and  vice  versa.  The 
area  where  the  distributions  decrease  toward  0  are  regions  where  the  data  is  biased. 


All  that  was  left  was  to  align  the  FPP  relative  to  the  focal  plane  of  the 
VDCX.  Again,  this  was  done  by  taking  DIF  spectra  between  FPP  positions 
projected  to  the  focal  plane  and  what  the  VDCX,  which  is  accurate  in  position 
to  0.1  mm,  measured.  In  this  case,  “integerization”  is  not  a  problem  since 
a  drift-chamber  is  involved.  We  determined  that  the  FPP’s  x-axis  had  to  be 
shifted  by  —0.160  cm  and  the  y-axis  by  0.154  cm.  These  offsets  are  reflected  in 
the  “external”  offsets  of  table  C.l. 


The  most  problematic  aspect  of  misaligned  chambers  is  not  the  subsequent 
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Y  IND  246  NO  TITLE  IN  FILE 
YL0=  0  YHI=  0 


TH=  45  PHI=  30 


Figure  C.3  Event  distribution  after  software  6 -cut  for  the  same  data  as  fig, 
C,1  (but  at  a  different  scale).  Events  whose  polar  scattering  was  less  than  7°  have 
been  excluded,  as  have  those  above  20°,  where  the  analyzing  power  is  not  well  known. 

need  to  determine  software  offsets.  Rather,  it  is  its  effects  on  the  FPP’s  small- 
angle  rejection  system  (described  in  section  3.6.2.2).  Some  events  that  should 
have  been  rejected  are  accepted  and  others  that  should  have  been  accepted  are 
not.  The  result  of  this  is  shown  in  figure  C,l,  where  x  =  0cos{(f>)  is  plotted 
versus  y  =  0  sin(^)  (these  are  second  scattering  angles),  This  effectively  shows 
the  events  that  passed  the  MLU.  The  “dead”  region  in  the  center  corresponds 
to  rejected  events.  Instead  of  being  centered  on  the  small-0  peak,  it  is  obviously 
off  to  one  side. 
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The  above  observation  necessitated  finding  the  smallest  Value  of  6  that 
produced  an  unbiased  sample.  Notice  that  even  if  the  FPP  had  been  perfectly 
aligned,  we  still  would  have  had  to  reject  events  that  fell  within  a  circle  whose 
radius  was  the  distance  from  the  center  of  fig.  C.l’s  dead  square  to  its  corners. 
This  is  because  the  hardware  electronics  implemented  a  box  cut  whereas  we 
require  a  circular  one  to  avoid  azimuthal  biasing.  For  the  early  part  of  our  data 
taking,  the  MLU  was  set  to  reject  events  with  0  below  5°.  Part  way  through 
the  experiment  (and  prior  to  the  figure  shown),  it  was  decided  that  3.5°  was  a 
more  optimal  setting.  These  values  correspond  to  radii  of  rejection  of  7.1°  and 
5°,  respectively.  To  determine  what  we  actually  needed  to  reject,  we  plotted 
x{=6  cos(^))  cut  on  —2  <  y{=  0  sin(^))  <  2  and  vice  versa.  These  are  shown  in 
figure  C.2.  Those  areas  where  the  distributions  of  events  are  decreasing  toward 
zero  are  where  bias  exists  in  the  data  due  to  both  the  chamber  misalignment 
and  the  MLU’s  square  “angular”  rejection.  This  figure  is  from  runs  where  the 
MLU  was  set  at  3.5°.  It  implies  that  data  below  \/5°^  +  4°^  =  6.4°  is  biased. 
We  only  analyzed  data  above  7°  (9°  for  the  early  data).  Recall  from  fig.  2.3 
that  this  is  stiU  not  in  the  region  where  Ac  pealcs.  Figure  C.3  shows  the  data 
of  fig.  C.l  with  this  software  fl-cut. 
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As  alluded  to  in  the  first  paragraph  of  chapter  one,  the  measurement 
described  in  this  dissertation  was  done  in  conjunction  with  several  others,  all 
using  the  FPP.  The  FPP  was  moved  to  Bates  in  the  spring  of  1993  following  its 
calibration  run  at  lUCF  and  installed  in  the  South  Hall  that  fall.  After  this, 
one  bad  thing  after  another  happened:  months  had  to  be  spent  resurveying  and 
relearning  the  recirculator  tune,  the  polarized  source  befell  one  disaster  after 
another,  the  target  ruptured,  etc.  (A  favorite  quote  of  the  graduate  students 
was  “One  woe  be  passed.  Lo,  here  cometh  two  more.”).  As  the  need  to  produce 
results  and  theses  became  more  and  more  desperate,  the  experiment  grew  to 
include  measurements  that  would  have  otherwise  been  done  as  part  of  separate 
experiments,  adding  students  and  collaborators  as  it  grew.  That  which  had 
been  a  big  experiment  already  became  a  humongous  experiment  which  started 
in  January  of  1995  and  did  not  end  until  early  July  of  that  year.  The  many 
measurements  resulted  in  eight  Ph.D.  dissertation  topics,  including  mine. 

Taken  concurrently  with  my  hydrogen  data  in  January  and  February  of 
1995  was  quasielastic  deuterium  data  at  the  same  two  values.  This  was  the 
dissertation  data  for  Justin  I.  McIntyre  of  William  &  Mary.  In  April  we  again 
took  quasielastic  =  0.38  GeV^  deuterium  data,  but  this  time  in  nonparaUel 
kinematics  =  19®  and  Precoii  =  100  Mev/c).  This  dissertation  data  was 
for  fellow  University  of  Virginia  graduate  student  David  H.  Barkhuff.  A  May 
deuterium  measurement  (also  at  Precoil  =  100  Mev/c)  in  the  “dip”  region 
beyond  the  quasielastic  peak  provided  data  for  Kyungseon  Joo  of  Massachusetts 
Institute  of  Technology.  In  February  gind  March,  while  the  polarized  source 
was  down,  ^^C(e,e'p)  data  was  taken  for  Rhett  J.  Woo  of  W&M  to  study  the 
induced  normal  polarization  as  a  function  of  recoil  momentum  in  order  to  test 
for  non-PWIA  nucleon-nucleus  dynamics.  Induced  normal  polarizations  were 
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also  looked  for  in  the  June  ^H(e,e'p)7r® experiment,  whose  data  comprises  the 
dissertation  of  M.I.T.’s  Glen  A.  Warren.  This  measurement  at  a>  =  1232  MeV 
was  designed  to  help  constrain  background  models  in  the  N  A  transition. 
^H(e,e'p)7r°measurements  from  May  to  early  July  of  the  response  functions  Rt 
and  Rlt  over  a  range  of  invariant  masses  were  taken  to  study  the  C2/M1  ratio 
in  the  N  A  transition.  Christophe  Mertz  of  Arizona  State  University  and 
Costas  VeUidis  of  the  University  of  Athens  analyzed  this  data. 

It  should  be  pointed  out  that  many  of  these  other  measurements  used  the 
hydrogen  results.  Since  there  are  no  final  state  interactions  for  an  outgoing 
proton  in  hydrogen,  p„  is  very  nearly  equal  to  0,  and  there  are  no  physical 
asymmetries  in  the  unpolarized  proton  data.  This  allowed  the  FPP’s  systematic 
asymmetries  to  be  determined  (as  described  in  section  4.6)  and  used  in  the 
hydrogen  and  non-hydrogen  data  alike.  Since  the  PWIA  can  be  considered 
valid  for  hydrogen,  the  measured  polarizations  of  the  other  experiments  could 
be  normalized  with  respect  to  hydrogen’s. 

Given  the  size  of  this  project,  there  are  a  lot  of  people  who  deserve  my 
thanks  and  acknowledgement.  Foremost  among  them  would  be  our  group 
leader  and  my  advisor,  Robert  Lourie.  He’s  a  brilliant  physicist  and  it  has 
been  a  pleasure  to  learn  under  him.  He  took  a  very  active  role  in  the  early 
planning  and  development  of  the  experiment’s  code  and  electronics  and  was 
more  involved  in  the  day-to-day  management  of  keeping  the  experiment  running 
successfully  than  I’ve  ever  seen  from  a  professor.  He  was  also  far-and-above  the 
best  classroom  teacher  I  had  at  UVa. 

The  other  spokesmen  for  my  experiment  were  Mike  Finn  and  Charles 
Perdrisat  of  William  &  Mary  and  Paul  Ulmer,  formerly  at  CEBAF  and  now 
at  Old  Donunion  University.  One  couldn’t  ask  for  a  nicer  group  of  experiment 
leaders.  Dr.  Finn  was  quite  involved  with  the  experimental  trigger  electronics 
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and  served  as  a  good  counter-balance  to  Robert.  Dr.  Perdrisat  provided  the 
FPP’s  scintillators  and  Dr.  Ulmer  was  quite  involved  with  determining  the 
kinematics,  beam  times,  etc.,  needed  for  making  our  experiment  a  success. 

Scott  Van  Verst,  now  at  the  Washington  Dept,  of  Health’s  Division  of 
Radiation  Protection,  was  our  group’s  UVa  postdoc  for  most  of  my  time  with 
the  group,  though  he  became  an  M.I.T.  research  scientist  right  before  the 
data  started  coming  in.  He  oversaw  the  building  of  the  FPP  and  was  chiefly 
responsible  for  the  experiment’s  complicated  electronics.  Scott  was  very  good 
at  explaining  things  in  simple  terms  and  1  enjoyed  our  many  lunches  together 
at  the  Pav. 

David  Barkhuff  was  my  friend  and  classmate  even  before  I  joined  this 
experimental  group  in  the  spring  of  1991  and  it’s  too  bad  that  we  were  not 
able  to  work  together  more,  since  he  was  at  Bates  and  I  was  primarily  in 
Charlottesville.  David  built  our  wire-chambers  pretty  much  all  by  himself,  was 
responsible  for  our  experimental  control  electronics,  and  developed  many  of  the 
key  parts  of  our  analyzer  code,  such  as  the  FPP  and  polarization-determination 
codes  -  much  more  than  any  one  student  should  have  to  do  to  get  a  Ph.D. 

Justin  McIntyre  grew  from  a  fresh,  rash-decision-making  new  student  to  a 
dependable,  take-charge  compatriot  and  also  my  best  friend  in  the  experiment. 
I  enjoyed  our  many  shifts  together  and  our  many  outings  together,  especially 
those  many  walks  at  Crane’s  Beach  and  conversations  in  the  Danvers  YMCA 
sauna.  I  think  those  more  than  anything  got  us  through  those  agonizing 
setbacks  of  1994.  Justin  was  deeply  involved  in  the  experiment’s  electronics, 
particularly  the  extensive  trigger  electronics. 

M.I.T.  was  well  represented  by  Glen  Warren,  Kyungseon  Joo,  and  Dan 
Dale  -  a  postdoc  who  later  became  a  professor  at  the  University  of  Kentucky. 
Glen  was  our  man  in  charge  of  the  target  and  also  developed  many  valuable 
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tools  for  some  of  the  not-so-fun  parts  of  our  analysis,  such  as  chamber  alignment 
and  vdcx  drift-time  to  drift-distance  conversion.  Kyungseon  analyzed  aU  of  the 
M0ller  data  for  us.  Dan  Dale  worked  extensively  on  getting  the  wire  chambers 
to  work  and  getting  our  lUCF  calibration  run  to  be  successful. 

Though  our  run  at  Bates  was  not  a  smooth  one,  it  was  certainly  no  fault 
of  the  people  there,  who  axe  top-notch  and  deserve  thanks  for  trying  so  hard  to 
make  our  experiment  a  success.  I  especially  enjoyed  working  with  Karen  Dow  on 
data  acquisition  matters.  Manousch  Faxakhondeh  and  George  Dodson  worked 
very  hard  to  get  us  polarized  beam,  as  did  the  accelerator  guys  in  general, 
particularly  Steve  Bradley,  the  tune-master.  Dan  Tieger  was  the  spectrometer- 
dude.  The  target  group  and  crane-operators  (especially  Stevie)  worked  very 
hard  for  us,  as  did  everyone  in  general.  I  appreciate  all  the  watching-out  for  us 
students  that  Anne,  Sheila,  Ginny,  and  Audrey  provided. 

Thanks  also  to  the  lUCF  staff  that  made  our  calibration  run  there  so 
successful:  Chuck  Foster,  Ed  Stephenson,  chamber-specialist  Keith  Solberg, 
computer  expert  N.  Yoder,  Jim  Grail  and  his  terrific  crew,  and  Theresa  May. 

Thanks  to  Christophe  and  Rhett  for  their  extensive  data-takdng  shift  work. 
Thanks  also  to  former  UVa  postdoc  Joe  Mitchell  (now  at  CEBAF)  for  his  work 
in  the  early  stages  and  to  Shalev  Gilad,  Bill  Bertozzi,  and  George  Sechen  of 
M.I.T.  for  their  work  in  getting  our  wire  chambers  built. 

Thanks  to  the  Cal  State-L.A.  group:  Marty  Epstein,  S.  Mukhopadhyay, 
and  D.  Liu  for  their  help  in  data  taking,  particularly  at  lUCF.  Old  Dominion 
Univ.’s  Larry  Weinstein  also  took  a  lot  of  shifts,  particularly  during  some  key 
check-out  runs.  Thanks  to  W&M  postdocs  R.  Pourang  for  his  scintillator 
work  and  Mark  Jones  for  his  data-taking  work.  W&M’s  Chris  Armstrong  also 
helped,  as  did  Pete  Markowitz  of  the  University  of  Mzuryland  (now  at  Florida 
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To  everyone  out  there  who  has  just  read  these  acknowledgements  without 
bothering  to  read  aJl  or  part  of  the  dissertation  itself,  I  say  this:  Have  you 
thought  about  getting  a  life?  This  is  not  the  most  important  part  of  the 
dissertation,  just  THE  END. 


